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I. INTRODUCTION 

A hundred years and more of fundamental research in 
atomic and nuclear physics has shown that all matter is 
corpuscular, with the atoms that comprise us, themselves 
containing a dense nuclear core. This core is composed of 
protons and neutrons, referred to collectively as nucleons, 
which arc members of a broader class of femtometrc-scale 
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particles, called hadrons. In working toward an under- 
standing of hadrons, we have discovered that they are 
complicated bound-states of quarks and gluons. These 
quarks and gluons are elementary, pointlike excitations, 
whose interactions are described by a Poincare invariant 
quantum non-Abelian gauge field theory; namely, quan- 
tum chromodynamics (QCD). The goal of hadron physics 
is the provision of a quantitative explanation of the prop- 
erties of hadrons through a solution of QCD. 

Quantum chromodynamics is the strong-interaction 
part of the Standard Model of Particle Physics and solv- 
ing QCD presents a fundamental problem that is unique 
in the history of science. Never before have we been 
confronted by a theory whose elementary excitations are 
not those degrees-of-frecdom readily accessible via exper- 
iment; i.e., whose elementary excitations are confined. 
Moreover, there are numerous reasons to believe that 
QCD generates forces which are so strong that less-than 
2% of a nucleon's mass can be attributed to the so-called 
current-quark masses that appear in QCD's Lagrangian; 
viz., forces capable of generating mass from nothing, a 
phenomenon known as dynamical chiral symmetry break- 
ing (DCSB). 

Neither confinement nor DCSB is apparent in QCD's 
Lagrangian and yet they play the dominant role in deter- 
mining the observable characteristics of real- world QCD. 
The physics of hadrons is ruled by emergent phenomena 
such as these, which can only be elucidated through the 
use of nonperturbative methods in quantum field the- 
ory. This is both the greatest novelty and the greatest 
challenge within the Standard Model. We must find es- 
sentially new ways and means to explain precisely via 
mathematics the observable content of QCD. 

The complex of Dyson-Schwingcr equations (DSEs) is 
a powerful tool, which has been employed with marked 
success to study confinement and DCSB, and their im- 
pact on hadron observables. This will be emphasised and 
exemplified in this synopsis of the KITPC lecture series, 
which describes selected progress in the study of mesons 
and baryons, and strongly coupled quarks in-medium. 
Our compilation complements and extends earlier and 
other reviews flT-fl3l]. 



II. HADRON PHYSICS 

The basic problem of hadron physics is to solve QCD. 
This inspiring goal will only be achieved through a joint 
effort from experiment and theory because it is the feed- 
back between them that leads most rapidly to improve- 
ments in understanding. The hadron physics commu- 
nity now has a range of major facilities that are ac- 
cumulating data, of unprecedented accuracy and preci- 
sion, which pose important challenges for theory. The 
opportunities for researchers in hadron physics promise 
to expand with the use of extant accelerators, and up- 
graded and new machines and detectors that will ap- 
pear on a five-to-ten-year time-scale, in China, Germany, 



Japan, Switzerland and the USA. A short list of facilities 
may readily be compiled: Beijing's electron-positron col- 
lider; in Germany - COSY (Jiilich Cooler Synchrotron), 
ELS A (Bonn Electron Stretcher and Accelerator), MAMI 
(Mainz Microtron), and FAIR (Facility for Antiproton 
and Ion Research) under construction near Darmstadt; 
in Japan - J-PARC (Japan Proton Accelerator Research 
Complex) under construction in Tokai-Mura, 150km NE 
of Tokyo, and KEK, Tsukuba; in Switzerland, the AL- 
ICE and COMPASS detectors at CERN; and in the USA, 
both the Thomas Jefferson National Accelerator Facility 
(JLab), currently being upgraded, with new generation 
experiments expected in 2016, and RHIC (Relativistic 
Heavy Ion Collider) at Brookhaven National Laboratory. 

Asymptotic coloured states have not been observed, 
but is it a cardinal fact that they cannot? No solution to 
QCD will be complete if it does not explain confinement. 
This means confinement in the real world, which contains 
quarks with light current-quark masses. This is distinct 
from the artificial universe of pure-gauge QCD without 
dynamical quarks, studies of which tend merely to focus 
on achieving an area law for a Wilson loop and hence are 
irrelevant to the question of light-quark confinement. 

In stepping toward an answer to the question of con- 
finement, it will likely be necessary to map out the long- 
range behaviour of the interaction between light-quarks; 
namely, QCD's /3-function at infrared momenta. In this 
connection it is noteworthy that the spectrum of meson 
and baryon excited states, hadron elastic and transition 
form factors, and the phase structure of dense and hot 
QCD, all contribute information that is critical to elu- 
cidating the long-range interaction between light-quarks. 
In addition, by studying such quantities one may expose 
the distribution of a hadron's characterising properties 
- such as mass and momentum, linear and angular - 
amongst its QCD constituents; and illuminate strong in- 
teraction properties that were critical in the universe's 
evolution. The upgraded and promised future facilities 
will provide data that should guide the charting process. 
However, to make full use of that data, it will be nec- 
essary to have Poincare covariant theoretical tools that 
enable the reliable study of hadrons in the mass range 1- 
2 GeV and in the vicinity of phase boundaries. Crucially, 
on these domains both confinement and DCSB are ger- 
mane. 

It is known that DCSB; namely, the generation of mass 
from nothing, does take place in QCD. It arises primar- 
ily because a dense cloud of gluons comes to clothe a 
low- momentum quark [l7j . This is readily seen by solv- 
ing the DSE for the dressed-quark propagator; i.e., the 
gap equation, Eq. (|IV.1[) P1 which yields the result illus- 
trated in Fig. III. ll However, the origin of the interaction 
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FIG. II. 1. Dressed-quark mass function, M(p) in Eq. (|IV.2ll : 
solid curves - DSE results, explained in Refs. [lJ,[lB|, "data" - 
numerical simulations of lattice-QCD (NB. m = 70 MeV 
is the uppermost curve and current-quark mass decreases from 
top to bottom.) One observes the current-quark of perturba- 
tive QCD evolving into a constituent-quark as its momentum 
becomes smaller. The constituent-quark mass arises from a 
cloud of low-momentum gluons attaching themselves to the 
current-quark. This is dynamical chiral symmetry breaking 
(DCSB): an essentially nonperturbative effect that generates 
a quark mass from nothing; namely, it occurs even in the 
chiral limit. 

strength at infrared momenta, which guarantees DCSB 
through the gap equation, is currently unknown. This 
relationship ties confinement to DCSB. The reality of 
DCSB means that the Higgs mechanism is largely irrele- 
vant to the bulk of normal matter in the universe. Instead 
the single most important mass generating mechanism 
for light-quark hadrons is the strong interaction effect of 
DCSB0c.g., one can identify it as being responsible for 
98% of a proton's mass QJ, 

There is a caveat; namely, as so often, the pion is ex- 
ceptional. Its mass is given by the simple product of 
two terms, one of which is the ratio of two order param- 
eters for DCSB whilst the other is determined by the 
current-quark mass (Sec. IV A|) . Hence the pion would be 
massless in the absence of a mechanism that can generate 
a current-mass for at least one light-quark. The impact 
of a massless, strongly-interacting particle on the physics 
of the Universe would be dramatic. 



2 In announcing hints for the Higgs boson at CERN, the following 
observations was made: "The Higgs field is often said to give 
mass to everything. That is wrong. The Higgs field only gives 
mass to some very simple particles. The field accounts for only 
one or two percent of the mass of more complex things like atoms, 
molecules and everyday objects, from your mobile phone to your 
pet llama. The vast majority of mass comes from the energy 
needed to hold quarks together inside atoms." This recognises 
the fundamental role of confinement and dynamical chiral sym- 
metry breaking, key emergent phenomena within the strongly 
interacting part of the Standard Model. 



It is natural to ask whether the connection between 
confinement and DCSB is accidental or causal. There 
are models with DCSB but not confinement, however, 
an in-vacuum model with confinement but lacking DCSB 
has not yet been identified (see, e.g., Sees. 2.1 and 2.2 
of Ref. Q). This leads to a conjecture that DCSB is a 
necessary consequence of confinement. It is interesting 
that there are numerous models and theories which ex- 
hibit both confinement and DCSB, and possess an exter- 
nal control parameter such that deconfinement and chi- 
ral symmetry restoration occur simultaneously at some 
critical value of this parameter; e.g., quantum electro- 
dynamics in three dimensions with Nf electrons (2p| - [22l ] , 
and models of QCD at nonzero temperature and chemical 
potential [23l - l32l |. Whether this simultaneity is a prop- 
erty possessed by QCD, and/or some broader class of 
theories, in response to changes in: the number of light- 
quark flavours; temperature; or chemical potential, is a 
longstanding question. 

The momentum-dependence of the quark mass, illus- 
trated in Fig. III. 11 is an essentially quantum field the- 
oretic effect, unrealisable in quantum mechanics, and a 
fundamental feature of QCD. This single curve connects 
the infrared and ultraviolet regimes of the theory, and 
establishes that the constituent-quark and current-quark 
masses are simply two connected points separated by a 
large momentum interval. The curve shows that QCD's 
dressed-quark behaves as a constituent-quark, a current- 
quark, or something in between, depending on the mo- 
mentum of the probe which explores the bound-state con- 
taining the dressed-quark. It follows that calculations ad- 
dressing momentum transfers Q 2 > M , where M is the 
mass of the hadron involved, require a Poincare-covariant 
approach that can veraciously realise quantum field the- 
oretical effects [H] . Owing to the vector-exchange char- 
acter of QCD, covariance also guarantees the existence of 
nonzero quark orbital angular momentum in a hadron's 
rest-frame @,HI|-[36|. 

The dressed-quark mass function has a remarkable 
capacity to correlate and to contribute significantly in 
explaining a wide range of diverse phenomena. This 
brings urgency to the need to understand the relation- 
ship between parton properties in the light-front frame, 
whose peculiar properties simplify some theoretical anal- 
yses, and the structure of hadrons as measured in the 
rest frame or other smoothly related frames. This is a 
problem because, e.g., DCSB, an established keystone of 
low-energy QCD, has not explicitly been realised in the 
light-front formulation. The obstacle is the constraint 
k + := k° + k 3 > for massive quanta on the light front 
[37j ■ It is therefore impossible to make zero momentum 
Fock states that contain particles and hence the vacuum 
is "trivial". On the other hand, it is conceivable that 
DCSB is inextricably tied with the formation and struc- 
ture of Goldstone modes and not otherwise a measurable 
property of the vacuum. This conjecture is being ex- 
plored [38l44l| and is something about which more will 
be written herein (Sec. IV 51) . In addition, parton distri- 
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bution functions, which have a probability interpretation 
in the infinite momentum frame, must be calculated in 
order to comprehend their content: parametrisation is 
insufficient. It would be very interesting to know, e.g., 
how, if at all, the distribution functions of a Goldstone 
mode differ from those of other hadrons HI. 



III. CONFINEMENT 

It is worth stating plainly that the potential between 
infinitely-heavy quarks measured in numerical simula- 
tions of quenched lattice-regularised QCD - the so-called 
static potential - is simply irrelevant to the question of 
confinement in the real world, in which light quarks are 
ubiquitous. In fact, it is a basic feature of QCD that 
light-particle creation and annihilation effects are essen- 
tially nonperturbative and therefore it is impossible in 
prin ciple to compute a potential between two light quarks 

Drawing on a long list of sources; e.g., Refs. [44|- l 47 l . 
a perspective on confinement was laid out in Ref. [48j . 
Confinement can be related to the analytic properties 
of QCD's Schwingcr functions, which are often called 
Euclidean-space Green functions. For example, it can be 
read from the reconstruction theorem [49l. |5(| that the 
only Schwinger functions which can be associated with 
expectation values in the Hilbert space of observables; 
namely, the set of measurable expectation values, are 
those that satisfy the axiom of reflection positivity. This 
is an extremely tight constraint. It can be shown to re- 
quire as a necessary condition that the Fourier transform 
of the momentum-space Schwinger function is a positive- 
definite function of its arguments. This condition sug- 
gests a practical confinement test, which can be used with 
numerical solutions of the DSEs (see, e.g., Sec. III. C of 
Ref. H3 and Sec. IV of Ref. [H). The implications and 
use of reflection positivity are discussed and illustrated 
in Sec. 2 of Ref. [|. 

It is noteworthy that any 2-point Schwinger function 
with an inflexion point at p > must breach the axiom 
of reflection positivity, so that a violation of positivity can 
be determined by inspection of the pointwise behaviour of 
the Schwinger function in momentum space (Sec. IV. B of 
Ref. [I(|). Consider then A(fc 2 ), which is the single scalar 
function that describes the dressing of a Landau-gauge 
gluon propagator. A large body of work has focused on 
exposing the behaviour of A(fc 2 ) in the pure Yang-Mills 
sector of QCD. These studies are reviewed in Ref. [13[. A 
connection with the expression and nature of confinement 
in the Yang-Mills sector is indicated in Fig. IIII.il The ap- 
pearance of an inflexion point in the two-point function 
generated by the gluon's momentum-dependent mass- 
function is impossible to overlook. Hence this gluon can- 
not appear in the Hilbert space of observable states. The 
inflexion point possessed by M(p 2 ), visible in Fig. III.ll 
conveys the same properties on the dressed-quark prop- 
agator. 
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FIG. III.l. A(fc 2 ), the function that describes dressing of 
a Landau-gauge gluon propagator, plotted for three distinct 
cases. A bare gluon is described by A(/c 2 ) — 1/k 2 (the dashed 
line), which is plainly convex on k' 2 £ (0, oo). Such a propa- 
gator has a representation in terms of a non-negative spectral 
density. In some theories, interactions generate a mass in 
the transverse part of the gauge-boson propagator, so that 
A(fc 2 ) = l/(fc 2 +m 2 ), which can also be represented in terms 
of a non-negative spectral density. In QCD, however, self- 
interactions generate a momentum-dependent mass for the 
gluon, which is large at infrared momenta but vanishes in the 
ultraviolet [l3|]. This is illustrated by the curve labelled "IR- 
massive but UV-massless." With the generation of a mass- 
function, A(fc 2 ) exhibits an inflexion point and hence cannot 
be expressed in terms of a non-negative spectral density. 

Numerical simulations of lattice-QCD confirm the ap- 
pearance of an inflexion point in both the dressed-gluon 
and -quark propagators; e.g., see Fig. III.ll and Ref. [l]|. 
The signal is clearest for the gluon owing to the greater 
simplicity of simulations in the pure Yang-Mills sector 
[53l - [55| . We emphasise that this sense of confinement is 
essentially quantum field theoretical in nature. Amongst 
its many consequences is that light-quark confinement 
in QCD cannot veraciously be represented in potential 
models. 

From the perspective that confinement can be related 
to the analytic properties of QCD's Schwinger functions, 
the question of light-quark confinement can be trans- 
lated into the challenge of charting the infrared behav- 
ior of QCD's universal /3-function. (Although this func- 
tion may depend on the scheme chosen to renormalise 
the theory, it is unique within a given scheme Of 
course, the behaviour of the /3-function on the perturba- 
tive domain is well known.) This is a well-posed problem 
whose solution is an elemental goal of modern hadron 
physics (e.g., Refs. [57ri59l ]) and which can be addressed 
in any framework enabling the nonperturbative evalua- 
tion of rcnormalisation constants. It is the /3-function 
that is responsible for the behaviour evident in Fig. III. Il 
and Fig. IIII. 21 below; and one of the more interesting of 
contemporary questions is whether it is possible to re- 
construct the /3-function, or at least constrain it tightly, 
given empirical information on the gluon and quark mass 
functions. 

In this connection we note that the DSEs connect the 
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FIG. III. 2. Upper panel - Rainbow-ladder gluon running- 
mass; and lower panel - rainbow-ladder effective running- 
coupling. Matching related curves in each panel (solid to 
solid or dashed to dashed) gives an equivalent description of 
observables within the rainbow-ladder truncation. 



/3-function to experimental observables. Hence, compar- 
ison between DSE computations and observations of the 
hadron mass spectrum, and elastic and transition form 
factors, can be used to chart /3-function's long-range be- 
haviour. Extant studies show that the properties of 
hadron excited states are a great deal more sensitive to 
the long-range behaviour of the /3-function than those of 
ground states. This is illustrated in Refs. [Ht], H3|, which 
through a study of ground-state, radially-excited and ex- 
otic scalar-, vector- and flavoured-pseudoscalar- mesons 
in rainbow-ladder truncation, which is leading order in 
the most widely used nonperturbative scheme [6ll . l62j , 
produced the effective coupling and running gluon mass 
depicted in Fig. IIII.2l 



IV. GAP AND BETHE-SALPETER EQUATIONS 

In order to proceed it is necessary to describe explicitly 
the best known and simplest DSE. The Dyson or gap 
equation determines how quark propagation is influenced 
by interactions; viz., for a quark of flavour /, 

Sf(p)- 1 =Z 2 (i 1 -p + m h ™) 

+Z 1 J ^2V(p-ff)y7M£/(?)yr£(g, P ), (iv.i) 



where: D^ u is the gluon propagator; r£, the quark-gluon 
vertex; , a symbol that represents a Poincare invari- 
ant regularization of the four-dimensional Euclidean in- 
tegral, with A the regularization mass-scale; m^ m (A), the 
current-quark bare mass; and Zi i 2(C 2 ,A 2 ), respectively, 
the vertex and quark wave-function renormalisation con- 
stants, with £ the renormalisation point - dependence 
upon which is not usually made explicit. 

The gap equation's solution is the dressed-quark prop- 
agator, 

S(p) = -ij -pa v {p 2 X 2 ) + vs{p 2 ,( 2 ) ■ 
1 



z 7 -M(p2,C 2 ) + W,C 2 ) 

z(p 2 X 2 ) 

i-y-p+M(p 2 ) ' 



(IV.2a) 
(IV.2b) 



(IV.2c) 



which is obtained from Eq. pV.l|) augmented by a renor- 
malisation condition. A mass-independent scheme is a 
useful choice and can be implemented by fixing all renor- 
malisation constants in the chiral limitjj 

The mass function, M(p 2 ) = B(p 2 , ( 2 )/A(p 2 , ( 2 ), is in- 
dependent of the renormalisation point, £; and the renor- 
malised current-quark mass, 



m) = Z m (C A)m bm (A) = Z^Z 2 m h f 



bm 



(IV.3) 



wherein Z4 is the renormalisation constant associated 
with the Lagrangian's mass-term. Like the running cou- 
pling constant, this "running mass" is familiar from text- 
books. However, it is not commonly appreciated that 
is simply the dressed-quark mass function evaluated at 
one particular deep spacelike point; viz, 



m C f = M/(C 2 ) 
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(IV.4) 



The renormalisation-group invariant current-quark mass 
may be inferred via 



TO/ 



lim 



1 p 

- In — g- 

2 A 2 



QCD 



Mf{p 2 



(IV.5) 



where j m = 12/(33 — 2Nf): Nf is the number 
of quark flavours employed in computing the running 
coupling; and Aqcd is QCD's dynamically-generated 
renormalisation-group-invariant mass-scale. The chiral 
limit is expressed by 



TO/ = . 



Moreover, 



VC > Aqcd, 



in 



h 



TO/ 2 



See, e.g., Ref. 1631 and references therein; or Ref. [64 
tailed discussion of renormalisation. 



(IV.6) 



(IV.7) 



for a de- 
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However, we would like to emphasise that in the presence 
of DCSB the ratio 



1 > f 



MfAP 2 
M h (p 2 



(IV. 



is not independent of p 2 : in the infrared; i.e., Vp 2 < 
Aq CD , it then expresses a ratio of constituent-like quark 
masses, which, for light quarks, are two orders-of- 
magnitude larger than their current-masses and nonlin- 
ear ly related to them [3 . 

The gap equation illustrates the features and flaws of 
each DSE. It is a nonlinear integral equation for the 
dressed-quark propagator and hence can yield much- 
needed nonperturbative information. However, the ker- 
nel involves the two-point function Z? M „ and the three- 
point function T[ . The gap equation is therefore coupled 
to the DSEs satisfied by these functions, which in turn 
involve higher n-point functions. Hence the DSEs are 
a tower of coupled integral equations, with a tractable 
problem obtained only once a truncation scheme is spec- 
ified. It is unsurprising that the best known truncation 
scheme is the weak coupling expansion, which reproduces 
every diagram in perturbation theory. This scheme is sys- 
tematic and valuable in the analysis of large momentum 
transfer phenomena because QCD is asymptotically free 
but it precludes any possibility of obtaining nonpertur- 
bative information. 

Given the importance of DCSB in QCD, it is significant 
that the dressed-quark propagator features in the axial- 
vector Ward-Takahashi identity, which expresses chiral 
symmetry and its breaking pattern|f| 



P^llfc p ) + * K(0 + m B (Q] i1 9 (fc; P) 



fg I 



= S f 1 (fc+)i75 + ij5S g 1 (k-) 



(IV.9) 



where P = p\ + p-2 is the total-momentum entering the 
vertex and k is the relative-momentum between the am- 
putated quark legsQ In this equation, T^ 9 and T^ 9 
are, respectively, the amputated axial-vector and pseu- 
doscalar vertices. They are both obtained from an in- 
homogeneous Bcthc-Salpeter equation (BSE), which is 
exemplified here using a textbook expression [65|: 

[T 5tl (k;P)] tu = Z 2 [iBltj\tu 



+ [S(q + )T 5fl (q;P)S(q.)] sr K^(q,k;P), (IV.IO) 

Jq 

in which K is the fully-amputated quark-antiquark scat- 
tering kernel, and the colour-, Dirac- and flavour- matrix 



structure of the elements in the equation is denoted by 
the indices r, s, t, u. N.B. By definition, K does not con- 
tain quark-antiquark to single gauge-boson annihilation 
diagrams, nor diagrams that become disconnected by 
cutting one quark and one antiquark line. 

The Ward-Takahashi identity, Eq. (|TV.9|) . entails that 
an intimate relation exists between the kernel in the gap 
equation and that in the BSE. (This is another example of 
the coupling between DSEs.) Therefore an understand- 
ing of chiral symmetry and its dynamical breaking can 
only be obtained with a truncation scheme that preserves 
this relation, and hence guarantees Eq. (|IV.9[) without a 
fine-tuning of model-dependent parameters. 



V. NONPERTURBATIVE TRUNCATION 

Through the gap and Bethe-Salpeter equations the 
pointwise behaviour of the /3-function determines the pat- 
tern of chiral symmetry breaking; e.g., the behaviour in 
Fig. III. II Moreover, the fact that these and other DSEs 
connect the /3-function to experimental obscrvablcs en- 
tails that comparison between computations and obser- 
vations of the hadron mass spectrum, and hadron clastic 
and transition form factors, can be used to constrain the 
/3-function's long-range behaviour. 

In order to realise this goal, a nonperturba- 
tive symmetry-preserving DSE truncation is necessary. 
Steady quantitative progress can be made with a scheme 
that is systematically improvable [(3l|, EH- In fact, the 
mere existence of such a scheme has enabled the proof of 
exact nonperturbative results in QCD. 

Before describing a number of these in some detail, it is 
worth explicating the range of applications. For example, 
there are: veracious statements about the pion tr-term 
18 1, radially-excited an d hybrid pseudoscalar mesons 



66l |67| . heavy- light [68( and heavy-heavy mesons [34| 
novel results for the pion susceptibility obtained via anal- 
ysis of the isovector-pseudoscalar vacuum polarisation 
|69j, which bear upon the essential content of the so- 
called "Mexican hat" potential, which is used in building 
models for QCD; and a derivation [7(| of the Weinberg 
sum rule 711. 



Pseudoscalar meson mass formula 



4 Section IV CI discusses the important differences encountered in 
treating flavourless pseudoscalar mesons. 

5 To be explicit, k = (1 — r])pt + r?P2, with r\ G [0, 1], and hence 
fc+ = pi = k + rjP, k- = P2 = k — (1 — rj)P. In a Poincarc 
covariant approach, such as presented by a proper use of DSEs, 
no observable can depend on rj; i.e., the definition of the relative 
momentum. 



Turning now to a fuller illustration, the first of the 
results was introduced in Ref. [72j ; namely, a mass for- 
mula that is exact for flavour non-diagonal pseudoscalar 
mesons: 
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where: rnj arc the current-masses of the quarks consti- 
tuting the meson; and 



(V.2) 



ipn - = -(0|g/ 2 i75<?/il-ffo 



^2tr C D/ i^l^Sf^q+^H ■ (q;P)S h (q-), (V.3) 



(V.4) 



= Z 4 tr C D l5 S fl (q + )r Ho Aq;P)S f2 (q-), (V.5) 

where Fh _ is the pseudoscalar meson's bound-state 
Bcthc-Salpctcr amplitude: 

T Hq _ (k; P) = 75 [iE Ho _ (fc: P) + 7 • PFh _ (k; P) 
+7 • fc G ffo _ (fc; P) - a^P v H Ha _ (k; P)] , (V.6) 

which is determined from the associated homogeneous 
BSE. N.B. An analogous formula for scalar mesons is 
presented in Ref. [ioj . 

It is worth emphasising that the quark wave- 
function and Lagrangian mass rcnormalisation constants, 
2^2,4 (C, A), respectively, depend on the gauge parame- 
ter in precisely the manner needed to ensure that the 
right-hand sides of Eqs. (|V.3[) , (|V.5[) are gauge- invariant. 
Moreover, Z2(£, A) ensures that the right-hand side of 
Eq. (|V.3|) is independent of both £ and A, so that /h _ 

is truly an observable; and Z^C,, A) ensures that p^ H 
is independent of A and evolves with £ in just the 
way necessary to guarantee that the product p* H 
is renormalisation-point-indcpcndent. In addition, it 
should be noted that Eq. (jV.lj) is valid for every pseu- 
doscalar meson and for any value of the current-quark 
masses; viz., rhf i € [0, 00), i = 1,2. This includes ar- 
bitrarily large values and also the chiral limit, in whose 
neighbourhood Eq. (|V.1[) can be shown [72| to reproduce 
the familiar Gell-Mann-Oakes-Renner relation. 

The axial- vector Ward-Takahashi identity, Eq. PV.9P , 
is a crucial bridge to Eqs. (|V.1|) - (|V.5|) ; and on the way 
one can also prove the following Goldbcrgcr-Treiman-likc 
relations 



72J: 



f O Ho _E Ho _(k;0)=B°(k 2 



(V.7) 



F R (k; 0) + 2/« Q _ F Hq _ (k; 0) = A°(k 2 ) , (V.8) 

G R (k; 0) + 2f Hg _ G Ho _ (k; 0) = ^A°(P) , (V.9) 
H R (k; 0) + 2/° q _ H Ho _ (k; 0) = , (V.10) 

wherein the superscript indicates that the associated 
quantity is evaluated in the chiral limit, and F R , G R , H R 
are analogues in the inhomogeneous axial-vector vertex 
of the scalar functions in the H - -meson's Bethe-Salpeter 
amplitude. 

These identities are of critical importance in QCD. The 
first, Eq. (|V.7|) . can be used to prove that a masslcss 
pseudoscalar meson appears in the chiral-limit spectrum 
if, and only if, chiral symmetry is dynamically broken. 



Moreover, it exposes the fascinating consequence that the 
solution of the two-body pseudoscalar bound-state prob- 
lem is almost completely known once the one-body prob- 
lem is solved for the dressed-quark propagator, with the 
relative momentum within the bound-state identified un- 
ambiguously with the momentum of the dressed-quark. 
This latter emphasises that Goldstonc's theorem has a 
pointwise expression in QCD. The remaining three iden- 
tities are also important because they show that a pseu- 
doscalar meson must contain components of pseudovec- 
tor origin. This result overturned a misapprehension of 
twenty-years standing; namely, that only Er _ (k; 0) is 
nonzero [73|]. These pseudovector components materi- 
ally influence the observable properties of pseudoscalar 
mesons 74-73], Sec. lVIIl as do their analogues in other 
mesons 73- 8U . 



It is natural to reiterate here a prediction for the prop- 
erties of non-ground-state pseudoscalar mesons, which 
follows from the exact results described above; namely, 
in the chiral limit 16a. [671 



, Vn > 1 



(V.ll) 



where n is a state label and n = denotes the ground 
state. This is the statement that Goldstone modes are 
the only pseudoscalar mesons to possess a nonzero lep- 
tonic decay constant in the chiral limit when chiral sym- 
metry is dynamically broken. The decay constants of all 
other pseudoscalar mesons on this trajectory, e.g., radial 
excitations, vanish. On the other hand, in the absence 
of DCSB the leptonic decay constant of each such pseu- 
doscalar meson vanishes in the chiral limit; i.e, Eq. (|V.11|) 
is true Vn > 0. 

From the perspective of quantum mechanics, 
Eq. (|V.11[) is a surprising fact. The leptonic decay 
constant for 5-wave states is typically proportional 
to the wave-function at the origin. Compared with 
the ground state, this is smaller for an excited state 
because the wave-function is broader in configuration 
space and wave-functions are normalised. However, 
it is a modest effect. For example, consider the e + e~ 
decay of vector mesons. A calculation in relativistic 
quantum mechanics based on light-front dynamics [83[ 
yields \ f Pl /f Pa \ = 0.5, consistent with the value inferred 
from experiment and DSEs in rainbow-ladder truncation 
[57I ] : \f Pl /f Po \ = 0.45. Thus, it is not uncommon for 
Eq. (jV.lip to be perceived as "remarkable" or "unbe- 
lievable." Notwithstanding this, in connection with the 
pion's first radial excitation, the value of f Wl = — 2MeV 
predicted in Ref. [6(| is consistent with experiment [84[ 
and simulations of lattice-QCD [82[, as illustrated in 
Fig. IV. 11 It is now recognised that the suppression of / Wl 
is a useful benchmark, which can be used to tune and 
validate lattice QCD techniques that try to determine 
the properties of excited states mesons. 
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Luminiferous Ether 



(V»t) 2 

FIG. V.l. Lattice-QCD results for the ratio of the decay 
constants for the first-excited- and ground-state pseudoscalar 
mesons as a function of the pion mass squared. (Lattice pa- 
rameters: volume = 16 3 x 32; /3 = 5.2, spacing a ~ 0.1 fm, 
two flavours of degenerate sea quarks; Wilson gauge action 
and clover fermions.) The "not improved" results were ob- 
tained from a fermion action with poor chiral symmetry prop- 
erties. In this case \ f-n 1 /f-n \ ~ 0.4, consistent with expecta- 
tions based on quantum mechanics. The "improved" results 
were obtained through implementation of the full ALPHA 
method for the nonperturbative improvement of the fermion 
action, which greatly improves the simulation's chiral sym- 



metry properties. In this case, \f^ 1 /f-n 



0.01. (NB. The 



sign of the ratio was not determined in the simulation but 
must be negative on general grounds [6f|. Figure adapted 
from Ref. fgl.) 



B. Condensates are confined within hadrons 

Dynamical chiral symmetry breaking and its connec- 
tion with the generation of hadron masses was first con- 
sidered in Ref. [85[. The effect was represented as a vac- 
uum phenomenon. Two essentially incquivalent classes of 
ground-state were identified in the mean-field treatment 
of a meson-nucleon field theory: symmetry preserving 
(Wigner phase); and symmetry breaking (Nambu phase). 
Notably, within the symmetry breaking class, each of an 
uncountable infinity of distinct configurations is related 
to every other by a chiral rotation. This is arguably the 
origin of the concept that strongly-interacting quantum 
field theories possess a nontrivial vacuum. 

With the introduction of the parton model for the de- 
scription of deep inelastic scatte ring (DIS), this notion 

that DCSB can be 



was challenged via an argument 
realised as an intrinsic property of hadrons, instead of 
through a nontrivial vacuum exterior to the observable 
degrees of freedom. Such a perspective is tenable because 
the essential ingredient required for dynamical symme- 
try breaking in a composite system is the existence of 
a divergent number of constituents and DIS provided 
evidence for the existence within every hadron of a sea 
of low-momentum partons. This view has, however, re- 
ceived scant attention. On the contrary, the introduction 




Since the Earth is in motion 
the flow of aether across 
the Earth should producea 
detectable "aether wind" 



fffffffffttff 



FIG. V.2. Physics theories of the late 19th century postulated 
that, just as water waves must have a medium to move across 
(water), and audible sound waves require a medium to move 
through (such as air or water), so also light waves require 
a medium, the "luminiferous aether." This was apparently 
unassailable logic until, of course, one of the most famous 
failed experiments in the history of science to date [87i |. 



of QCD sum rules theoretical artifice to estimate 
nonperturbative strong-interaction matrix elements en- 
trenched the belief that the QCD vacuum is characterised 
by numerous distinct, spacetimc- independent conden- 
sates. Faith in empirical vacuum condensates might be 
compared with an earlier misguided conviction that the 
universe was filled with a luminiferous aether, Fig. IV.2l 

Notwithstanding the prevalence of the belief in empir- 
ical vacuum condensates, it does lead to problems; e.g., 
entailing, as explained below, a cosmological constant 
that is 10 46 -times greater than that which is observed 
[39l [H| . This unwelcome consequence is partly respon- 
sible for reconsideration of the possibility that the so- 
called vacuum condensates are in fact an intrinsic prop- 
erty of hadrons. Namely, in a confining theory - and 
confinement is essential to this view - condensates are 
not constant, physical mass-scales that fill all spacetimc; 
instead, they are merely mass-dimensioned parameters 
that serve a practical purpose in some theoretical trun- 
cation schemes but otherwise do not have an existence 
independent of hadrons [H-Sl HI, HI- 

To account for the gross overestimate, recall that the 
expansion of the universe is accelerating (88|. One piece 
of evidence is provided by the observations of type la su- 
pernovae reported in Refs. [HI, [92|]; and another by mea- 
surements of the composition of the universe, which point 
to a missing energy component with negative pressure. 
To explain the latter, the curvature of the universe may 
be characterised by a mass-energy density, pu . There is 
a critical value of this density for which the universe is 
flat: pf- Observations of the cosmic microwave back- 
ground's power spectrum, e.g., Ref. [93j . indicate that 
Oo = pu/ pf = 1 ± 0.04. In a flat universe, the mat- 
ter density, fl m , and the energy density must sum to the 
critical density. However, matter only contributes about 
one-third of the critical density, ttnj = 0.33 ±0.04. Thus 
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two-thirds of the critical density is missing. 

In order to have escaped detection, the missing energy 
must be smoothly distributed. In order not to interfere 
with the formation of structure (by inhibiting the growth 
of density perturbations) the energy density in this com- 
ponent must change more slowly than matter (so that it 
was subdominant in the past). The universe's acceler- 
ated expansion can be accommodated in general relativ- 
ity through the cosmological constant, A, and observa- 
tions determine an associated density, p\ ■ (Recall that 
Einstein introduced the repulsive effect of the cosmolog- 
ical constant in order to balance the attractive gravity 
of matter so that a static universe was possible. How- 
ever, he promptly discarded it after the discovery of the 
expansion of the Universe.) 

It has been suggested that the advent of quantum field 
theory makes consideration of the cosmological constant 
obligatory not optional [88]. Indeed, the only possible 
covariant form for the energy of the (quantum) vacuum; 
viz., 

T^ C = Pvac^ (V.12) 

is mathematically equivalent to the cosmological con- 
stant. The vacuum is [H| "... a perfect fluid and pre- 
cisely spatially uniform ..." so that "Vacuum energy is 
almost the perfect candidate for dark energy." Now, if 
the ground state of QCD is truly expressed in a nonzero 
spacetime-independent expectation value (qq), then the 
energy difference between the symmetric and broken 
phases is of order ATqcd ~ 0.3 GeV, as indicated by 
Fig. III. II One obtains therefrom: 

p QCD = 10 46 p ota fV.13) 

In fact, the discrepancy is far greater if the Higgs vacuum 
expectation value is treated in a similar manner. 

This mismatch has been called "the greatest embar- 
rassment in theoretical physics." However, it vanishes if 
one discards the notion that condensates have a physical 
existence, which is independent of the hadrons that ex- 
press QCD's asymptotically realisable degrees of freedom 
[39( 1 ; namely, if one accepts that such condensates are 
merely mass-dimensioned parameters in one or another 
theoretical truncation scheme. This appears mandatory 
in a confining theory [38], [4^] , a perspective one may em- 
bed in a broader context by considering just what is ob- 
servable in quantum field theory [94| : "... although indi- 
vidual quantum field theories have of course a good deal 
of content, quantum field theory itself has no content 
beyond analyticity, unitarity, cluster decomposition and 
symmetry." If QCD is a confining theory, then the prin- 
ciple of cluster decomposition is only realised for colour 
singlet states [48| and all observable consequences of the 
theory, including its ground state, can be expressed via 
an hadronic basis. This is quark-hadron duality. 

It is worthwhile to recapitulate upon the arguments in 
Refs. [H, E3|. To begin, note that Eq. ((V73l) is the exact 
expression in QCD for the leptonic decay constant of a 



pscudoscalar meson. It is a property of the pion and, 
as consideration of the integral expression reveals, it can 
be described as the pseudovector projection of the pion's 
Bethe-Salpeter wave-function onto the origin in configu- 
ration space. Note that the product ip = STS is called 
the Bethe-Salpeter wave-function because, when a non- 
relativistic limit can validly be performed, the quantity 
tjj at fixed time becomes the quantum mechanical wave- 
function for the system under consideration. (N.B. In the 
neighborhood of the chiral limit, a value for fu can be 

. . . . . O" 1 — I 

estimated via either of two approximation formulae [951 - 
W\ . These formulae both illustrate and emphasize the 
role of /if as an order parameter for DCSB.) 

If chiral symmetry were not dynamically broken, then 
in the neighborhood of the chiral limit /if oc rh [66| . Of 
course, chiral symmetry is dynamically broken in QCD 
p^ - [l7j and so for the ground-state pseudoscalar 

Jim /ff o _(m)=/£ ^0. (V.14) 

Taken together, these last two observations express the 
fact that /if , which is an intrinsic property of the pseu- 
doscalar meson, is a bona fide order parameter for DCSB. 
An analysis within chiral perturbation theory [98| sug- 
gests that the chiral limit value, fjj _ , is ~ 5% below the 
measured value of 92.4 MeV; and efficacious DSE studies 
give a 3% chiral-limit reduction [99} . 

Now, Eq. (jV.5j) is kindred to Eq. (|V.3|I : it is the expres- 
sion in quantum field theory which describes the pseu- 
doscalar projection of the pseudoscalar meson's Bethe- 
Salpeter wave-function onto the origin in configuration 
space. It is thus truly just another type of pscudoscalar 
meson decay constant. 

In this connection it is therefore notable that one may 
rigorously define an "in-meson" condensate; viz. 

- (QhQh) C Ho _ = (01^759/, \H -) (V.15) 

= /h _/4 _ =:4 _(™>- (V.16) 
Now, using Eq. (fV?f|) . one finds [zH 

Jim 4 (m) = Z 4 tr CD [ J^SPfaC) = -<99>? • 

(V.17) 

Hence the so-called vacuum quark condensate is, in fact, 
the chiral-limit value of the in-meson condensate; i.e., it 
describes a property of the chiral-limit pseudoscalar me- 
son. One can therefore argue that this condensate is no 
more a property of the "vacuum" than the pseudoscalar 
mes on's chiral-limit leptonic decay constant. Moreover, 
Ref. [lOCf establishes the equivalence of all three defini- 
tions of the so-called vacuum quark cond ensa t e: a con- 
stant in the operator produ ct ex pansion [lOll Il02| ; via 
the Banks-Casher formula [103]; and the trace of the 
chiral-limit dressed-quark propagator. Hence, they are 
all related to the in-meson condensate via Eq. ([V.17P and 
none is defined essentially in connection with the vacuum. 

It is worth remarking that in the presence of confine- 
ment it is impossible to write a valid nonperturbativc 
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definition of a single quark or gluon annihilation opera- 
tor; and therefore impossible to rigorously define a second 
quantised vacuum (ground state) for QCD upon a foun- 
dation of gluon and quark (quasiparticle) operators. To 
do so would be to answer the question: What is the state 
that is annihilated by an operator which is - as appears 
at present - unknowable? However, with the assumptions 
that confinement is absolute and that it entails quark- 
hadron duality, the question changes completely. In this 
case, the nonperturbative Hamiltonian of observable phe- 
nomena in QCD is diagonalised by colour-singlet states 
alone. The ground state of this nonperturbative strong- 
interaction Hamiltonian is the state with zero hadrons. 
One may picture the creation and annihilation operators 
for such states as rigorously defined via smeared sources 
on a spacctimc lattice. The ground-state is defined with 
reference to such oper ator s, employing, e.g., the Gell- 
Mann - Low theorem [104{ . which is applicable in this 
case because there are well-defined asymptotic states and 
associated annihilation and creation operators. 

In learning that the so-called vacuum quark conden- 
sate is actually the chiral-limit value of an in-pion prop- 
erty, some respond as follows. The electromagnetic ra- 
dius of any hadron which couples to pseudoscalar mesons 
must diverge in the chiral limit. This long- known ef- 
fect arises because the propagation of massless on -shel l 
colo ur-singlet pseudoscalar mesons is undamped |l05l - 
Il08l |. Therefore, does not each pion grow to fill the 
universe; so that, in this limit, the in-pion condensate 
reproduces the conventional paradigm? 

Confinement, again, vitiates this objection. Both DSE- 
and latticc-QCD studies indicate that confinement en- 
tails dynamical mass generation for both gluons and 
quarks, see Sec. lIIII The dynamical gluon and quark 
masses remain large in the limit of vanishing current- 
quark mass. In fact, the dynamical masses are almost 
independent of the current-quark mass in the neighbour- 
hood of the chiral limit. It follows that for any hadron the 
quark-gluon containment-radius docs not diverge in the 
chiral limit. Instead, it is almost insensitive to the mag- 
nitude of the current-quark mass because the dynamical 
masses of the hadron's constituents are frozen at large 
values; viz., 2 — 3 Aqcd- These considerations show that 
the divergence of the electromagnetic radius does not cor- 
respond to expansion of a condensate from within the 
pion but rather to the copious production and subsequent 
propagation of composite pions, each of which contains a 
condensate whose value is essentially unchanged from its 
nonzero current-quark mass value within a containment- 
domain whose size is similarly unaffected. 

There is more to be said in connection with the defini- 
tion and consequences of a chiral limit. Plainly, the exis- 
tence of strongly-interacting massless composites would 
have an enormous impact on the evolution of the uni- 
verse; and it is naive to imagine that one can simply set 
TTiu,d — and consider a circumscribed range of manage- 
able consequences whilst ignoring the wider implications 
for hadrons, the Standard Model and beyond. For exam- 



ple, with all else held constant, Big Bang Nucleosynthe- 
sis is very sensitive to the value of the pion- mass jl09| . 
We are fortunate that the absence of quarks with zero 
current-quark mass has produced a universe in which we 
exist so that we may carefully ponder the alternative. 

The discussion of Ref. 38[ was restricted to pseu- 
doscalar mesons. It is expanded in Ref. [40j via a demon- 
stration that the in-pseudoscalar-meson condensate can 
be represented through the pseudoscalar-meson's scalar 
form factor at zero momentum transfer. With the aid 
of a mass formula for scalar mesons, revealed therein, 
the analogue was shown to be true for in-scalar-meson 
condensates. As argued, the concept is readily extended 
to all hadrons so that, via the zero momentum transfer 
value of any hadron's scalar form factor, one can readily 
extract the value for a quark condensate in that hadron 
which is a measure of dynamical chiral symmetry break- 
ing. 

Given that quark condensates are an intrinsic property 
of hadrons, one arrives at a new paradigm, as observed 
in the popular science press [llOj : "EMPTY space may 
really be empty. Though quantum theory suggests that 
a vacuum should be fizzing with particle activity, it turns 
out that this paradoxical picture of nothingness may not 
be needed. A calmer view of the vacuum would also 
help resolve a nagging inconsistency with dark energy, 
the elusive force thought to be speeding up the expansion 
of the universe." In connection with the cosmological 
constant, putting QCD condensates back into hadrons 
reduces the mismatch between experiment and theory 
by a factor of 10 46 . If tcchnicolour-like theories ar e th e 
correct scheme for extending the Standard Model [l ll| . 
then the impact of the notion of in-hadron condensates 
is far greater still. 



C. Flavourless pseudoscalar mesons 

In connection with electric-charge-neutral pseu- 
doscalar mesons, Eq. (|V.1[) is strongly modified owing to 
the non-Abelian anomaly. This entails that whilst the 
classical action associated with QCD is invariant under 
UA{Nf) (non-Abelian axial transformations generated by 
^o75i where Ao oc diag[l, . . . , Ijv,]), the quantum field 
theory is not. The modification is particularly important 
to properties of 77 and rf mesons. The latter is obviously 
a peculiar pseudoscalar meson because its mass is far 
greater than that of any other light-quark pseudoscalar 
meson; e.g., m v ' = 1.75 m v . We note that the diagram 
depicted in Fig. IV. 31 is often cited as central to a solution 
of the 77-77' puzzle. However, as will become clear below, 
whilst it does contribute to flavour-mixing, the process 
is immaterial in resolving the 77-77' conundrum, as is any 
collection of processes for which the figure may serve as 
a skeleton diagram. 

The correct mass formula for flavourless pseudoscalars 
follows from consideration of the complete UA{Nf) Ward- 
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FIG. V.3. This simple flavour-mixing diagram is immaterial 
to the resolution of the r)-v{ conundrum, as is any collection of 
processes for which the figure may serve as a skeleton diagram. 
(Straight lines denote quarks and springs denote gluons.) 

Takahashi identity: 

~2iM ab T b 5 (k;P) ~ A a (k;P) , (V.18) 

which generalises Eg. (|IV.9|) . In Eg. (|V.18|) . {j a \a = 
0, . . . , Nf — 1} are the generators of U(Nf) in the fun- 
damental representation, orthonormalised according to 
tiT a T b = \8 ab ] the dressed-quark propagator S = 
diag[S' tl , Sd, S s , S C1 S{,, ■ ■ ■} is matrix- valued; and 

M ah =tr F [{J a ,rw c }j h ] , (V.19) 

where 5W is a matrix of renormalised current-quark 
masses and the trace is over flavour indices. 

The final term in the last line of Eq. (|V.18|) expresses 
the non-Abelian axial anomaly. It can be written 

A a (k;P) ^ S^iK) 8 a0 A u (k- 1 P)s- 1 (k^) , (V.20) 

with 

Au{k;P) = Jd^ye^+ ^-^Nf (f°q(x) 0.(0) g(y)) , 

(V.21) 

and since A a ~°(k; P) is a pscudoscalar, it has the general 
form 

A °(fc; P) = 7 ° 75 [*£* (k; P) + 7 • Pf* 0; P) 

+ -/-kk-Pg a (k; P) + a^PvX* (k; P)} . (V.22) 

The matrix element in Eq. (|V.21|) represents an oper- 
ator expectation value in full QCD; the operation in 
Eq. (|V.20[) amputates the external quark lines; and 

O(x) = Ai c [e^F^Fp^x)) = dpKp(x) (V.23) 

is the topological charge density operator, where the 
trace is over colour indices and F^ v = \X a F^ v is the 
matrix- valued gluon field strength tensor. It is plain and 
important that only i? a=0 is nonzero. NB. While Q(x) 
is gauge invariant, the associated Chcrn-Simons current, 
Kp,, is not. Thus in QCD no physical state can couple to 
and hence no state which appears in the observable 
spectrum can contribute to a resolution of the so-called 
C/4(l)-problem; i.e., physical states cannot play any role 
in ensuring that the rf is not a Goldstone mode. 



As described in Scc. IV A[ if one imagines there are Nf 
massless quarks, then DCSB is a necessary and suffi- 
cient condition for the a ^ components of Eq. (|V.18[) to 
guarantee the existence of N 2 — 1 massless bound-states 
of a dressed-quark and -antiquark. However, owing to 
Eq. (|V.20p . a = in Eq. (|V.18[) requires special consider- 
ation. One case is easily covered; viz., it is clear that if 
A = 0, then the a = component of Eq. (|V.18[) is no 
different to the others and there is an additional massless 
bound-state in the chiral limit. 

On the other hand, the large disparity between the 
mass of the Ty'-meson and the octet pseudoscalars sug- 
gests that J4° in real-world QCD. If one carefully 
considers that possibility, then the Goldberger-Treiman 
relations in Eqs. (|V7f|) - ((VTlO]) become fill} 

2fH o .E BS (k;0) = 2B°(k 2 ) - MM), 

(V.24) 

F°(fc;0) + 2f° Ho _F BS (k;0) = A°(k 2 ) - ^(fc;0), 

(V.25) 

G° R (k;0) +2f° Io _G BS (k;0) = 2^A°(fc 2 ) - §„(k;0), 

(V.26) 

H R (k;0)+2f Ho _H BS (k;0) = 0), (V.27) 

It follows that the relationship 

■£x(k;0) = 2B°(k 2 ) , (V.28) 

is necessary and sufficient to guarantee that T^(k;P), 
the flavourless pseudoscalar vertex, does not possess a 
massless pole in the chiral limit; i.e., that there are only 
Nf — 1 massless Goldstone bosons. Now, in the chiral 
limit, B°(k 2 ) 7^ if, and only if, chiral symmetry is dy- 
namically broken. Hence, the absence of an additional 
massless bound-state is only assured through the exis- 
tence of an intimate connection between DCSB and an 
expectation value involving the topological charge den- 
sity. 

This critical connection is further highlighted by the 
following result, obtained through a few straightforward 
manipulations of Eqs. (jV7T8]) . (pT^O]) and (fOTt : 

(qq)°C = - lim K% _ (to) (V.29) 

rh— s-0 

,A 

= - lim Z 4 (C 2 ,A 2 )tr CD / S°(q,0 (V.30) 

A^oo J q 

= ^Jd 4 x (q(x^ l5 q(x)d(0))°. (V.31) 

The absence of a Goldstone boson in the a = channel 
is only guaranteed if this explicit identity between the 
chiral-limit in-meson condensate and a mixed vacuum 
polarisation involving the topological charge density is 
satisfied. 

Mass formulae v alid f or all pseudoscalar mesons have 
also been obtained |l!2j 

fH a - m k- = 2 M ab PH - + HH o- > ( V - 32 ) 
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e.g. 



FIG. V.4. An illustration of the nature of the contribution 
to the Bethe-Salpeter kernel associated with the non-Abelian 
anomaly. All terms have the "hairpin" structure illustrated 
in the lower panel. No finite sum of such intermediate states 
is sufficient. (Straight lines denote quarks, with /i and /2 
independent, and springs denote gluons.) 



where 



/V P M = ^tr/ 7 a ~f^XH Aq;P), (V.33) 

J a 



i/>S!r _(0 = Zi tr / 7 ° 75 XH o- (95 P ) ' ^ V ' 34 ) 

Jq 



Nf 
2 V ' 



I Q,| H -)- (V.35) 



For charged pseudoscalar mesons, Eq. (|V.32j) is equiva- 
lent to Eq. (|V.ip . but the novelty of Eq. (|V.32j) is what 
it expresses for flavourless pseudoscalars. To illustrate, 
consider the case of a U(Nf = 3)-symmetric mass ma- 
trix, in which all Nf = 3 current-quark masses assume 
the single value m£, then this formula yields: 



(V.36) 



Plainly, the rj is split from the Goldstone modes so long 
as n„i 7^ 0. Numerical simulatio ns of lattice-QCD have 
confirmed this identity }l!3l Ill4| . 

It is important to elucidate the physical content of n n i . 
Returning to the definition; viz., 



(o\o.W) 



(V.37) 



it is readily seen to be another type of in-meson conden- 
sate. It is analogous to those discussed in Sec. lVBI but 
in this case the hadron-to-vacuum transition amplitude 
measures the topological content of the rj . One may 
therefore state that the rj is split from the Goldstone 
modes so long as its wavefunction possesses nonzero topo- 
logical content. This is plainly very different to requiring 
the QCD vacuum be topologically nontrivial. 

Within QCD the properties of the rj can be com- 
puted via the BSE, just like other mesons. A nonzero 
value of rin' can be achieved with a Bethe-Salpeter ker- 
nel that contains the contribution depicted in Fig. IV.4l 
because one may argue from Eqs. (|V.21[) and (jV.23l) that 



an anomaly-related contribution to a meson's Bethc- 
Salpctcr kernel cannot contain external quark or anti- 
quark lines that are connected to the incoming lines: 
purely gluonic configurations must mediate, as illustrated 
in Fig. IV. 41 Furthermore, it is straightforward to see that 
no finite sum of gluon exchanges can serve this purpose. 
Indeed, consider any one such single contribution in the 
chiral limit. It will be proportional to the total mo- 
mentum and hence vanish for P = 0, in conflict with 
Eq. (|V.36[) . This lies behind the need for something like 
the Kogut-Susskind ghost; i.e., the coupling of a massless 
axial- vector gauge-like field to the Chern-Simons current, 
which does not appear in the particle spectrum of QCD 
because the current is no t gauge invariant. (See Ref. jl!5| 
and Sec. 5.1 of Ref. fill .) 

It is argued jll7l Ill8| that in QCD with N c colours, 



1 



(V.38) 



and it can be seen to follow from the gap equation, the 
homogeneous BSE and Eqs. ([V33"]) . ()V34|) that 



One thus obtains 



* + 2m(0- 

J r j' J -qi 



#(0 



(V.39) 



(V.40) 



The first term vanishes in the limit N c — > 00 while the 
second remains finite. Subsequently taking the chiral 
limit, the 77' mass approaches zero in the manner char- 
acteristic of all Goldstone modes. (N.B. One must take 
the limit N c — > 00 before t he c hiral limit because the 
procedures do not commute |l!9j |-) The se re sults are re- 
alised in the effective Lagrangian of Ref. mi in a fashion 

that is consistent with all the constraints of the anoma- 
lous Ward identity. N.B. This i s not true of the so-called 

't Hooft determinant [HI [ml HH . 

Implications of t he m ass formula in Eq. (|V.32[) were 
exemplified in Ref. [ 1 1 2| using an elementary dynamical 
model that includes a one-parameter Ansatz for that 
part of the Bethe-Salpeter kernel related to the non- 
Abelian anomaly, an illustration of which is provided in 
Fig. IV. 41 The study compares ground-state pseudoscalar- 
and vector-mesons constituted from all known quarks, 
excluding the i-quark. Amongst the notable results is a 
prediction for the mixing angles between neutral mesons; 
e-g-, 



-15.4°, e, n 



-15.7°. 



(V.41) 



N.B. There are necessarily two mixing angles, with each 
determined at the appropriate pole position in the inho- 
mogeneous vertex. It is interesting that the angles are 
approximately equal and compare well with the value 
inferred from a single mixing angle analysis [123| 9 = 
-13.3° ± 1.0°. 
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It is worth explicating the nature of the flavour-induced 
difference between the 7r° and ir^ masses. If one ignores 
mixing with mesons containing other than u, d-quarks; 
viz., works solely within SU(Nf = 2), then m^o — to t + = 
— 0.04MeV. On the other hand, the full calculation yields 
m^o — m w + = — 0.4MeV, a factor of ten greater, and one 
obtains a ir°-ri mixing angle, whose value at the neutral 
pion mass shell is 

Mmjo) = L2 °- ( V - 42 ) 

For comparison, Ref. |l24j infers a mixing angle of 0.6° ± 
0.3° from a AT-matrix analysis of the process pd — > 
3 He7r°. Plainly, mixing with the 77-meson is the domi- 
nant non-electromagnetic contribution to the n^-ir mass 
splitting. The analogous angle at the rj mass-shell is 

^r,K) = 1-3°. (V.43) 

The angles in Eq. (|V.41[) correspond to 

|tt°) - 0.72 uu - 0.69 dd - 0.013 ss , (V.44) 
[77) ~ 0.53 uu -f- 0.57 dd — 0.63 ss , (V.45) 
|r/) - 0.44 uu + 0.45 dd + 0.78 ss . (V.46) 

Evidently, in the presence of a sensible amount of isospin 
breaking, the ir° is still predominantly characterised by 
J 3 but t here is a small admixture of ss. It is found in 
Ref. [l 12j | that mixing with the 7T° has a similarly modest 
impact on the flavour content of the 77 and 77'. It's effect 
on their masses is far less. 



VI. MANY FACETS OF DCSB 

The importance and interconnection of confinement 
and DCSB are summarised in Sccs. lIIIl IIV1 and some 
of the profound implications of DCSB for pscudoscalar 
mesons are detailed in Sec.fV] The latter could be 
proved owing to the existence of at least one system- 
atic nonperturbative symmetry-preserving DSE trunca- 
tion scheme [6l], [62|]. On the other hand, the practi- 
cal application of this particular scheme has numerous 
shortcomings. For example, at leading-order (rainbow- 
ladder) the truncation is accurate for ground-state 
vector- and electrically-charged pseudoscalar-mesons be- 
cause corrections in these channels largely cancel, ow- 
ing to parameter-free preservation of the Ward-Takahashi 
ident ities . Ho wever, they do not cancel in other chan- 
nels [I25l4i28j . Hence studies based on the rainbow- 
ladder truncation, or low-order improvements thereof, 
have usually provi ded poor results for scalar- and axial- 
vector-mesons [HI, Il29l4l33| , produced masses for exotic 
states that are to o low in comparison with other esti- 
mates [M M M EM Il33| . and exhibit gross sensitiv- 
ity to model parameter s for tensor-mesons [l34| and ex- 
cited states [57], H3, H3, Il35| . In these circumstances one 
must conclude that physics important to these states is 
omitted. One anticipates therefore that significant qual- 
itative advances in understanding the essence of QCD 



could be made with symmetry-preserving kernel Ansatze 
that express important additional nonperturbative ef- 
fects, which are impossible to capture in any finite sum 
of contributions. Such an approach has recently become 
available [93 and will be summarised below. 



A. DCSB in the Bethe-Salpeter kernel 

In order to illustrate the decisive importance of DCSB 
in the Bethe-Salpeter kernel, consider, e.g., flavoured 
pscudoscalar and axial- vector mesons, which appear as 
poles in the inhomogeneous BSE for the axial- vector ver- 
tex, where f,g are flavour labels. An exact form of 
that equation is (k, q are relative momenta, P is the total 
momentum flowing into the vertex, and q± = q ± P/2, 
etc.) 

rg(*5 P) = Z 2l5l(1 - / g 2 D a p{k - q) 

Jq 

+ J g*D a p{k-q)^ la Sf(q + )^-k{lp(k,q-P), 

(VI.l) 

where Ag?g is a 4-point Schwinger function. [The pseu- 
doscalar vertex satisfies an analogue of Eg. (IVI.1[) .] This 
form of the BSE was first written in Ref. |l27| and is illus- 
trated in the lower-panel of Fig. IVI.il The diagrammatic 
content of the right-hand-side is completely equivalent to 
that of Eq. pV.10|) . which is depicted in the upper-panel 
of the figure. However, in striking qualitative opposi- 
tion to that textbook equation, Eq. (|VI.1|) partly em- 
beds the solution vertex in the four-point function, A, 
whilst simultaneously explicating a part of the effect of 
the dressed-quark-gluon vertex. This has the invaluable 
consequence of enabling the derivation of both an inte- 
gral equation for the new Bethe-Salpeter kernel, A, in 
whic h the driving term is the dressed-quark-gluon vertex 
jl27l |. and a Ward-Takahashi identity relating A to that 
vertex (97| . No similar equations have yet been found 
for K and hence the textbook form of the BSE, whilst 
tidy, is very limited in its capacity to expose the effects 
of DCSB in bound-state physics. 

As emphasised above, no study of light-quark hadrons 
is dependable if it fails to comply with the axial-vector 
Ward-Takahashi identity, Eq. pV.9|) . The condition 

P,K{U T, P) + * K (0 + m g (C)] Ag (k, q- P) 
= T f (q + , k+) i l5 + i l5 I*(g_, , (VI.2) 

where A5J is the analogue of Agfg in the pseudoscalar 
equation, is necessary and sufficient to ensure the 
Ward-Takahashi identity is satisfied by the solution of 
Eqs. pVTl) and (1VLT1 
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FIG. VI. 1. Omitting the inhomogeneity, the upper panel il- 
lustrates the textbook form of the Bethe-Salpeter equation, 
Eq. (|IV,10[) . whereas the lower panel depicts the form ex- 
pressed in Eq. (|VI.1|) . The reversal of the total- momentum's 
flow is immaterial here. N.B. In any symmetry-preserving 
truncation, beyond the leading-order identified in Ref. [p^ |; 
i.e., rainbow-ladder, the Bethe-Salpeter kernel is n onplanar 
even if the vertex in the gap equation is planar [l27l | . 



Consider Eq. (|VI.2|) . Rainbow-ladder is the leading- 
order term in the systematic DSE truncation scheme of 
Refs. [6l], HH. It corresponds to T[ = 7„, in which case 
Eq. (fVT2|t is solved by = = A£|. This is the 

solution that indeed provides the rainbow-ladder forms of 
Eq. (|VI.1[) . Such consistency will be apparent in any valid 
systematic term-by-term improvement of the rainbow- 
ladder truncation. 

However, since the two-point functions of elementary 
excitations are strongly modified in the infrared, one 
must accept that the same is generally true for three- 
point functions; i.e., the vertices. Hence the bare vertex 
will be a poor approximation to the complete result un- 
less there are extenuating circumstances. This is readily 
made apparent, for with a dressed-quark propagator of 
the form in Eq. pV.2p . one finds immediately that the 
Ward-Takahashi identity is breached; viz., 

P^7„ + S-\k + P/2) - S~\k - P/2) , (VI.3) 

and the violation is significant whenever and wherever 
the mass function in Fig. III.ll is large. This was actu- 
ally realised early on, with studies of the fermion-gaugc- 
boson vertex in Abelian gauge theories |l36l | that have 
inspired numerous ensuing analyses. The importance 
of this dressing to the reliable computat ion o f ha dron 
physics observables was exposed in Refs. [l37l Il38| . in- 
sights from which have subsequen tly been exploited ef- 
fectively; e.g., Refs. P. I33l.l72l.l8ll fl39l4l43i . 

The most important feature of the perturbative or 
bare vertex is that it cannot cause spin-flip transitions; 
namely, it is an helicity conserving interaction. How- 
ever, one must expect that DCSB introduces nonpertur- 
batively generated structures that very strongly break 
helicity conservation. These contributions will be large 
when the dressed-quark mass-function is large. Con- 
versely, they will vanish in the ultraviolet; i.e., on the 
perturbative domain. The exact form of the vertex con- 



tributions is still the subject of study but their existence 
is model-independent. 

Critical now is a realisation that Eq. (|VI.2[) is far more 
than just a device for checking a truncation's consistency. 
For, just as the vector Ward-Takahashi identity has long 
been used to build Ansat ze for the dressed-quark-photon 
vertex [J [HI EH, [HI, Eq. (fVT2|) provides a tool for 
constructing a symmetry-preserving kernel of the BSE 
that is matched to any reasonable form for the dressed- 
quark-gluon vertex which appears in the gap equation. 
With this powerful capacity, Eq. (|VI.2p achieves a goal 
that has been sought ever since the Bethe-Salpeter equa- 
tion was introduced [65j]. As will become apparent, it 
produces a symmetry-preserving kernel that promises to 
enable the first reliable Poincarc invariant calculation of 
the spectrum of mesons with masses larger than 1 GeV. 

The utility of Eq. ()VT2|) was illustrated in Ref. [97| 
through an application to ground state pseudoscalar and 
scalar mesons composed of equal-mass u- and d-quarks. 
To this end, it was supposed that in Eq. pV.ip one em- 
ploys an Ansatz for the quark-gluon vertex which satisfies 

P ll iTl(k+,k-) = «(P 2 ) [sj 1 (k+) - Sj\k-)] , (VI.4) 

with ® flavour-independent. (N.B. While the true quark- 
gluon vertex does not satisfy this identity, owing to the 
form of the Slavnov- Taylor identity which it does satisfy, 
it is plausible that a solution of Eq. (|VI.4[) can provide 
a rea sonable pointwise approximation to the true vertex 
[H|.) Given Eq. (IVLl) . then Eq. (TVT21 entails (I = q - 
k) 



U p A f 5 l(k,q;P) = 'B(l 2 ) T f 5 9 (q;P) 



r f 5 9 (k;P) 



(VI.5) 



with an analogous equation for P^lpiKt^^k, q; P). This 
identity can be solved to obtain 

Ag(fc, q; P) := * ((k - qf) 75 A^ (fc, q; P) , (VI.6) 

with, using an obvious analogue of Eq. (|V.6|) , 

~K f /(k, q; P) = 2£ [iA E , (q, k; P) + 7 • PA Fs (q, k; P)} 
+ 7/3 S G5 (q, k- P) + 2£p 7 ■£ A Gs (q, k; P) + [ 7/3 , 7 • P] 
x S ff5 (q, k; P) + 2£f,[<y ■ £, 7 • P] A H , (q, k; P) , (VI.7) 

where t = (q + k)/2, E»(g, k; P) = [$(<?; P) + $(fc; P)}/2 
and A*(ff, k\ P) = [$(g; P) - $(fc; P)]/[q 2 - k 2 ]. 

Now, given any Ansatz for the quark-gluon vertex 
that satisfies Eq. (|VI.4I) , then the pseudoscalar analogue 
of Ea.lfVLT]). and Eas. ffVUl) . dVT6|) . WH\ provide a 
symmetry-preserving closed system whose solution pre- 
dicts the properties of pseudoscalar mesons. The relevant 
scalar meson equations are readily derived. With these 
systems one can anticipate, elucidate and understand the 
influence on hadron properties of the rich nonperturba- 
tive structure expected of the fully-dressed quark-gluon 
vertex in QCD: in particular, that of the dynamically 
generated dressed-quark mass function, whose impact is 
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FIG. VI. 2. Dependence of pseudoscalar (upper panel) and 
scalar (lower) meson masses on the current-quark mass, m. 
The Ball-Chiu vertex (BC) result is compared with the 
rainbow-ladder (RL) result. (Figure adapted from Ref. [97f .) 



quashed at any finite order in the truncation scheme of 
Ref. [62j , or any kindred scheme. 

To proceed one need only specify the gap equation's 
kernel because, as noted above, the BSEs are completely 
defined therefrom. To complete the illustration [97fl a 
simplified form of the effective interaction in Ref. [99] 
was employed and two vertex Ansdtze were compared; 
viz., the bare vertex = 7^, which defines the rainbow- 
ladder truncation of the DSEs and o mits vertex dressing; 
and the Ball-Chiu (BC) vertex [l36| |. which nonperturba- 
tively incorporates some of the vertex dressing associated 
with DCSB: 



* r ^(9,fc) = iZ Ag (q 2 ,k 2 )~f^ 

+ 2£„ [»7 • £ A Aa (q 2 , k 2 ) + A Bg (q 2 , fc 2 )]. (VI.8) 

A particular novelty of the study is that one can calcu- 
late the current-quark-mass-dependence of meson masses 
using a symmetry-preserving DSE truncation whose dia- 
grammatic content is unknown. That dependence is de- 
picted in Fig. lVI.2l and compared with the rainbow-ladder 
result. The m-dependence of the pseudoscalar meson's 
mass provides numerical confirmation of the algebraic 
fact that the axial- vector Ward-Takahashi identity is pre- 
served by both the rainbow-ladder truncation and the 



Constituent-quark model 




Pseudoscalar meson: L=0 & S=0 



Scalar meson: L=l & S=l, & L+S=0 



FIG. VI. 3. In constituent-quark-like models, pseudoscalar 
mesons are 1 So states - constituent spins antiparallel and zero 
orbital angular momentum; and scalar mesons are 3 Po states 
- constituent spins parallel and one unit of orbital angular 
momentum. Hence a scalar is a spin and orbital excitation of 
a pseudoscalar meson. 



BC-consistent Ansatz for the Bethe-Salpeter kernel. The 
figure also shows that the axial-vector Ward-Takahashi 
identity and DCSB conspire to shield the pion's mass 
from material v ariation in response to dressing the quark- 
gluon vertex (7I [l27|Jl28 1 . 

As noted in Ref. [9 V| .a rainbow-ladder kernel with re- 
alistic interaction strength yields 



.RL 



2M(0) - m c 



2M(0) 



(0.3 ±0.1) 



(VI.9) 



RL 



which can be contrasted with the value obtained using 
the BC-consistent Bethe-Salpeter kernel; viz., 



£ G <o.i 



(VI.10) 



Plainly, significant additional repulsion is present in the 
BC-consistent truncation of the scalar BSE. 

Scalar mesons arc commonly identified as 3 Po states, 
see Fig. IVI.3l This assignment expresses a constituent- 
quark-model perspective, from which a J PC = ++ 
fermion-antifermion bound-state must have the con- 
stituents' spins aligned and one unit of constituent or- 
bital angular momentum. Hence a scalar is a spin and 
orbital excitation of a pseudoscalar meson. Of course, 
no constitucnt-quark-modcl can be connected system- 
atically with QCD. Nevertheless, the presence of or- 
bital angular momentum in a hadron's rest frame is a 
necessary consequence of Poincare covariance and the 
momentum-dependent vector-boson-exchange character 
of QCD [H, [34|, HH, so there is a realisation in QCD of 
the quark-model anticipation. 

Extant studies of realistic corrections to the rainbow- 
ladder tr unca tion show that they reduce hyperfinc 
splitting [l28j |. Hence, with the comparison between 
Eqs. (|VI.9|) and (|VI. 10[) one has a clear indication that in 
a Poincare covariant treatment the BC-consistent trun- 
cation magnifies spin-orbit splitting, an effect which can 
be attributed to the influence of the quark's dynamically- 
enhanced scalar self-energy Q in the Bethe-Salpeter ker- 
nel. 
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B. Quark anomalous magnetic moments 

It was conjectured in Ref. |)7j that a full realisation of 
DCSB in the Bethe-Salpeter kernel will have a big im- 
pact on mesons with mass greater than 1 GeV. Moreover, 
that it can overcome a longstanding failure of theoretical 
hadron physics. Namely, no extant continuum hadron 
spectrum calculation is believable because all symmetry 
preserving studies produce a splitting between vector and 
axial- vector mesons that is far too small: just one-quart er 
of the experimental value (see, e.g., Refs. [36l . Il30l4l32| ). 
Sign i fican t developments have followed that conjecture 
Il40l |142;1 and will now be related. 



In Dirac's relativistic quantum mechanics, a fermion 
with charge q and mass m, interacting with an elec- 
tromagnetic field, has a magnetic moment [i = q/[2m] 
1471 ]. For the electron, this prediction held true for 



twenty years, until improvements in experimental tech- 
niques enabled the discovery of a small deviation [1481 ] . 
with the moment increased by a multiplicative factor: 
1.00119 ± 0.00005. This correction was explained by the 
first systematic computation usin g renormalised quan- 
tum electrodynamics (QED) jl49j : 



q 

2m 



1 



g 

2?rJ 2m ' 



(VI.ll) 



where a is QED's fine structure constant. The agreement 
with experiment established quantum electrodynamics as 
a valid tool. The correction defines the electron's anoma- 
lous magnetic moment, which is now known with extraor- 
dinary precision and agrees with theory at 0(a 5 ) [l50j . 
The fcrmion-photon coupling in QED is described by: 



/ 



(VI.12) 



where ipi x )y ^{x) describe the fermion field and A^(x) de- 
scribes the photon. This interaction generates the follow- 
ing electromagnetic current for an on-shell Dirac fermion 

(k =p f -p^, 



iqu{p f ) 



7^1 (fc 2 ) + ^-a^k u F 2 {k 2 ) 
2m 



i(Pi), (VI.13) 



where: Fi(k 2 ), F 2 {k 2 ) are form factors; and u(p), u(p) 
are electron spinors. Using their Euclidean space def- 
inition, one can derive a Gordon- identity; viz., with 
21 = Pf+Pi, 

2mu(pf)i~f fl u(p i ) = u(p f ) [21^ + ia^ku] u{pi) . (VI. 14) 

With this rearrangement one sees that for massive 
fcrmions the interaction can be decomposed into two 
terms: the first describes the spin-independent part of 
the fermion-photon interaction, and is common to spin- 
zero and spin-half particles, whilst the second expresses 
the spin-dependent, helicity flipping part. Moreover, one 
reads from Eqs. (|VI. 13|) and (|VI.14[) that a point-particle 
in the absence of radiative corrections has F\ = 1 and 



F 2 = 0, and hence Dirac's value for the magnetic mo- 
ment. The anomalous magnetic moment in Eq. (j VI. lip 
corresponds to F 2 (0) = a/2ir. 

One infers from Eq. (|VI. 14[) that an anomalous contri- 
bution to the magnetic moment can be associated with 
an additional interaction term: 



d 4 x \qijj(x)a l _ LV il}(x)F l _ Ly (x) 



(VI.15) 



where F^{x) is the gauge-boson field strength tensor. 
This term is invariant under local U(l) gauge transfor- 
mations but is not generated by minimal substitution in 
the action for a free Dirac field. 

Consider the effect of the global chiral transformation 
?p(x) — > cxp(id~/ 5 )?p(x). The term in Eq. (|VI.12[) is in- 
variant. However, the interaction of Eq. (|VI.15[) is not. 
These observations facilitate the understanding of a gen- 
eral result: F 2 = for a massless fermion in a quantum 
field theory with chiral symmetry realized in the Wigncr 
mode; i.e., when the symmetry is not dynamically bro- 
ken. A firmer conclusion can be drawn. For m = it 
follows from Eq. (|VI.14|) that Eq. (|VI.12j) does not mix 
with the helicity-flipping interaction of Eq. (|VI. 15[) and 
hence a massless fermion does not possess a measurable 
magnetic moment. 

A reconsideration of Ref. [l49j reveals no manifest con- 
flict with these facts. The perturbative expression for 
F 2 (0) contains a multiplicative numerator factor of m 
and the usual analysis of the denominator involves steps 
that are only valid for m 7^ 0. Fundamentally, there 
is no conundrum because QED is not an asymptotically 
free theory and hence alone does not have a well-defined 
nonperturbativc chiral limit. (N.B. Four-fermion opera- 
tors become relevant in strong-coupling QED and must 
be in cluded in order to obtain a well-defined continuum 
limit fl5ll fl52l ]^ 

On the other hand, in QCD the chiral limit is rig- 
orously defined nonperturbatively J9i|. (It remains to 
be seen whether the theory thus obtained is meaning- 
ful, as indicated in the antepenultimate paragraph of 
Sec. lVBl ) The analogue of Schwinger's one-loop calcu- 
lation can then be carried out to find an anomalous 
c/iromo-magnetic moment for the quark. There are two 
diagrams in this case: one similar in form to that in 
QED; and another o wing to the gluon self- interaction. 
One reads from Ref. jl53| | that the perturbative result 
vanishes in the chiral limit. However, Fig. III. ll demon- 
strates that chiral symmetry is dynamically broken in 
QCD and one must therefore ask whether this affects the 
chromomagnetic moment. 

Of course, it does; and it is now known that this is 
signalled by the appearance of A# B in Eq. (|VI.8[) . If one 
writes the quark-gluon vertex as 



* r ^(p/,K; k) = Ai(p/,pj; k) i^^ 
+ 2£fj, ■ (.\ 2 {pf,pi\k) + \ 3 (p f ,pi;k) 



(VI.16) 



then contemporary simulations of lattice-regularised 
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FIG. VI. 4. Direct anomalous chromomagnetic moment con- 
tribution to the dressed-quark-gluon vertex comput ed i n 
quenched-QCD with current-quark mass m ~ 100 MeV |l54| . 
The one-loop perturbative result is shown for comparison. 
Plainly, the nonperturbative result is two orders-of-magnitude 
larger than the perturbative computation. This level of 
magnification is typical of DCSB; e.g., with a current-quark 
mass of 4 MeV, one obtains M(p 2 = 0) ~ 400 MeV. (Figure 
adapted from Ref. [154| -) 



QCD [Toil and DSE studies [146( agree that 
A 3 (p,P;0)«^i?b 2 ,O 



(VI.17) 



and also on the form of \i(p,p; 0), which is functionally 
similar to A(p 2 , £). However, owing to non-orthogonality 
of the tensors accompanying Ai and A2, it is difficult to 
obtain a lattice signal for A2 . One must therefore consider 
the DSE prediction for A2 in Ref. |l46| more reliable. 

As pointed out above, perturbative massless-QCD con- 
serves helicity so the quark-gluon vertex cannot pcrtur- 
batively have a term with the helicity-flipping character- 
istics of A3. Equation (|VI.17[) is thus remarkable, show- 
ing that the dressed-quark-gluon vertex contains at least 
one chirally-asymmetric component whose origin and size 
owe solely to DCSB; and Sec. lVI Al illustrates that A3 has 
a material impact on the hadron spectr um. 

This reasoning is extended in Ref. jl40| : massless 
fermions in gauge field theories cannot possess an anoma- 
lous chromo/electro-magnetic moment because the term 
that describes it couples left- and right-handed fermions; 
however, if chiral symmetry is strongly broken dynami- 
cally, then the fermions should also posses large anoma- 
lous magnetic moments. Such an effect is expressed in 
the dressed-quark-gluon vertex via a term 



T^(p f , Pi ;k) = a 



r 5 (p f ,Pi,k). (VI. 18) 



That QCD generates a strongly momentum-dependent 
chromomagnetic form factor in the quark-gluon ver- 
tex, T5, w ith a large DCSB-component, is confirmed in 
Ref. [154| . Only a particular kinematic arrangement was 
readily accessible in that lattice simulation but this is 



enough to learn that, at the current-quark mass consid- 
ered: T5 is roughly two orders-of-magnitude larger than 
the perturbative form (see Fig. lVI.4[) ; and 

V>0: hi(p,-p;2p)|>|A3(p,p;0)|. (VI. 19) 

The magnitude of the lattice re sult is co nsistent with 
instanton- liquid model estimates [l55l Il56| . 

This large chromomagnetic moment is likely to have 
a bro ad impact on the properties of light-quark sys- 
tems |l56l Il57l |. In particular, as will be illustrated 
in Sec. lVI Cl it can explain the longstanding puzzle of 
the mass splitting betw een the a\- and p-mesons in the 
hadron spectrum [142f . Here a different novelty will be 
elucidated; viz., the manner in which the quark's chro- 
momagnetic moment generates a quark anomalous elec- 
tromagnetic moment. This demonstration is only possi- 
ble now that t he m ethod of Ref. [9?} is available. It was 
accomplished jl40j | using the same simplification of the 
effective interaction in Ref. |)9| that produced Figs. lVI.2] 

In ord er to understand the vertex Ansatz used in 
Ref. [l40f . it is necessary to return to perturbation the- 
ory. As m entio ned above Eq. (|VI.17|) , one can determine 
from Ref. }l53j | that at leading-order in the coupling, a s , 
the three-gluon vertex does not contribute to the QCD 
analogue of Eq. (|VI.11[) and the Abelian-like diagram pro- 
duces the finite and negative correction (— cv s /[127r]). The 
complete cancellation of ultraviolet divergences occurs 
only because of the dynamical generation of another term 
in the transverse part of the quark-gluon vertex; namely, 

^ cm4 (p/,Pi) = [il7-k + i^a vp £ v k p ]n(p f , Pi ) , (VI.20) 

with Tfj, u = 5^ - k^k u /k 2 , aJ L := T^ u a„. (N.B. The 

tensor denominated Tf, here is labelled T 8 in Refs. [144, 

1 1 

11451 ].) Cognisant of this, one may use a simple Ansatz to 

express the dynamical generation of an anomalous chro- 
momagnetic moment via the dressed-quark gluon vertex; 
viz., 



pacm 



(Pf,Pi) 



r BC 

-pacm4 



(Pf,Pi) + 
(Pf,Pi) 



pacm 
M 

1 pacm 
fJi 



(p f , Pi ), (vi.21) 

5 (p f , Pi ), (VL22) 



with T 5 (pf, Pi ) = (-7/4)A B (p 2 p 2 ) and 



Ti{Pf,Pi) = 7 (z) 



1 - 2r] 

Me 



A B {p 2 ,p 2 )- A A {p),p 2 



The damping factor J [z) 
- 2 ' - 2 2M 2 )/A 2 r , Af = 



(VI.23) 

(1 - cxp(-z))/z, z = 
(pf +pj - 2M|)/A 2 , A, f = lGeV, simplifies numeri- 
cal analysis but is otherwise irrelevant; and Me = { P \p > 
0,p 2 = M 2 (p 2 )} is the Euclidean constituent-quark mass. 
A confined quark does not possess a mass-shell (Sec.lIIip. 
Hence, one cannot unambiguously assign a single value 
to its anomalous magnetic moment. One can nonetheless 
compute a magnetic moment distribution. At each value 
of p 2 , spinors can be defined to satisfy the free-particle 
Euclidean Dirac equation with mass m — > M(p 2 ) =: q, 
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so that 



= u{p f )[F x {k 2 ) lp + ^<J^KF 2 (k 2 )}u( Pt ) (VI.24) 



and then, from Eqs. (|VI.21|t - (|VI.23|t . 



K acm (c) 



<ri - 2^ + 2 ^ R ' 



(VI.25) 



where a\ 



Ea(?,?), = (?,?)> etc - Tne numera- 
tor's simplicity owes to a premeditated cancellation be- 
tween T4 and T5 , which replicates the one at leading-order 
in perturbation theory. Where a comparison of terms 
is possible, this vertex Ansatz i s semi-quantitatively in 
agreement with Refs. |l46l Il54j . However, the pres- 
ence and understanding 
contribution by Ref. [14] 



of the role of rf, cm4 is a novel 
N.B. It is apparent from 
2 in the absence of DCSB, 



Eg. (|VI.25|) that K a 
so that K acm /[2m] —J- in the chiral limit. 

The BSE for the quark-photon vertex can be written 
following the method of Ref. [§7j . Since the method guar- 
antees preservation of the Ward-Takahashi identities, the 
general form of the solution is 



r;i(p/,Pi) = r^ c (p f , Pl ) + Tl(p f , Pl ) , 



(VI.26) 



r J (PfiPi) = tJ-Fi + Oni>k v F 2 + T pp a pv i v £-kF 3 
+ [#y • k + i"ila vp t v k p \F± - UjF 5 



+ £^.k£-kF 6 -e^-£F r 



+l p \a vp l v k p F s , , 



(VI.27) 



where {Fi\i = 1, ... ,8} are scalar functions of Lorentz- 
invariants constructed from pf, pt, k. The Ward- 
Takahashi identity is plainly satisfied; viz., 




p/M E 

FIG. VI.5. Left panel - / 7 (GeV -1 ) in Eq. (|VI.29|l cf. 
(-7/4)A B (p 2 ,p 2 ), both computed using Eg. flVI.21j) and the 
same simplification of the interaction in Ref. [99|]. Right panel 
- Anomalous chromo- and electro-magnetic moment distribu- 
tions for a dressed-quark, computed using Eq. (|VI.30|) . The 
dashed-curve in both panels is the rainbow-ladder (RL) trun- 
cation result. 



k^iT p (p f ,p t ) = k^iT^ (pf,pi) = S (p/) - S (p t ) . 

(VI.28) 

Figure lVL5l depicts the results obtained for the quark's 
anomalous electromagnetic moment form factor 



flip) 



i->f 



(VI.29) 

The result is evidently sizable. It is worth reiterating 
that f-y is completely nonperturbative: in the chiral limit, 
at any finite order in perturbation theory, / 7 = 0. For 
contrast the figure also displays the result computed in 
the rainbow-ladder truncation of QCD's DSEs. As the 
leading-order in a systematic but stepwise symmetry- 
preserving scheme [62| , this truncation only partially ex- 
presses DCSB: it is exhibited by the dressed-quark prop- 
agator but not present in the quark-gluon vertex. In this 
case / 7 is nonzero but small. These are artefacts of the 
truncation that cannot be remedied at any finite order of 
the procedure in Ref. [62j or a kindred scheme. 

Employing Eq. (|VI.24[) , in connection with the dresscd- 
quark-photon vertex, one can write an expression for the 



quark's anomalous electromagnetic moment distribution 



«(0 



2^F 2 + 2^F 4 + A.(0 
a\ + A - A K (c) 



(VI.30) 



where: A„(?) = 2<j 2 <5J 1 - 2qd% - <;F 5 - <j 2 F 7 ; and the 



M(p 2 ) 2 



0. 



Fi are evaluated at p 2 = p 

Plainly, n(s) = in the chiral limit when chiral symmetry 
is not dynamically broken. Moreover, as a consequence 
of asymptotic freedom, k(<;) — > rapidly with increasing 
momentum. The computed distribution is depicted in 
Fig. IVI.5l It yields Euclidean mass-shell values: 



0.44 GeV, 
0.35 GeV, 



"full 



-0.22. 



= °- 45 

-RL 1 = °' 048 - 



cf. M£ L = 0.35 GeV, = . K SL 

(VI.31) 

It is thus apparent that DCSB produces a dressed light- 
quark with a momentum-dependent anomalous chromo- 
magnetic moment, which is large at infrared momenta. 
Significant amongst the consequences is the generation 
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of an anomalous electromagnetic moment for the dressed 
light-quark with commensurate size but o ppos ite sign. 
(N.B. This result was anticipated in Ref. |l58l |. which 
argued that DCSB usually triggers the generation of 
a measurable anomalous magnetic moment for light- 
quarks.) The infrared dimension of both moments is 
determined by the Euclidean constituent-quark mass. 
This is two orders-of-magnitudc greater than the physi- 
cal light-quark current-mass, which sets the scale of the 
perturbative result for both these quantities. 

There are two more notable features; namely, the 
rainbow-ladder truncation, and low-order stepwise im- 
provements thereof, underestimate these effects by an or- 
der of magnitude; and both the T4 and r 5 terms in the 
dressed-quark-gluon vertex are indispensable for a real- 
istic description of hadron phenomena. Whilst a simple 
interaction was used to illustrate these outcomes, they 
are robust. 

These results are stimulating a reanalysis of ha dron 
elas tic and transition electromagnetic form factors jl43l 
Il59j . and the hadron spectrum, results of which will be 
described below. Furthermore, given the magnitude of 
the muon "g^ — 2 anomaly" and its assumed importance 
as an harbinger of physics beyond the Standard Model 
jlOOj ] . it might also be worthwhile to make a quanti- 
tative estimate of the contribution to — 2 from the 
quark's DCSB-induced anomalous moments following, 
e.g., the computational pattern for the hadronic light- 
by-light scattering comp onen t of the photon polarization 
tensor indicated in Ref. 



1.5 
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C. a±-p mass splitting 

The analysis in Ref. |)7j enables one to construct a 
symmetry-preserving kernel for the BSE given any form 
for IV Owing to the importance of symmetries in form- 
ing the spectrum of a quantum field theory, this is a piv- 
otal advance. One may now use all information available, 
from any reliable source, to construct the best possible 
vertex Ansatz. The last two subsections illustrated that 
this enables one to incorporate crucial nonperturbative 
effects, which any finite sum of contributions is incapable 
of capturing, and thereby prove that DCSB generates 
material, momentum-dependent anomalous chromo- and 
electro-magnetic moments for dressed light-quarks. 

The vertex described in Sec. lVI B1 contains a great deal 
of information about DCSB. It is the best motivated 
Ansatz to date, has stimul ated a detailed reanalysis of 
the quark-photon coupling [l45j , and may be used in the 
calculation of the masses of ground-state spin-zero and 
-one light-quark mesons in order to illuminate the impact 
of DCSB on the hadron spectrum. This analysis expands 
significantly on the discussion of scalar and pseudoscalar 
mesons in Sec. lVI Al 

Given a vertex, a prediction for the spectrum fol- 
lows once the gap equation's kernel is specified and the 
Ward-Identity solved for A;^. In the pseudoscalar and 
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FIG. VI. 6. Illustration of the procedure use d to determine 
meson masses, which is fully described in Ref. [l62l ] and anal- 
ogous to the method used in lattice-QCD. Solid curve - ai- 
meson, nonperturbative kernel; dot-dash-dash - ai, kernel de- 
rived from Eq. (|VI.8[) only (Ball-Chiu, BC); and dash - ai, 
kernel derived from just the first term in E q. (|VI.8[) (IBC, 
a minimal renormalization improvement [l63t ] of the leading- 
order - RL, rainbow-ladder - kernel [13]). Dot-dash curve - p- 
meson, nonperturbative kernel; Dot-dash-dot - p, BC-kernel; 
and dotted - p, lBC-kernel. Points - values of 1/T(k = 0; P 2 ) 
in the given channel computed with the kernel described. 
Pade approximants are constructed in each case; and the lo- 
cation of the zero is identified with (— m.^ eaon ). 



axial- vector channels the Ward-Takahashi identity for the 
Bcthe-Salpeter kernel is solved by 

2A 5| g( M) = [fp(q + ,k+) +7 5 f 1 g(g_,fc_)7 5 ] 
1 

•5(|0 



S- 1 {k+) + S-^-k-) 



(k;P) 



5-i(-g + ) + 5-i(g-) 
x [ 75 f p (q+ ,k+) l5 +f p (q_ ,k-)], ( VI.32) 

where T is the chosen Ansatz for the quark-gluon vertex. 
Kernels for ot her channels are readily constructed. 

Reference [142| computes ground-state masses us- 
ing the interaction described in Ref. (57j . which pro- 
duced Fig. IIII. 2| the vertex model explicated in Sec. lVI Bl 
and the method for solving the inhomogeneous Bethe- 
Salpeter equation that is detailed in Ref. \162\ , which en- 
sures one need only solve the gap and Bethe-Salpeter 
equations at spacelike momenta. This simplifies the nu- 
merical problem. To explain, the inhomogeneous BSE 
is solved for the complete Bcthe-Salpeter amplitude in 
a particular channel on a domain of spacelike total- 
momenta, P- > 0. Any bound-state in that channel 
appears as a pole in the solution at P 2 = — JTi^ eson . De- 
noting the leading Chebyshev moment of the amplitude's 
dominant Dirac structure by T(k;P), where k is the rel- 
ative momentum, then 1/T(fc = 0; P 2 ) exhibits a zero at 
(— TO^ esoI1 ). The location of that zero is determined via 
extrapolation of a Pade approximant to the spacclike- 
behavior of 1/T(k = 0; P 2 ). This is illustrated for the p- 
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TABLE VI. 1. Col. 1: Spectrum obtained with the full non- 
perturbative Bethe-Salpeter kernels describe d her ein, which 
express effects of DCSB. The method of Ref. [l62| was used: 
the error reveals the sensitivity to varying the order of Pade 
approximant. Col. 2 - Experimental val ues; c omputed, except 
m CT , from isospin mass-squared averages |l68i ]. Col. 3 - Masses 
determined from the inhomogeneous BSE at leading-order in 
the DSE truncation scheme of Ref. [13] using the interaction 
in Ref. [97j (with this simple ker nel, t he Pade error is negligi- 
ble); and Col. 4 - results in Ref. |l69l |. obtained directly from 
the homogeneous BSE at the same order of truncation. 





Ref. [142] 


Expt. 


RL-Pade 


RL-direct 




0.138 


0.138 


0.138 


0.137 


m p 


0.817 ±0.016 


0.777 


0.754 


0.758 


m a 


0.90 ±0.05 


0.4 - 1.2 


0.645 


0.645 


m ai 


1.30 ±0.11 


1.24 ±0.04 


0.938 


0.927 


m bl 


1.15 ±0.07 


1.21 ±0.02 


0.904 


0.912 


m ai — m p 


0.48 ±0.12 


0.46 ± 0.04 


0.18 


0.17 


rrib 1 — m p 


0.33 ±0.09 


0.43 ± 0.02 


0.15 


0.15 



and ai-channels in Fig. IVI.6l 

A full set of results is listed in Table I VI. 11 wherein 
the level of agreement between Cols. 3 and 4 illustrates 
the efficacy of the method used to compute masses: no 
difference is greater than 1%. Next consider m a and 
compare Cols. 1-3. It is an algebraic result that in the 
RL-truncation of QCD's DSE s, m „ « 2M, where M is 
a constituent-like quark mass 164j . On the other hand, 
incorporating the quark mass function into the Bethe- 
Salpcter kernel via. T^f generates a strong spin-orbit in- 
teraction, which significantly boosts m a [93]. This fea- 
ture is evidently unaffected by the inclusion of r^ cm ; 
i.e., those terms associated with a dressed-quark anoma- 
lous chromomagnctic moment. Since terms associated 
with pion final-state interactions were deliberately omit- 
ted from the nonperturbative kernel derived in Ref. |l42| , 
it is noteworthy that m a in Col. 1 matches estimates for 
the mass of the dressed-quark-core component of the a- 
mes on obtaine d using unitarised chiral perturbation the- 
ory [l65l Il66| . N.B. In addition to providing a width, 
such final-state interacti ons necessarily reduce the real 
part of the mass [H 0, [H3I- 

Now compare the entries in Rows 2, 4-6. The p- 
and ai-mesons have been known for more than thirty 
years and are typically judged to be parity-partners; 
i.e., they would be degenerate if chiral symmetry were 
manifest in QCD. Plainly, they are not, being split by 
roughly 450 McV (i.e., > m p /2). It is suspected that 
this large splitting owes to DCSB. Hitherto, however, 
no symmetry-preserving bound-state treatment could ex- 
plain it. This is illustrated by Cols. 3, 4, which show 
that whilst a good estimate of m p is readily obtained at 
leading-order in the systematic DSE truncation scheme 
of Ref. [Ill, the axial- vector masses are much under- 
estimated. The flaw persists at next-to-leading-order 



p0Hl32l. 

The analysis in Ref. [l42j points to a remedy for this 
longstanding failure. Using the Poincarc-covariant, sym- 
metry preserving formulation of the meson bound-state 
problem enabled by Ref. (97j . with nonperturbative ker- 
nels for the gap and Bethe-Salpeter equations, which in- 
corporate effects of DCSB that are impossible to cap- 
ture in any step-by-step procedure for improving upon 
the rainbow-ladder truncation, it provides realistic esti- 
mates of axi al- ve ctor meson masses. In obtaining these 
results, Ref. |l42| showed that the vertex Ansatz used 
most widely in studies of DCSB, T p c , is inadequate as 
a tool in hadron physics. Used alone, it increases both 



and m ai but yields m ai 



0.21 GeV, qualita- 



tively unchanged from the rainbow-ladder-like result (see 
Fig. lVI.6| ). A good description of axial- vector mesons 
is only achieved by including interactions derived from 
r^ cm ; i.e., connected with the dressed-quark anomalous 
chromomagnetic moment |l40| . Moreover, used alone, 
neither term in r^ cm , Eq. (|VI.22|) . can produce a satis- 
factory result. The full vertex Ansatz and the associated 
gap and Bethe-Salpeter kernels described in Sec. lVI Bl are 
the minimum required. 

Row 5 contains additional information. The leading- 
covariant in the &i-meson channel is 75^^. The appear- 
ance of k p suggests that dressed-quark orbital angular 
momentum will play a significant role in this meson's 
structure, even more so than in the ai-channel for which 
the dominant covariant is 757^. (N.B. In a simple quark- 
model, constituent spins are parallel within the a\ but 
antiparallel within the 61. Constituents of the b\ may 
therefore become closer, so that spin-orbit repulsion can 
exert a greater influence.) This expectation is borne out 
by the following: with the full kernel, is far more 
sensitive to the interaction's momentum-space range pa- 
rameter than any other state, decreasing rapidly as the 
interaction's spatial- variation is increasingly suppressed. 

The results reviewed in this subsection rest on an 
Ansatz for the quark-gluon vertex and whilst the best 
available information was used in its construction, im- 
provement is nonetheless possible. That will involve 
elucidating the role of Dirac covariants in the quark- 
gluo n ver tex which have not yet been considered, as in 
Ref. [l45| . and of resonant contributions; viz., meson loop 
effects that give widths to some of the states consid- 
ered. In cases for which empirical width-to-mass ratios 
arc < 25%, one might judge that such contributions can 
relia bly be obtained via bound-state perturbation the- 
ory [l7fj |. Contemporary studies indicate that these ef- 
fects reduce bound-state masses but the reduction can 
uniformly be compensat ed by a m odest inflation of the 
interaction's mass-scale jl64l Il7 1| . so that the masses in 
Table IVI.ll are semiquantitativcly unchanged. The case 
of the tr-meson is more complicated. However, the predic- 
tion of a large mass for this meson's dressed-quark core 
leaves sufficie nt room fo r a strong reduction by resonant 
contributions [l65l Il66| . 

This section reviewed a continuum framework for com- 
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puting and explaining the meson spectrum, which com- 
bines a veracious description of pion properties with esti- 
mates for masses of light-quark mesons heavier than m p . 
(A contemporary lattice-QCD perspectiv e on this prob- 
lem may be drawn from Refs. jl72l Il73l |.) The method 
therefore offers the promise of a first reliable Poincare- 
invariant, symmetry-preserving computation of the spec- 
trum of light-quark hybrids and exotics; i.e., those puta- 
tive states which are impossible to construct in a quan- 
tum mechanics based upon constituent-quark degrees-of- 
freedom. So long as the promise is promptly fulfilled, the 
approach will provide predictions to guide the forthcom- 
ing generation of facilities. 



VII. PION ELECTROMAGNETIC FORM 
FACTORS 

In charting the long-range interaction between light- 
quarks via the feedback between experiment and theory, 
hadron elastic and transition form factors can provide 
unique information, beyond that obtained through stud- 
ies of the hadron spectrum. 



A. Charged pion 

This is demonstrated very clearly by an analysis of 
the electromagnetic pion form factor, F£ m (Q 2 ), because 
the pion has a unique place in the Standard Model. It 
is a bound-state of a dressed-quark and -antiquark, and 
also that almost-massless collective excitation which is 
the Goldstone mode arising from the dynamical breaking 
of chiral symmetry. This dichotomy can only be under- 
stood through the symmetry-preserving analysis of two- 
body bound-states [72j . Furthermore, the possibility that 
this dichotomous nature could have wide-ranging effects 
on pion properties has made the empirical investigation 
of these properties highly desirable, despite the difficulty 
in preparing a system that can act as a pion target and 
the concomitant c omplexiti es in the interpretation of the 



experiments; e.g., |l74H177j 



The merit of using F° m (Q 2 ) to elucidate the poten- 
tial of an interplay between experiment and nonperturba- 
tive theory as a means of constraining the long-range be- 
haviour of QCD's /3-function is amplified by the existence 
of a prediction that Q 2 F 7T (Q 2 ) « constant for Q 2 ^> m 2 
in a theory whose interactio n is medi ated by massless 
vector-bosons. To be explicit jl78l4l8C| : 



Q 2 FAQ 2 



167r/ 2 a(Q 2 ), 



which takes the value 0.13 GcV 2 at Q 2 



(VII. 1) 

lOGeV 2 if 



one uses the one- loop result a(Q = 10 GcV ) 



0.3. 



The verification of this pre diction i s a strong motiva- 
tion for modern experiment [l74Ul7(| , which can also be 
viewed as an attempt to constrain and map experimen- 
tally the pointwise behaviour of the exchange interaction 
that binds the pion. 



TA BLE V II. 1. Results obtained with [73] oir/tt = 0.93 in 
Eq. dVlL2]|, and (in GeV): m = 0.007, A ir = 0.24, A uv = 
0.905 in Eg. ()VII.3I) . (Dimensioned quantities are listed in 
GeV.) 



m 


E n E p 


M kJ 3, TUtc m p /tt f p 




0.007 


3.568 0.459 1.520 
3.639 0.481 1.531 


0.358 0.241 0.919 0.100 0.130 
0.368 0.243 0.140 0.928 0.101 0.129 



Section IV Al details some extraordinary consequences 
of DCSB, amongst them the Goldberger-Treiman re- 
lations of Eos. (|VT7)) - (|VT0)) . Of these, Eos. (fV8|) 
and (|V.9[) entail that the pion possesses components of 
pscudovector origin which alter the asymptotic form of 
F° m (Q 2 ) by a multiplicative factor of Q 2 cf. the result 
obtained in their absence [74j |. 

QCD-based DSE calculations of F° m (Q 2 ) exist [zllll, 
the most sys tem atic of which [8l| predicted the measured 
form factor |l74| . Germane to this discourse, however, is 
an elucidation of the sensitivity of F° m (Q 2 ) to the point- 
wise behaviour of the interaction between quarks. We 
therefore recapitulate on Refs. [75|, [TtJ, which explored 
how predictions for pion properties change if quarks in- 
teract not via massless vector-boson exchange but instead 
through a contact interaction; viz., 



2n , s x 1 x 47rQ; IR 

9 iVu 5 — q) = V' - 2~ = <V 2~~ 

m G m G 



(vn.2) 



wherein tuq = 0.8 GeV is a gluon mass- scale (such a scale 
is generated dynamically in QCD [l3|, Il8l| ) and aiR is 
a parameter that specifies the interaction strength, and 
proceeded by embedding this interaction in a rainbow- 
ladder truncation of the DSEs. 

In this case, using a confining regularisation scheme 
}l82| . the gap equation, which determines this interac- 
tion's momentum-independent dressed-quark mass, can 
be written 



M = m + M 



4 »" r 

iura 2 - 



C(M 2 ;r ir ,T uv ). 



(VII.3) 



where m is the current-quark mass, r; r = 1/Aj r , r uv = 
1/A UV , and 

C{M 2 - r a , r uv ) =il/ 2 [r(-l, M 2 r 2 v ) - T(-l, M 2 r 2 )], 

(VII.4) 

with T(a, y) being the incomplete gamma- function. Re- 
sults are presented in Table lVILT] 

With a symmetry-preserving regularisation of the in- 
teraction in Eq. (|VIL2[) , the Bethe-Salpeter amplitude 
cannot depend on relative momentum. Hence Eq. (|V.6[) 
becomes 



r ff (P) = 7s 



(VII.5) 
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and the explicit form of the model's ladder BSE is 
4a m 



FAP) 



37T?T1q 



%.ee Kef Ek{P) 

X.FE X.FF \ [ Fir{P) 

(VII.6) 

where, with m — = P 2 , anticipating the Goldstone 
character of the pion, 



Xee = C{M 2 ;t?,tI 



■Kef = 0, 



1%-fe = Ci{M 2 ] t£, ri v ) , Kff = -2Kfe ■ 



(VII.7) 



and C\{z) = —zc'(z), where we have suppressed the 
dependence on T; r uv . The solution of Eg. (|VII.6|) gives 
the pion's chiral-limit Bethe-Salpeter amplitude, which, 
for the computation of observables, should be normalised 
canonically; viz., 



P^ = N c tv 



^ r ^- p )^^ + P)r.(P)S(<z). 

(VII.8) 
(VII.9) 



Hence, in the chiral limit, 

1 = Si Ci(M 2 ;T*,r* v )iM^ - 2F n ], 
and the pion's leptonic decay constant is 

/° = ^—CiCM'jt&t*,)^ - 2F n ] . 



(VII. 10) 



If one has preserved Eq. pV.9|) , then, for m = in the 
neighbourhood of P 2 — 0, the solution of the axial- vector 
BSE has the form: 

r 5/ ,(fc+, k) = ^ 2f° r w (P) + muMP) (vn.ii) 

and the following subset of Eqs. (|V77|> - (|V.10|) will hold: 



f°E^ = M, 2^ + F R = l 



(VII.12) 



Hence F„ (P) is necessarily nonzero in a vector exchange 
theory, irrespective of the pointwise behaviour of the in- 
teraction. It has a measurable impact on the value of /V 
and on the form factor, as we shall see. 

Based upon these results, one can proceed to com- 
pute the electromagnetic pion form factor in the gen- 
eralised impulse approximation [75l l77j: i.e., at leading- 
order in a sy mmetry-preserving DSE truncation scheme 
[74l , l8ll Il38| |. Namely, for an incoming pion with mo- 
mentum pi = K — Q/2, which absorbs a photon with 
space- like momentum Q, so that the outgoing pion has 
momentum P2 = K + Q/2, 



2K,Fr(Q 2 ) = 2N c J ^tr D 



iT n (—p 2 )S(t +P2) 



x*7 M i, r(Q 2 )'S(t + P 1 )*T w (pi)5(t) 



(VII. 13) 
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FIG. VII. 1. Dressing function for the transverse piece of the 
quark-photon vertex, computed using the parameter values 
described in connection with Table lVII.il The pole associated 
with the ground-state vector meson is clear. 



where Pt(Q 2 ), depicted in Fig. IVII.ll describes the full 
extent of dressing on the quark-photon vertex produced 
by a contact interaction in the rainbow-ladder truncation 
[77j . The form factor is expressible as follows: 

Fr(Q 2 ) = Pt(Q 2 ) ^ cn ;(Q 2 ), (vn.i4) 

KjiQ 2 ) = F Zee + K m EF + F™ FF , (VII.15) 
= E c n E<T EE (Q 2 ) + E^T EF (Q 2 ) 
+F^T FF (Q 2 ), (VH.16) 

where F° m ^(Q 2 ) is that part of the form factor produced 
by the undressed quark-photon vertex and the functions 
T 71 " have simple algebraic forms in this calculation ffl\ . 

In the upper panel of Fig. IVII.2l we present F^(Q 2 ) 
and the three separate contributions defined in 
Eq. (|VII.16|) . It is evident that, in magnitude, F™ EF 
contributes roughly one-third of the pion's unit charge. 
This could have been anticipated from Eq. (|VII.9|) . More 
dramatically, perhaps: the interaction in Eq. (|VII.2I) gen- 
erates 



F° m (Q 2 -> 00) = constant. 



(VII.17) 



Both results originate in the nonzero value of F„(P), 
which is a straightforward consequence of the symmetry- 
preserving treatment of a vector exchange theory [72j. 
Equation (|VII.17|) should not come as a surprise: with a 
symmetry-preserving regularisation of the interaction in 
Eq. (|VII.2[) . the pion's Bethc-Salpctcr amplitude cannot 
depend on the constituent's relative momentum. This is 
characteristic of a pointlikc particle, which must have a 
hard form factor. The lower panel of the figure illustrates 
that the necessary inclusion of Pt(Q 2 ) is critical in the 
timelike region and has a measurable quantitative impact 
for a significant range of spacelike momenta. It does not, 
however, affect the truly ultraviolet behaviour. 

In Fig. IVII.l3l we compare the form factor computed 
from Eq. (|VIL2[) with contemporary experimental data 
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FIG. VII. 2. Upper panel. F"^(Q 2 ) and the separate contribu- 
tions introduced in Eq. (|VII.16[) . F° m (Q 2 = 0) = 1, without 
fine-tuning, because a symmetry-preserving regularisation of 
the interaction in Eq. (|VII.2|) was implemented. Lower panel. 
F% m (Q 2 ) computed in rainbow-ladder truncation from the in- 
teraction in Eq. (|VII,2[) : solid curve - fully consistent, i.e., 
with a dressed-quark-photon vertex so that the p-pole ap- 
pears; and dashed curve - computed using a bare quark- 
photon vertex, namely, F^(Q 2 ). Dotted curve - fit to the 
result in Ref. [8l|, which also included a consistently-dressed 
quark-photon vertex and serves to illustrate the trend of con- 
temporary data. 



[174^176} and a QCD-based DSE prediction jUJ. Both 
the QCD-based result and that obtained from the 
momentum-independent interaction yield the same val- 
ues for the pion's static properties [7a - [77j |. However, for 
Q 2 > the form factor computed using ~ l/fc 2 -vector- 
boson exchange is immediately distinguishable empiri- 
cally from that produced by a momentum-independent 
interaction. Indeed, the figure shows that for F% m , exist- 
ing experiments can already distinguish between different 
possibilities for the quark-quark interaction. 

Combining Figs. lVlOl and IVII.3I it becomes apparent 
that F^ n EE is only a good approximation to the net pion 
form factor for Q 2 < M 2 . F°™ EE and F° m EF evolve with 
equal rapidity - there is no reason for this to be otherwise, 
as they are determined by the same mass-scales - but a 
nonzero constant comes quickly to dominate over a form 
factor that falls swiftly to zero. 



FIG VII 3. Solid curve : Q 2 F 7V , I /{Q 2 ) obtained with 
Eq. (|VII.2I) . Dashed curve : DSE prediction which 
employed a momentum-dependent renormalisation-group- 
improved gluon exchange interaction. For Q 2 > 0.17 GeV 2 w 
M 2 , marked by the vertical dotted line, the contact interac- 
tion result for F^(Q 2 ) differs from that in Ref. [HI by more 
than 20%. The data are from Refs. [l7^17r3 ]. 



It is plain now that when a momentum-independent 
vector-exchange interaction is regularised in a symmetry- 
preserving manner, the results are directly comparable 
with experiment, computations based on well-defined and 
systematically- improvable truncations of QCD's DSEs 
[81j . and perturbative QCD. In this context it will be 
apparent that a contact interaction, whilst capable of 
describing pion static properties well, Table rVILll gen- 
erates a form factor whose evolution with Q 2 devi- 
ates markedly from experiment for Q 2 > 0.17 GeV 2 ps 
M 2 and produces asymptotic power-la w behav iour, 
Eq. (|VII.17|) . in serious conflict with QCD fl78l - fl80l ] . 

In that connection Fig. III.ll and Eqs. (|V.7|t - (|V.10|) 
are relevant again. In the electromagnetic elastic scat- 
tering process, the momentum transfer, Q, is primarily 
shared equally between the pion's constituents because 
the bound-state Bethe-Salpeter amplitude is peaked at 
zero relative momentum. Thus, one can consider k ~ 
Q/2. The Goldberger-Treiman-like relations express a 
mapping between the relative momentum of the pion's 
constituents and the one-body momentum of dressed- 
quark; and the momentum dependence of the dressed- 
quark mass function is well-described by perturbation 



theory when k 2 > 2 GeV 2 . Hence, one should expect 
a perturbative-QCD analysis of the pion form factor to 
be valid for k 2 = Q 2 /4 > 2 GeV 2 ; i.e., 



f: 



cmpQCD 



(Q 2 ) for Q 2 > 8 GeV 2 . (VII.18) 



This explains the result in Ref. }74| . A similar argument 
for baryons suggests that the nuclcon form factors should 
be perturbative for Q 2 > 18 GeV 2 . 

It is worth reiterating that the contact interaction 
produces a momentum-independent drcsscd-quark mass 
function, in contrast to QCD-based DSE studies @, UH 
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and latticc-QCD [lyj. This is fundamentally the origin of 
the marked discrepancy between the form factor it pro- 
duces and extant experiment. Hence Refs. [7514771 ] high- 
light that form factor observables, measured at an up- 
graded JLab, e.g., are capable of mapping the running of 
the dressed-quark mass function. Efforts are underway 
to establish the signals of the running mass in baryon 
clastic and transition form factors. They are reviewed in 
SecKl 



B. Neutral pion 

The process 7*7 — s- 7r° is also of great interest because 
in order to explain the associated transition form factor 
within the Standard Model on the full domain of momen- 
tum transfer, one must combine, using a single internally- 
consistent framework, an explanation of the essentially 
nonpcrturbative Abelian anomaly with the features of 
perturbative QCD. The case for attempting this received 
a significant boost wi th th e publication of data from the 
BaBar Collaboration |l83| because, while they agree with 
earlier experiments o n th e ir co mmon domain of squared- 
momentum-transfer [184 Il85| . the BaBar data arc un- 
expectedly far above the prediction of perturbative QCD 
at larger values of Q 2 . 

This so-called "Babar anomaly" was considered in 
Ref. [76| , wherein it is argued that in fully-self-consistent 
treatments of pion: static properties; and clastic and 
transition form factors, the asymptotic limit of the prod- 
uct <3 2 G 7 » 77r o (Q 2 ), which is determined a priori by 
the interaction employed, is not exceeded at any finite 
value of spacelike momentum transfer: the product is a 
monotonically-incrcasing concave function. A consistent 
approach is one in which: a given quark-quark scattering 
kernel is specified and solved in a well-defined, symmetry- 
preserving truncation scheme; the interaction's parame- 
ter (s) are fixed by requiring a uniformly good description 
of the pion's static properties; and relationships between 
computed quantities are faithfully maintained. 

Within such an approach it is nevertheless possible for 
Q 2 F° m (Q 2 ) to exceed its perturbativc-QCD asymptotic 
limit because the leading-order matrix-element involves 
two Bethe-Salpcter amplitudes. This permits an interfer- 
ence between dynamically-generated infrared mass-scales 
in the computation. Moreover, for F° m (Q 2 ) the pertur- 
bative QCD limit is more than an order-of-magnitudc 
smaller than m 2 . Owing to the proximity of the p- meson 
pole to Q 2 = 0, the latter mass-scale must provide a fair 
first-estimate for the small-Q 2 evolution of F° m (Q 2 ). A 
monopole based on this mass-scale exceeds the pQCD 
limit VQ 2 > 0. For the transition form factor, however, 
the opposi te is true because m 2 is less-than the pQCD 



0.30 



limit; viz. [18C 

lim Q 2 2G(Q 2 ) 



(VII. 19) 



0.05 



0.00 



^(x) = constant 




The vector current-current contact-interaction can- 
vassed in this Section may be described as a vector-boson 
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FIG. VII. 4. Q 2 -weighted_7*7 — ► 71-0 transition fo rm f actor. 
Data: red circles, Ref. [18^ ]: green diamonds, Ref. [l84l ]: and 
blue squares, Ref. [l85( |. Solid curve - Q 2 G(Q 2 ) computed 
using the symmetry-preserving, fully-self-consistent rainbow- 
ladder treatment of the contact interaction in Eq. (IVII.2II . 
which produces a pion distribution amplitude 4>tv(x) = 
constant; and dot-dashed curve - fit to the 7*7 — > tv° transi- 
tion form factor comp uted in a QCD-based rainbow-ladder- 
truncation DSE study [l86t | . Both curves have been divided 
by (27r 2 /^-) in order to match the data's normalisation. 



exchange theory with vector-field propagator (l/fc 2 ) K , 
k — Q. It was shown [76[ that the consistent treatment 
of such an interaction produces a 7*7 — > ir° transition 
form factor that disagrees with all available data. On 
the other hand, precisely the same treatment of an in- 
teraction which preserves the one-loop renormalisation 
group behaviour of QCD, produces a form factor in 
good agreement with all b ut the largc-Q 2 data from the 
BaBar Collaboration [l83| . These points are illustrated 
in Fig. lVTL4l 

Studies exist which interpret the BaBar data as an 
indication that the pion's distribution amplitude, 4>tt(%), 
deviates dramatically from its QCD asymptotic form, in- 
deed, that (j>Tr(x) = constant , or is at least flat and nonva- 
nishing at x = 0, 1 fl87l . I l88j |. However, it has often been 
explained 0, [75H77 . Il89l ] that such a distribution ampli- 
tude characterises an essentially-pointlikc pion; and, as 
we have seen, when used in a fully-consistent treatment, 
it produces results for pion clastic and transition form 
factors that are in striking disagreement with experiment. 
Reiterating, a bound-state pion with a pointlike compo- 
nent will produce the hardest possible form factors; i.e., 
form factors which become constant at large-Q 2 . 

On the other hand, QCD-based studies produce soft 
pions, a valence-quark distribution amplitude for the pion 
that vanishes as ~ (1 — x) 2 for x ~ 1, and results that 
agree well with the bulk of existing data. We will return 
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to this in the next section. 

The analysis in Ref. [76| shows that the large-Q 2 BaBar 
data is inconsistent with QCD and also inconsistent with 
a vector current-current contact interaction. It supports 
a conclusion that the large-Q 2 data reported by BaBar 
is not a true representation of the 7*7 — > ir° tran sition 
form factor, a perspective also developed elsewhere [l90l — 
Il93j . There is experimental evidence in support of this 
view; namely, the 7* —> 777 and 7* — > 7/7 transition form 
factors have also been measured by the BaBar Collabo- 
ration [HH, at Q 2 = 11 2 GeV 2 , and in these cases the 
results from CLEO jl85| and BaBar arc fully consistent 
with perturbative-QCD expectations. 



VIII. PION AND KAON VALENCE-QUARK 
DISTRIBUTIONS 



Quarks were discovered in deep inelastic scattering 
(DIS) experiments at SLAC, performed during the period 
1966-1978. DIS is completely different to clastic scatter- 
ing. In this process, one disintegrates the target instead 
of keeping only those events in which it remains intact. 
On a well-defined kinematic domain; namely, the Bjorken 
limit: 



p 



•^Bjorken •- 



fixed, 



2P- q 

(VIII. 1) 

where P is the target's four-momentum and q is the mo- 
mentum transfer, the cross-section can rigorously be in- 
terpreted as a measurement of the momentum-fraction 
probability distribution for quarks and gluons within the 
target hadron: q(x), g{x). These quantities describe the 
probability that a quark/gluon within the target will 
carry a fraction x of the bound-state's momentum, as 
defined in the infinite-momentum or light-front frame. 
(The light-front formulation of quantum field theory is 
built up on D irac's front form of relativistic quantum dy- 
namics |l95| .) 

The past forty years have seen a tremendous effort to 
deduce the parton distribution functions (PDFs) of the 
most accessible hadrons - the proton, neutron and pion. 
There are many reasons for this long sustained and thriv- 
ing interest @ but in large part it is motivated by the 
suspected process-independence of the usual parton dis- 
tribution functions and hence an ability to unify many 
hadronic processes through their computation. In con- 
nection with uncovering the essence of the strong inter- 
action, the behaviour of the valence-quark distribution 
functions at large Bjorken- x is most relevant. Further- 
more, an accurate determination of the behavior of distri- 
bution functions in the valence region is also important to 
high-energy physics. Particle discovery experiments and 
Standard Model tests with colliders are only possible if 
the QCD background is completely understood. QCD 
evolution, apparent in the so-called scaling violations by 




Pion 



Proton 



FIG. VIII. 1. ttN Drell-Yan process, in which, e.g., a valence 
antiquark from the pion annihilates with a valence quark in 
the nucleon to produce a fi + fi~ pair. 



parton distribution functions^ entails that with increas- 
ing center-of-mass energy, the support at large- a; in the 
distributions evolves to small- a; and thereby contributes 
materially to the collider background. N.B. The nucleon 
PDFs are now fairly well determined for x < 0.8 but the 
pion and kaon PDFs remain poorly known on the entire 
domain of x. 

Owing to the dichotomous nature of Goldstone bosons, 
understanding the valence-quark distribution functions 
in the pion and kaon is of great importance. Moreover, 
given the large value of the ratio of s-to-w current-quark 
masses, a comparison between the pion and kaon struc- 
ture functions offers the chance to chart effects of explicit 
chiral symmetry breaking on the structure of wou l d-be 
Goldstone modes. There is also the prediction [l96l Il97j 
that a theory in which the quarks interact via 1/fc 2 - 
vector-boson exchange will produce valence-quark distri- 
bution functions for which 



q v (x) cx {l-x) 2+ ~* , x > 0.85. 



(VIII.2) 



where 7 > is an anomalous dimension that grows 
with increasing momentum transfer. (See Sec. VLB. 3 of 
Ref. [1] for a detailed discussion.) 

Owing to the absence of pseudoscalar meson targets, 
experimental knowledge of the parton structure of the 
pion and kaon arises primarily from pionic or kaonic 
Drell-Yan processes, illust rated in Fig. | yiII.l[ involving 
nucleons in heavy nuclei [l77l Il98l420l| . Theoretically, 
given that DCSB plays a crucial role in connection with 
pseudoscalar mesons, one must employ an approach that 
realistically expresses this phenomenon. The DSEs there- 
fore provide a natural framework: studies of the pion 
and kaon exist and will be reviewed here. The first 
[l89j computed pion PDFs, using efficacious parametrisa- 
tions of both the B ethe-Salp e ter a mplitude and dressed- 
quark propagators [HI HqI HqI . The second [zl em- 
ployed direct, numerical DSE solutions in the compu- 
tation of the pion and kaon PDFs, adapting the ap- 
proach employed in successful pred ictions of electromag- 
netic form factors (6?i H3, HH, Il86| ; and also studied the 



6 DGLAP evolution is described in Sec.IID of Ref. The evolu- 
tion equations are derived in perturbative QCD and determine 
the rate of change of parton densities when the energy-scale cho- 
sen for their definition is varied. 
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FIG. VIII. 2. Pion valence quark distribution function 
evolved to (5.2 GeV) 2 . Solid curve - full DSE calculation [H; 
dot- dashed cu rve - semi-phenomenological DSE-based calcu- 
latio n in Ref. [l89( |; filled circles - experimental data from 
Ref. HqJ, at scale (4.05 GeV) 2 ; dashed curve - NLO reanal- 
ysis of the experimental data [lTTt ] : and dot- dot- dashed curve 
- NLO reanalysis of exp erimental data with inclusion of soft- 
gluon resummation [2041 ]. 

ratio UK{x)/u 7r (x) in order to elucidate aspects of the 
influence of an hadronic environment. 

In rainbow-ladder truncation, one obtains the pion's 
valence-quark distribution from 

1 f d i l 

= -2 y (^}i trcd F-^-p) 

x S u (£) r n (£; x) S u (£) T„(£, P) S d (£ - P)} , (VIII.3) 

wherein the Bethe-Salpeter amplitudes and dressed- 
quark propagators are discussed above and T n (£;x) is 
a generalization of the dressed-quark-photon vertex, de- 
scribing a drcssed-quark scattering from a zero momen- 
tum photon. It satisfies a BSE (here with a rainbow- 
ladder kernel) with the inhomogeneous term ■ n6(£ ■ 
n — xP ■ n). In Eq. (|VIII.3|) . n M is a light-like vector 
satisfying n 2 = 0. In choosing rainbow-ladder trunca- 
tion one implements a precise parallel to the symmetry- 
preserving treatment of the pion charge form factor at 
Q 2 = 0, wherein the vector current is conserved by use 
of ladder dynamics at all three vertices and rainbo w dy - 
namics for all three quark propagators [lH [8(1 HH Il38| . 
Equation (|VIH.3[) ensures automatically that 

(x° f }:=[ dxq}(x) = 1 for / = u, d, (VIII.4) 
Jo 

since / dx T n (£; x) gives the Ward-identity vertex and the 
Bcthc-Salpeter amplitudes arc canonically normalised. 

Figure IVIII.2I displays the DSE results for the pion's 
valence it-quark distribution, evolved from a resolving 
scale Ql = (0.6 GeV) 2 to Q 2 = (5.2 GeV) 2 using leading- 
order DGLAP evolution (Sec.IID of Ref. @), and a com- 
parison with ttN Drell-Yan data [20 lj at a scale Q 2 ~ 
(4.05 GeV) 2 , inferred via a leading-order analysis. The 



computation's resolving scale, Qq, was fixed by match- 
ing the {x 11 )^ momen ts fo r n = 1, 2, 3 to an experimental 
analysis at (2 GeV) 2 [2051 ]. 

It is notable that at Qo the DSE results yield 

2(xy Qo =0.7, 2<*>g =0.8. (VIII.5) 

(For comparison, the par amet rised valence-like pion par- 
ton distributions of Ref. [206| yield a gluon momentum 
fraction of (x g )Q a=0 51 = 0.3.) In each case the remain- 
der of the hadron's momentum is carried by gluons, which 
effect binding of the meson bound state and are invisi- 
ble to the electromagnetic probe. Some fraction of the 
hadron's momentum is carried by gluons at all resolving 
scales unless the hadron is a point particle Indeed, 
it is a simple algebraic exercise to demonstrate that the 
only non-increasing, convex function which can produce 
(x°) = 1 and (x) = i, is the distribution u(x) = 1, 
which is uniquely connected with a pointlike meson; viz., 
a meson whose Bcthc-Salpctcr amplitude is momentum- 
independent. Thus Eqs. ([VIII. 5[) are an essential conse- 
quence of momentum conservation. 

Whilst the DSE results in Fig. lVIII.^I are both consis- 
tent with Eq. (|VIII.2|) ; i.e., they produce algebraically the 
precise behaviour predicted by perturbative QCD, on the 
valence-quark domain - which is uniquely sensitive to the 
behaviour of the dressed-quark mass- function, M{p 2 ) - it 
is evident that they di sagr ee markedly with the Drell-Yan 
data reported in Ref. |20l| . This tension was long seen as 
a crucial mystery for a QCD description of the lig htest 
and subtlest hadron Q . Its re-emergence with Ref. Il89ll 
motivated a NLO reanalysis of the Drell-Yan data |l77| , 
the result of which is also displayed in Fig. lVIlO] At 
NLO the extracted PDF is softer at high- a; but the dis- 
crepancy nevertheless remains. To be precise, Ref. [l77j ] 
determined a high-x exponent of j3 ~ 1.5 whe reas the 
exponents produced by the DSE studies (78l . ll89l ] are, re- 
spectively, 2.1 and 2.4 at the common model scale. They 
do not allow much room for a harder PDF at high-x. 

Following the highlighting of this discrepancy in 
Ref. Q, a resolution of the conflict between da ta a nd 
well-constrained theory was proposed. In Ref. [204| a 
long-overlooked effect was incorporated; namely, "soft- 
gluon resummation." With the inclusion of this next-to- 
leading-logarithmic threshold resummation effect in the 
calculation of the Drell-Yan cross section, a considerably 
softer valence-quark distribution was obtained at high-x. 
This is readily understood. The Drell-Yan cross-section 
factorises into two pieces: one hard and the other soft. 
The soft piece involves the PDF and the hard piece is cal- 
culable in perturbation theory. Adding additional inter- 
actions to the latter, which are important at large- a;; viz., 
soft gluons, provides greater strength in the hard piece 
on the valence-quark domain. Hence a description of the 
data is obtained with a softer PDF. Indeed, the distribu- 
tion obtained thereby matches precisely the expectations 
based on perturbative-QCD and obtained using DSEs. 
This is evident in a comparison between the dash-dot and 
dash-dot-dot curves in Fig. lVIlPl This outcome again 
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Jt-N Drell-Yan, <Q> = 27 GeV" 
Q 2 = 27 GeV 2 , Full BSE 
Q 2 = 27 GeV 2 (r = y , y 5 y.P; q.P=0) 
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FIG. VIII. 3. DSE prediction for the ratio of w-quark distri- 
butions in the kaon and pion 

H,H1. The full Bethe-Salpeter 
amplitude produces the solid curve; a reduced BSE vertex 
produces the dashed curve. The reduced amplitude retains 
only the invariants and amplitudes involving pseudoscalar and 
axial vector Dirac matrices, and ignores dependence on the 
variable q ■ P. These are part of the reductions that occur 
in a pointlike treatm ent of pse udoscalar mesons. The experi- 
mental data is from [1981 , 1 1991 ] . 



emphasises the predictive power and strength of using a 
single internally-consistent, well-constrained framework 
to correlate and unify the description of hadron observ- 
ables. 

The ratio uk/u^ measures the effect of the local 
hadronic environment. In the kaon, the u-quark is bound 
with a heavier partner than in the pion (m s w 25md) 
and this should cause uk{x) to peak at lower- a; than 
Unix). In fact, one finds that the s-quark distribution 
peaks at a value of x which is just 15% larger than that 
of the w-quark. Hence, even though DIS is a high-Q 2 pro- 
cess, constituent-quark- like mass-scales explain this shift: 
M s /M u w 1. 25. The DSE prediction ofii^K @,[zi is 
shown in Fig. IVIII. 31 along with available Drell-Yan data 
[l98l 1 199| | . The parameter- free DSE result agrees well 
with the data. We note that 



U K (0) DGLAP:Q 2 - 



1 



(VTII.6) 



viz, the ratio approaches one under evolution to larger 
resolving scales owing to the increas ingly large popula- 
tion of sea-quarks produced thereby 207]. On the other 
hand, the value at x = 1 is a fixed-point under evolution: 
VQ? > Ql 
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i.e., it is the same at every value of the resolving scale 
Q 2 , and is therefore a persistent probe of nonperturbativc 
dynamics |{|. 



With Ref. [78| a significant milestone was achieved; 
viz., unification of the computation of distribution func- 
tions that arise in analyses of deep inelastic scattering 
with that of numerous other properties o f pseudos calar 
mesons, including meson-meson scattering j208l |209| and 
the successful prediction of electromagnetic clastic and 
transition form factors. The results confirm the large- 
x behavior of distribution functions predicted by the 
QCD parton model; provide a good account of the ir-N 
Drell-Yan data for (x) ; and a parameter- free prediction 
for the ratio uk(x)/u t (x) that agrees with extant data, 
showing a strong environment-dependence of the it-quark 
distribution. The new Drell-Yan experiment running at 
FNAL is capable of validating this comparison, as is the 
COMPASS II experiment at CERN. Such an experiment 
should be done so that complete understanding of QCD's 
Goldstonc modes can be claimed. 



IX. CHARM AND BEAUTY MESONS 
A. CP violation and strong phases 

The past two decades have brought important ad- 
vances in flavour physics and in particular the under- 
standing of weak decays of heavy-mesons. From the first 
observation of a B meson by the CLEO Coll abor ation in 
1981 at the Cornell Electron Storage Ring [210j . to the 
dedicated P-physics facilities at SLAC in California and 
KEK in Japan, much progress has been made. Although 
B physics was the main focus of the Belle and BaBar Col- 
laborations at KEK and SLAC, respectively, and of the 
CDF experiment at Fermilab, considerable efforts have 
also been devoted to studies of P-meson decays, charmo- 
nium and r physics. 

The driving force is the confirmation of the electroweak 
sector of the Standard Model, which has established it- 
self as the foremost paradigm to describe CP violation in 
terms of the Cabibbo-Kobayashi-Maskawa (CKM) mech- 
anism. The aim of experimental data analyses of branch- 
ing fractions, polarisations or CP-violating amplitudes 
in, today, a large variety of decay channels, is determina- 
tion of the precise area of the "CKM triangle" and the 
weak CP violating phase encoded within its angles. 

From a theoretical perspective, heavy mesons can be 
used to test simultaneously all manifestations of the Stan- 
dard Model, namely the interplay between electroweak 
and strong interactions. It is noteworthy that no CP- 
violating amplitude can be generated without relative 
strong phases, i.e. a single strong phase is insufficient 
to produce a CP asymmetry. This can be understood as 
follows: suppose a heavy-meson decays, H — > M, where 
M = M\,M%... denotes the final-state mesons, and that 
the Lagrangian of the Standard Model contributes two 
terms (e.g., two Feynman diagrams) to this process. The 
decay amplitude and its corresponding CP conjugate are, 
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most generally, 

A(H - 
k{H - 



M) 
M) 



\ 2 A 2 e ltp2 



(IX.l) 
(IX.2) 



The weak coupling Aj is a combination of possibly com- 
plex CKM matrix elements and Aje llp denote the strong 
component of the transition amplitude. We emphasise 
that they, too, can have both a real part, or magnitude, 
and a phase, or absorptive part, owing to final rcscatter- 
ing of quarks and mesons. These CP-related intermedi- 
ate states must contribute the same absorptive part to 
the two decays, therefore the strong phases <pi are the 
same in Eqs. (|TXTj) and ([TO]) . The difference of the 
absolute squares, known as direct CP violation, is given 
by 

\A | 2 - |i| 2 = 2A 1 A 2 Im(AiA^) sin(^ - V2 ) . (IX.3) 

It is plain that CP violation, \Si/si\ ^ 1, cannot occur 
if all weak couplings are real or if the decay amplitude 
contains no relative strong phases. Thus, it is crucial to 
evaluate reliably any strong contributions to the decay in 
order to extract weak CKM phases with precision. 



B. Flavourful hadrons 

Heavy mesons also provide an excellent opportunity to 
reassess the nonperturbative features of QCD discussed 
above for light hadrons. A daunting challenge, how- 
ever, is presented by the widely disparate energy scales 
involved in heavy-meson decays. Not surprisingly, fac- 
torisation theorems, which allow for a disentanglement 
of short-distance or hard physics from long-distance or 
soft physics, are a central aspect of heavy-meson phe- 
nomenology. 

Soft matrix elements, commonly expressed as hadronic 
form factors, couplings or decay constants, usually con- 
sist of multiple Green functions between two or more 
physical mass states, with one Green function describ- 
ing the propagation of the heavy quark. These matrix 
elements are formally obtained via the same expedients 
introduced for the flavoured light mesons: Eqs. (|V.3|) to 
(|V.5[) yield the heavy-pseudoscalar's mass formula and 
weak decay constant and, likewise, the related Bethe- 
Salpeter amplitude of a heavy pseudoscalar meson is 
given by Eq. (|V.6|) . 

To this end, the extension of the dynamical drcsscd- 
quark mass- function, M(p 2 ), to hea vy flavours via solu- 
tions of Eq. pV72]) is necessary [H, HI- This is 
depicted in Fig. IIX.1I for different current quark masses 
and the chiral limit. It is clear from this figure that 
whilst DCSB is at the origin of a rapid increase of a 
light quarks's mass-function to the order of several hun- 
dred MeV in the infrared, the effect of dressing the c 
quark is modest and barely noticeable for the b quark, 
whose large current-quark mass almost entirely sup- 
presses momentum-dependent dressing. Thus, M&(p ) 
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FIG. IX.l. Dynamical quark mass function obtained in 
rainbow-ladder truncation for the flavours qt_= u, d, s, c, b and 
in the chiral limit. See, for example, Refs. 0, [68| for details. 



is nearly constant on a large domain. This is true to a 
lesser extent for the charm quark. 

This can also be appreciated from the definition of a 
single quantitative measure, namely the renormalisation- 
point invariant ratio qf := af/M E , where er/ is a 
constituent-quark cr-term fl9| , 



°s : = m f(C) 



dMf 



(IX.4) 



which probes the impact of explicit chiral symme- 
try breaking on the mass function. The Euclidean 
constituent-quark mass is defined via 

(M E ) 2 := {s\s > 0, s = M 2 {s)}. (IX.5) 

For the solutions depicted in Fig. IIX.ll one finds: 



/ 


chiral u, d 


s 


c b 


Mf/GeV 


0.42 0.42 
0.02 


0.56 
0.23 


1.57 4.68 
0.65 0.8 



(IX.6) 



The ratio quantifies the effect of explicit chiral sym- 
metry breaking on the dressed-quark mass-function com- 
pared with the sum of the effects of explicit and dy- 
namical chiral symmetry breaking. Its values are readily 
understood. In the neighbourhood of the chiral limit, 
is small because the magnitude of the light-quark 
constituent-mass owes primarily to DCSB, whilst ap- 
proaches unity for heavy quarks because explicit chiral 
symmetry breaking becomes the dominant source of their 
mass. 



C. Quarkonia 

We first review some aspects of heavy equal-mass, 
, pairs before discussing heavy-light flavoured mesons, 
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in the following section, for which DCSB is decisively 
more relevant. First explorations for mass-symmetric 
QQ states combining a consistent use of the rainbow- 
ladder truncation in the kernels of the gap- and Bethe- 
Salpeter- equation s were undertaken for cc bound states 
in Refs . 134 . 1 2 1 1 j | . with extrapolations to the bb systems 
in Ref. |212| . A full numerical solution for flavour-singlet 
pscudoscalar mesons yielded charmonium and bottonium 
masses and decay constants in very good agreement with 
experimental data; predictions for states w ith exotic 
quantum numbers were also made jl3ll Il33l |. The ef- 
fect of the quark-gluon interaction in the gap equation 
and the vertex-con sisten t Bethe-Salpeter kernel was in- 
vestigated in Ref. [128]. (This is discussed further in 
Sec. UXDl ) 

For a quark flavour Q, a constituent-quark spectrum- 
mass may be defined through 



M, 



Q 



M Q (p 2 
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where Mq(p 2 ) is the renormalisation-point- invariant 
dressed-quark mass-function in Eq. (|IV.2[) . obtained as 
the solution of Eq. (|IV.1[) when a heavy-quark current- 
mass is used via Eq. pV.3[) . As the rcnormalisation-group 
invariant current-quark mass, rhq in Eq. pV.5[) . is in- 
creased, Mq becomes equivalent to the so-called pole- 
mass in non-relativistic QCD (NRQCD). For 0~ quarko- 
nia one can prove [3~ij 



M%® = 2M| 



(IX. 



where ej? is a binding-energy that does not grow with 
Mq. Here Hq denotes the lowest- mass pseudoscalar and 
increasing n labels bound-states of increasing mass. 

Using the renormalisation-group-invariance of the 
product m(( 2 )pH n (( 2 ), Eq. (|V.3|) predicts 
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which establishes an identity between the pseudoscalar 
and pseudovector projections of the meson's Bethc- 
Salpctcr wave function at the origin in configuration 
space. The elements in Eq. pX.9[) are each gauge in- 
variant and renormalisation point independent. In case 
of heavy-ligh t systems it can be shown algebraically 
[H Ell Ell that E q- PXJ]) is realised via 
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which follows from the expansion in 



41 /M S Q 



and 



w hI/ M q of thc integrands in Eqs. ((V73]) to (f\T5|) . Here, 

w® q is the width of the meson's Bethe-Salpeter wave 
function defined as the value of the relative momentum, 
k, whereat the first Chebyshev moment of the amplitude 
E^ 9 falls to one-half of its maximum value. In heavy- 
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FIG. IX. 2. Evolution of the spatial size, rg^, with 1/M|. 

Circles - r®® calculated using the interaction model of 
Ref. [2 1 5f | : solid curve, as described by Eq. (|IX.11I) ; dashed 
curve - linear fit to calculated result in neighbourhood of 
Mj? . Dashed vertical lines mark, from left, and 1/Mj? . 



Note that as Mq -> oo, 



0. 



light-quark. A remarkable feature is that u^ 9 reaches a 
finite nonzero value in the limit Mq — > oo which implies 
that a heavy-light meson always has a nonzero spatial 
extent. 

As shown in Ref. [215| . this is not thc case for QQ 
systems. Figure [TX. 21 depicts the evolution of the heavy- 
heavy meson's spatial size as a function of the inverse 
constituent-quark spectrum-mass, MqQ The curve in 
Fig. IIX.2l is the function 
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t m + -. ; 7m = 0.68, t m = 8.56, 

Aqcd 

(rx.il) 

with Aqcd = 0.234 GcV. One deduces that with increas- 
ing constituent-quark mass, the heavy-heavy system be- 
comes "point-like" in configuration space and hence de- 
localised in momentum space. 

Consequently, the evolution with Mq of an observable 

like the decay constant, f^®, is sensitive to the QQ inter- 
action over a wide range of momentum scales and there- 
fore a useful probe of that interaction. In this connec- 
tion we note that NRQCD predicts the matrix elements 
for various spin states of a given quarkonium system 
are equal up to corrections of order Vq w (w®® /Mq) 2 , 

where k ~ is the typical magnitude of the heavy- 

constituent's three-momentum in the meson's rest frame 
}216j . In this picture, ^ and 1~~ mesons, which differ 
because spins are anti-aligned in the pseudoscalar and 



7 The evolution was computed using thc rcnormalisation-group- 
improved rai nbow -ladder DSE truncation and the interaction 
model in Ref. I215H . As shown, e.g., in Ref. Il28ll , all corrections 
to this truncation vanish in the heavy-heavy limit. 



30 



TABLE IX. 1. Calculated masses and electroweak decay con- 
stants for ground state pseudoscalar and vector heavy-light 
mesons c ompa red with respective experimental data when 
available [168[ |. except for fs s for which the most recent 
num erical value from lattice-QCD is indicated by a dagger 
[21S| (all values are GeV). In the rows labelled "R.L.", the 
heavy quark in the rainbow-ladder truncation is described by 
a constituent-quark mass, which is fit to the lightest pseu- 
doscalar (the D and B mesons). In the rows labelled "fc m i n ", 
the heavy quark is dressed through the DSE with an infrared 
suppression of the gluon momentum as described in the text. 
No such infrared suppression is applied to the dressing of t he 
light quark or the binding kernel. (Adapted from Ref. [217t | 



M Qq 


D D* 


D s D* s 


B B* 


B s B* 


B c B* 


Exp. 


1.86 2.01 


1.97 2.11 


5.28 5.33 


5.37 5.41 


6.29 


R.L. 


1.85 2.04 


1.97 2.17 


5.27 5.32 


5.38 5.42 


6.36 6.44 


&min 


1.88 


1.90 


5.15 


4.75 


5.83 


f Qq 

•'o-(i-) 


D D* 


D a d: 


B B* 


B s B* a 


B c B* 


Exp. 


0.22 


0.29 


0.18 


0.23 1 " 




R.L. 


0.15 0.16 


0.20 0.18 


0.11 0.18 


0.14 0.20 


0.21 0.18 


n>min 


0.26 


0.28 


0.27 


0.16 


0.45 



aligned in the vector, become degenerate in the limit 
Mq — > oo and their leptonic decay constants become 

identical; i.e., f®® = f®®- It is noteworthy, however, 
that Eq. pX.llj) gives v\ w 0.27 and v\ w 0.18. More- 
over, Vq falls only as a^(Mg). Hence a quantitative dis- 
crepancy between fQQ and f®® can conceivably exist 
until rather large quark masses. 



D. Heavy- light mesons 

The mass asymmetry in flavoured Qq mesons leads to a 
diverse array of energy scales, a feature that leads to diffi- 
culties in DSE studies not encountered in the calculation 
of either light or heavy equal-mass systems within the 
rainbow-ladder truncation. For illustration, Table IIX.1I 
unmistakably states that while masses of flavoured pseu- 
doscalar mesons, such as D/ s ) and Br s \ mesons, are in 
good agreement with experimental data, this is n ot so 
for their respective weak decay constants |212l . |217| | . 

More precisely, if the masses and decay constants of 
heavy-light ground state pseudoscalars and vectors in- 
volving a c- or 6-quark are calculated using DSE so- 
lutions for the dressed light-quarks and a constituent- 
quark propagator for the heavy quark, the various me- 
son masses are readily reproduced (see values in the row 
labelled R.L.). This is expected since, as we saw in 
Sec. lIXBl heavy quarks, and in particular the b quark, 
have mass functions whose momentum dependence is 
nearly constant on a large domain. However, the weak 
decay constants obtained from the constituent-like mass 
approximation are 30-50% below the experimental val- 



ues. 

Using fully-dressed quark propagators, both heavy and 
light, in the rainbow-ladder model does not much im- 
prove the situation. Indeed, within the contextof fully- 
dressed propagators, it was numerically shown |212| that 
for small current-quark masses, weak decay constants 
of pseudoscalar and vector heavy-heavy and heavy-light 
mesons increase with the quark mass, yet tend to level 
off between the s- and c-quark mass. This behaviour 
is consistent with the heavy-quark limit; namely, a de- 
crease of the decay constant with increasing meson mass 
like ]%l cx l/^/m% q n , Eq. pX.lOj) . This asymptotic be- 
haviour might occur as low as Q = c for the Qu mesons. 
However, as in the case of constituent-quark propagators, 
fn and f p, are about 20% below their experimental val- 
ues [212| . Since a Qq decay constant depends on the norm 
of the Bethe-Salpetcr amplitude, which in turn depends 
on the derivative of the quark propagators, it might be 
anticipated that the decay constants are more sensitive 
to details of the model than meson masses and are thus 
better indicators of deficiencies in the modelling. 

These observations suggest strongly that the rainbow- 
ladder truncation is not a reliable tool for the study of 
heavy-light mesons the charm-quark region (or larger), 
despite the fact that experimental meson masses are well 
reproduced. In fact, in such systems cancellations, which 
largely mask the effect of dressing the quark-gluon ver- 
tices in light-light mesons, arc blocked by the dressed- 
propagator asymmetry. 

In this context, as discussed in Ref. jl28| . Fig. IIX.3l is 
instructive: it shows the evolution with current-quark 
mass of the difference between the qq meson mass calcu- 
lated in the rainbow-ladder truncation and the value ob- 
tained using a completely resummed dressed-quark-gluon 
vertex in the gap equation and vertex-consistent Bethe- 
Salpeter kernel. With growing current-quark mass, the 
rainbow-ladder truncation provides an increasingly reli- 
able estimate of the meson masses. The accuracy is best 
for pseudoscalar qq mesons; e.g., the absolute difference 
reaches its peak of ~ 60 MeV atm~ 4m s whereat the rel- 
ative error is only 3%. This behaviour in the pseudoscalar 
sector is fundamentally important because it expresses 
Goldstone's theorem: no truncation is valid unless it 
guarantees that the axial-vector Ward-Takahashi iden- 
tity be preserved, something which guarantees a mass- 
less pseudoscalar meson in the chiral limit. Of sign- 
ficance here is the observation that for large current- 
quark masses, say m ~ nib, the contributions from ver- 
tex corrections, which produce nonplanar diagrams in 
the Bethe-Salpeter kernel, are suppressed. Nonetheless, 
Fig. IIX.31 emphasises that the rainbow- ladder truncation 
is not appropriate for charmed mesons; moreover, even 
for heavier mesons, corrections to that truncation ought 
to be included in precision spectroscopic calculations. 

An ad hoc approach to alter th e inf rared sector of the 
model kernel was recently devised [2171 ] . In essence it sup- 
presses infrared contributions from a quark-gluon-vertex 
model originally introduced for light quarks. This intro- 
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E. Flavoured form factors and couplings 
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FIG. IX. 3. Evolution with current-quark mass of the differ- 
ence between the meson mass calculated in the rainbow-ladder 
truncation and the exact value; namely, that obtained using 
the completely resummed dressed-quark-gluon vertex in the 
gap equation and the vertex-consistent Bethe-Salpeter ker- 
nel. The upper panel depicts the absolute error and the lower 
panel, the relative error. Solid lines: vector meson trajecto- 
ries; and dashed-lines; pseudoscalar meson trajectories. The 
dotted vertical lines mark the m u ,d, rn s , m c and m.b current- 
quark masses. 



duces a sharp cut-off, fc m j n , for the gluon momentum in 
the DSE kernel of the heavy quark, while the original 
kernel is kept for the dressing of the light quark and the 
binding. The results are displayed in the rows of Ta- 
ble IIX.1I labelled by "fc m j n " . The decay constant values 
for pseudoscalar mesons (vector mesons were not inves- 
tigated in this scheme) increase dramatically for the D, 
D s and B meson, and are somewhat closer to experiment 
than the values obtained with the constituent-mass ap- 
proximation. The expedient has a somewhat detrimental 
effect on the computed masses. 

To summarise, it has become clear that realistic calcu- 
lations of heavy-light Bethe-Salpeter amplitudes must in- 
volve a truncation that goes beyond rainbow-ladder. The 
exact form of the Bethe-Salpeter equation developed in 
Ref. [97| might provide the advance needed. 



As a corollary of Eq. (|IX.3[) , complex weak phases are 
a necessary but not sufficient condition for CP- violating 
amplitudes. Any decay must involve non-vanishing rel- 
ative strong phases. Their contributions may originate 
at any stage of the heavy-meson decay and can arise 
through: hard radiative corrections to dimension-six and 
higher operators in heavy quark effective theory (HQET) ; 
hadronic transition matrix elements between initial- and 
final-state mesons; as well as final- state interactions be- 
tween daughter hadrons [219H222 [ . T he d evelopment 

of perturbative QCD factorization [223I4227 1 provides a 
means of simplifying this problem by enabling a system- 
atic approximation for a given process in terms of prod- 
ucts of soft and hard matrix elements. 

Transition form factors that characterise the decays 
of S-mesons into light pseudoscalar and vector mesons 
are basic to an understanding of B-meson exclusive non- 
leptonic, semi-leptonic and rare radiative decays. An un- 
derstanding of the entirety of decay processes is essential 
to the reliable determination of CKM matrix elements, 
which is also meant to provide a means of searching for 
non Standard Model effects and CP violation. Consider- 
ing all these factors, it is not surprising that heavy-light 
form factors are the subject of much experimental and 
theoretical scrutiny. 

An immediate consequence of QCD-factorisation the- 
orems for B — > M\M<x decays, mb ^> Aqcd, is the emer- 
gence of two hadronic matrix elements, 



(M 1 M 2 \O i \B) = (Mi|ii|fl)(Af 2 |ia|0) 



r n a n s + 0(A QC D/™.fc) 



, (IX.12) 



where j± and j 2 are bilinear currents. This facto risation 
has been verified to leading order in a s [223l |227| | and in- 
cluding the one- loop correction, a 2 s , to tree-diagram scat- 
tering between the emitted meson and the one cont aining 
the spectator- (light-) quark in the S-mcson [228]. It is 
evident from Eq. (|IX.12[) that higher orders in a s violate 
the factorisation but the corrections can systematically 
be included. The analogy with perturbative factorisa- 
tion for exclusive process es in QCD at large-momentum 
transfer is not accidental [l80j . Of course, semi-leptonic 
decays are much cleaner and, in the absence of inter- 
actions with a second meson, the transition form factor 
{M\\Oi\B) co mp l etely desc ribes the decay amplitude. 

References [H, l213l |229| | compute heavy-to-light tran- 
sition amplitudes, (M(p 2 )\qTjh\H(pi)), in generalised 
impulse approximation, which is the leading-order in the 
systematic, nonperturbative and s ymm etry-preserving 
DSE truncation scheme of Refs. [oil. |62| . The transition 
amplitude is given by 



d 4 k 



f^\k;~p 2 )S q (k+p 2 ) 



X-ipuPi) = trcD . 

J (^r 

x ^i(Pi,P2)S Q (k +p 1 )T H (k;p 1 )S q >(k), 



(IX. 13) 
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where: S(p) are dressed quark propagators, Q = c,b, 
q = q' = u,d,s; and M = S,P,V,A, with the index 
fj, indicating a possible vector structure in the final-state 
BSA. Ti is the interaction vertex whose Lorentz structure 
depends on the operator Oi in the HQET and Th is the 
heavy meson BSA. The trace is over Dirac and colour 
indices. For M = P,V and Tj = 7^(1 — 75), o^g", 
the amplitude !A (pi,P2) can be decomposed into Lorentz 
vectors and tensors as follows (qi = u,d,s,c), 

(P( P 2)\qn„b\B) = F+iq 2 ^ + F-{q 2 )q^ (LX.14) 
lF T (q 2 ) 



B-pi form factors 



(P(p2)\qi(J^aq a b\B) = 



mi + m 2 



{q 2 P fl -q-Pq fl }, 



(IX.15) 



{V{P2^2)\qil^l-lf > )b\B 



mi + VTL2 

{-V P-qA Q (q 2 ) + P^P„A + (q 2 ) (LX.16) 



-q^PvA-iq 2 



which defines the heavy-to-light transition form factors 
F±(q 2 ), F T (q 2 ), A (q 2 ), A ± (q 2 ) and V(q 2 ), and where 
P = pi + P2 , q = Pi — P2 j ^2 i s the polarisation vec- 
tor of the vector meson, and mi, TO2 are the respective 
initial- and final-state meson masses. Similar decomposi- 
tions hold for M = S, A an d oth er four-quark interaction 
operators Tj; see, e.g. Ref. |229j . 

As noted in Scc. lIX Bl the ratio q, indicates that the 
mass function of the b quark varies little on a large 
momentum-squared domain, whilst the alteration for the 
c quark is modest. It is therefore sensible to replace so- 
lutions for the dressed quark propagators from DSEs by 
constituent-mass propagators whose dressed mass is con- 
stant: 



Sn= 



1 



Q=c,b 



Z7 -p + Mq 



Furthermore, in Sec. lIX Dl we discussed obstacles to defin- 
ing a realistic symmetry-preserving truncation scheme for 
the interaction kernel of heavy-light mesons, which has 
so far impeded efforts to build adequate Bethe-Salpeter 
amplitudes for heavy- light ^ and mesons. Cal- 
culations arc for now limited to Bcthe-Salpctcr ampli- 
tude Ansdtze with a single width parameter that is fixed 
by a fit to extant hadronic data. Nonetheless, the form 
factors calculated with such a DSE-based model are ob- 
tained for the entire physical momentum domain without 
any e xtrapolati ons and the chiral limit is directly acces- 
sible [Ilg ill. Predictions for the B ->• tt, B -> K 
and B —> p transition form factors using Eqs. pX.131) . 
PX.14|) . (|LX.15|) and (|IX.16|) are shown in Figs.|IO]and 
ECU 

To make contact with HQET, the heavy-quark limit 
has also been investigated. In this limit, the heavy- 
quark carries all the momentum of the heavy meson, 
which constrains the momentum of the light-quark in the 
Bethe-Salpeter amplitude. The heavy-meson momentum 
is P^ = (Mq + E)vn, where E = M H - Mq, M h is the 
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FIG. IX.4. Results for the form factors F±{q 2 ) in Eq. (|IX.14|) 
(IX 17) an d Fr(q 2 ) in Eg. (|IX.15|l . Top panel, B — > n; and bottom 



panel, B — >■ K: solid F+, dot-dashed F-, dashed Ft 



mass of the heavy meson and is a time-like unit vector, 
v = — 1. A heavy quark has four- momentum P^ + fc M , 
where the residual momentum k sa Aqcd is associated 
with motion of the light-quark within the bound-state. 
In the limit Mq 00, the expansion of the propagator 
in Eq. PX.17P in powers of 1/ Mq yields 



S Q {k + P) 



1 1 — 17 ■ v 

2 k-v-E 



Mq Mq 



(IX.18) 



given that \k\/MQ <C 1. If we employ a two- 
covariant Ansatz for the heavy pseudoscalar meson's 
Bethe-Salpeter amplitude; viz., 

T Hg _(k;P)= l5 (l-± n -v)^-v(k 2 ) , (IX.19) 

where 9{h is the canonical normalisation constant, then 
the weak decay constant, Eq. (|V.3[) . becomes, with z = 
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B-rho+gamma form factors 
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FIG. IX.5. Results for the B -s> p form factors A (q 2 ), A±(q 2 ) 
and F(g 2 ) in Eg. (|IX.16|) : so/id curve, Aq; dashed, A+; dot- 
dashed, A-\ and dotted V. 
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Here we have introduced the My-independent canoni- 
cal normalisation, Kt, in order to make explicit the pre- 
viously mentioned relation, Eg. pX.10[) : /jj _ \/Mh = 
constant. 

Similarly, at leading order in 1/Mq, the B — > P tran- 
sition form factors in Eq. (|IX.14[) are 



2 1Mp±Mb 
F±(9)= 2Vi™ eM; 



(IX.21) 



i.e. the form factors depend on a single universal function 
(the so-called Isgur-Wise function): 



: 32tt- 



1 f 00 
dr— J du<p 2 (z w ) 



w 



4(z) 



(IX.22) 



with W = 1 + 2r(l - t)(w - 1), zv/ = « - 2Ey/u/W 
and w = -i>b • t'p = (M| + Mp - q 2 )/(2M B M P ). Ow- 
ing to the canonical normalisation of the Bcthc-Salpctcr 
amplitudes, £(w = 1) = 1. 

The DSEs plainly reproduce the results of heavy-quark 
symmetry and, in addition, provide a means by which 
to examine the fidelity of the formulae obtained in the 
heavy-quark limit. In other words: one may pose and 
answer the question: under which conditions is heavy- 
quark symmetry applicable? Based on the analysis in 
Ref. fo8j , which presents a unified and uniformly accurate 



description of a broad range of light- and heavy-meson 
observables, one concludes that corrections to the heavy- 
quark symmetry limit of < 30% are encountered in b — >• c 
transitions but that these corrections can be as large as 
a factor of 2 in c — > d transitions. 

Radiative and strong decay amplitudes have also 
been calculated in the impuls e approx imation, and their 
heavy-quark limits exposed [68l . l230j . The amplitude 
for the decay -> H(p 2 )~f(k), with p\ = -M%,, 



2 



-Mfj and k 2 = 0, is given by 

x[e Q MQ( Pl ,p 2 ) + e qf M*i( Pl , P2 )] , (IX.23) 

where £Q, qs is the fractional charge of the active quark in 
units of the positron charge. The sum indicates that the 
decay occurs via a spin-flip transition by either the heavy 
or light quark. The hadronic tensors in Eq. (|IX.23|) are: 



M%,(pi,P2) = tr CD 



d A £ 



r H (£;- P2 )s Q (e 2 ) 



M^( Pl ,p 2 ) 



(2tt)4 

xir^ 2 ,4)s Q (4)rf (t;Pi)s qf (£) 

(IX.24) 

d A e 



tr 



CD 



t H (i;-p 2 )S Q {h) 



(2tt)' 



x (£; Pl )S qf (£) iTl(£, £ + k)S qf {£ + k) , 

(IX.25) 

where Q = c,b, qj = u,s and £ 1>2 — t+Pi,2- Fl(£i,£ 2 ) 
is the dressed-quark-photon vertex, for which the BC 
Ansatz, Eq. (|VI.8|) . was used because it satisfies the vec- 
tor Ward-Takahashi identity and thus ensures current 
conservation. In numerical calculations it is helpful that, 
in case of a heavy constituent, T(,(£i,£ 2 ) = ~f v owing to 
Eq. (|rX~T7j> : A Q = 1 and B Q = M Q . 

Likewise, in the impulse approximation, the strong de- 
cay H*(pi) — > H(p 2 )Tr(q) is described by the invariant 
amplitude, 



H'Htt 
fit. 



(IX.26) 



The dimensionless coupling gn*Hixi which can in prin- 
ciple be determined via the decay width Th'Htt, is re- 
lated to the putative universal strong coupling, g, be- 
tween heavy-light-vector and -p seud oscalar mesons in a 
heavy-meson chiral Lagrangian [231]. At tree level, the 
couplings gii'Hix and 5 arc related as, 



gn-Hir = 



y/m H m H * 



-g- 



(IX.27) 



Practically, the matrix element in Eq. (|IX.26[) describes 
the physical process D* — > Dir, with both the final pseu- 
doscalar mesons on-shell. It also serves to compute the 
unphysical soft-pion emission amplitude B* —> Bit in the 
chiral limit (raj. — raj < m n ), which defines gB'B-n- A 



34 



comparison between these two couplings is an indication 
of the degree to which notions of heavy-quark symme- 
try can be applied in the charm sector. This was done 
in Rcf. 2321. wherein it was demonstrated that the dif- 



ference in either extracting g from gr>*Dix or qb'B-k is 
material: 



9d*d-k = 15.8 

ffB*B7T = 30.0 



g = 0.53 , 
g = 0.37 . 



(IX.28) 
(IX.29) 



Here again, the results emphasise that when the c-quark 
is a system's heaviest constituent, Aqcd /m c -corrections 
are not well controlled. Since the DSE approach lends it- 
self to the calculation of other amplitudes involving light 
and heavy mesons; e.g., couplings of phcnomenological 
inter est for charmonium production and D-mesic nuclei 
[2031 ] . it should help to overcome the limitations of an 
expansion in 1/M C . 



X. DESCRIBING BARYONS AND MESONS 
SIMULTANEOUSLY 

While a symmetry-preserving description of mesons is 
essential, it is only part of the physics that nonpertur- 
bative QCD must describe because Nature also presents 
us with baryons: light-quarks in three-particle compos- 
ites. An explanation of the spectrum of baryons and the 
nature of interactions between them is basic to under- 
standing the Standard Model. The present and planned 
experimental programmes at JLab, and other facilities 
worldwide, are critical elements in this effort. 

No approach to QCD is comprehensive if it can- 
not provide a unified explanation of both mesons and 
baryons. We have explained that DCSB is a key- 
stone of the Standard Model, which is evident in the 
momentum-dependence of the dressed-quark mass func- 
tion - Fig. III. ll it is just as important to baryons as 
it is to mesons. Since constituent-quark-like models can- 
not incorporate the momentum-dependent dressed-quark 
mass-function, they are not a viable tool for use in this 
programme. The DSEs furnish the only extant contin- 
uum framework that can simultaneously connect both 
meson and baryon obscrvables with this basic feature of 
QCD, having provided, e.g., a direct correlation of me- 
son and baryon properties via a single interaction kernel, 
which preserves QCD's one-loop renormalisation group 
behaviour and can systematically be improved. This is 
evident in the preceding secti ons and their combination 

with Refs. [m EH El mi . 

In order to illustrate this programme, we will re- 
view the computation of baryon masses, and nucleon 
elastic and transition electromagnetic form factors. In 
this connection, we recall that, for a structureless or 
simply-structured fermion Fi(Q 2 ) = 1 and F 2 (Q 2 ) 
I) [see after Eg. (IVLTl) ]. so that G E {Q 2 ) = G M {Q 2 ) 
and the distribution of charge and magnetisation is 
identical within that composite. This was believed to 



be the case for the proton until 1999. In that year, 
enabled by the high luminosity of the accelerator at 
JLab, a new meth od w as employed to measure the ratio 
G P E (Q 2 )/G P M (Q 2 ) (234|. The result astonished the com- 
munity: whilst G P E (Q 2 )/G P M (Q 2 ) M 1 for Q 2 < lGeV 2 , 
G P E (Q 2 ) / G P M (Q 2 ) is a rapidly decreasing function for 
Q 2 > lGeV 2 (see Fig.|XT2]). As one of the highlights 
herein, we will indicate how this may be understood. 



A. Faddeev equation 

In quantum field theory a baryon appears as a pole in 
a six-point quark Green function. The residue is propor- 
tional to the baryon's Faddeev amplitude, which is ob- 
tained from a Poincare covariant Faddeev equation that 
sums all possible exchanges and interactions that can 
take place between three dress ed-q uarks. A tractable 
Faddeev equation for baryons [235j is founded on the 
observation that an interaction which describes colour- 
singlet mesons also generates nonpointlike quark-quark 
(diquark) correlations in the colour-3 (antitriplet) chan- 
nel [236| . The dominant correlations for ground state 
octet and decuplet baryons are scalar (0 + ) and axial- 
vector (1 + ) diquarks because, for example, the associ- 
ated mas s-scales are smaller than the baryons' masses 
[l29ll237| and their parity matches that of these baryons. 
It follows that only they need be retained in approximat- 
ing the quark-quark scatteri ng matrix wh ich appears as 
part of the Faddeev equation jH, 1 139L Il64j . On the other 
hand, pscudoscalar (0~) and vector (1 _ ) diquarks domi- 
nate in the pari ty-pa rtners of ground state octet and de- 
cuplet baryons |l64j . The DSE approach treats mesons 
and baryons on the same footing and, in particular, en- 
ables the impact of DCSB to be expressed in the predic- 
tion of baryon properties. 

In quantum field theory a baryon appears as a pole in 
a six-point quark Green function. The residue is propor- 
tional to the baryon's Faddeev amplitude, which is ob- 
tained from a Poincare covariant Faddeev equation that 
sums all possible exchanges and interactions that can 
take place between three dress ed-q uarks. A tractable 
Faddeev equation for baryons [235| is founded on the 
observation that an interaction which describes colour- 
singlet mesons also generates nonpointlike quark-quark 
(diq uark ) correlations in the colour-3 (antitriplet) chan- 
nel |236| . The dominant correlations for ground state 
octet and decuplet baryons are scalar (0 + ) and axial- 
vector (1 + ) diquarks because, for example, the associ- 
ated mas s-scales are smaller than the baryons' masses 
[l29ll237| and their parity matches that of these baryons. 
It follows that only they need be retained in approximat- 
ing the quark-quark scatteri ng matrix wh ich appears as 
part of the Faddeev equation [H, Il39l Il64| . On the other 
hand, pseudoscalar (0~) and vector (1~) diquarks domi- 
nate in the pari ty-pa rtners of ground state octet and de- 
cuplet baryons |l64j . The DSE approach treats mesons 
and baryons on the same footing and, in particular, en- 
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quark exchanged between the diquarks as 



FIG. X.l. Diagrammatic representation of a Poincare covari- 
ant Faddeev equation for a baryon. is the Faddeev ampli- 
tude for a baryon of total momentum P = p q +Pd- It expresses 
the relative momentum correlation between the dressed-quark 
and -diquarks within the baryon. The shaded region demar- 
cates the kernel of the Faddeev equation, in which: the single 
line denotes the dressed-quark propagator, F is the diquark 
Bethe-Salpeter amplitude; and the double line is the diquark 
propagator. 



ables the impact of DCSB to be expressed in the predic- 
tion of baryon properties. 

It is important to appreciate that diquark s do not ap- 
pear in the strong interaction spectrum [62|, Il28| . How- 
ever, the attraction between quarks in this channel jus- 
tifies a picture of baryons in which two quarks within 
a baryon are always correlated as a colour-3 diquark 
pseudopartiele, and binding is effected by the iterated 
exchange of roles between the bystander and diquark- 
participant quarks. Here it is important to emphasise 
strongly th at Q CD supports nonpointlike diquark corre- 
lations j77l . l238| . as we shall see in Sec. lX B[ Hence models 
that employ pointlike diquark degrees of freedom have 
little connection with QCD. 

The Faddeev equation, illustrated in Fig. IX. 11 is a lin- 
ear homogeneous matrix equation, in many respects sim- 
ilar to a Bethe-Salpeter equation. Its solution is the 
nucleon's Poincare-covariant Faddeev amplitude, which 
describes quark-diquark relative motion within the nu- 
cleon. The composite nature of the diquark correlations, 
and their dynamical breakup and reformation through 
dressed-quark exchange, pictured in Fig. IX.ll is crucial 
to maintaining fermion statistics for the nucleon bound- 
state. Furthermore, owing to the critical importance of 
both scalar- and axial-vector-diquark correlations, the 
nucleon's rest- frame wave- function possesses S-, P- and 
D-wave correlations; i.e., a nucleon should a priori be ex- 
pected to contain significant dressed-quark orbital angu- 
lar momentum. This is verified in Ref. [36j . which shows 
that in the nucleon's rest frame just 37% of the total 
spin of the nucleon is contained within components of 
the Faddeev amplitude which possess zero quark orbital 
angular momentum. 

It is worthwhile expressing one Faddeev equation con- 
cretely and we choose that for the A-resonance, a J = 
3/2 bound-state of three valence light-quarks. It is possi- 
ble to obtain a simple yet extremely informative equation 
if one employs the interaction presented in Eq. (|VII.2[) 
and an additional simplification; i.e., representing the 



S T (k) 



M 



(X.l) 



where <?a = 1-56 [1641 ] . This is a variant of the so-called 
"sta tic a pproximation," which itself was introduced in 
Ref. [2391 ] and has subsequently been used in studying 
a range of nucleon properties ]240]. In combination 
with diquark correlations generated by Eq. ([VII.2[) , whose 
Bethe-Salpeter amplitudes are momentum-independent, 
Eq. (jX.l[) generates Faddeev equation kernels which 
themselves are independent of the external relative mo- 
mentum variable. The dramatic simplifications which 
this produces are the merit of Eq. ([X.l I) . 



Following this through, one derives 164 



1 = 8 si^W I dH 



Mm r 1 



Mm^J (2n fJ Q 
{m 2 qqi+ + (1 - a) 2 m 2 A )(am A + M) 
[I' 2 + cta(q;, M, m q q i+ , m A )] 2 

9 a 



(X.2) 



7 2 E 2 



M ml qi+ 2^ Jo 



da {m 2 qqi+ + (1 - afm\) 



x(am A + M)C™ {<TA(a, M,m qqi+ ,m A )), (X.3) 

where Eqq l+ is the canonically-normalised axial-vector 
diquark Bethc-Salpctcr amplitude, rn qqi+ is the com- 
puted mass for this correlation, 



o-a(oi, M, m qqi+ , tua) 
= (1 -a)M 2 4 



a(l-a)m 2 A , (X.4) 



C l (z) = C\{z)/z, w ith cj u Q) = -z(d/dz)c m {z) and 
C m {z) defined in Eq. ([VII. 4[) . Equation (jX.3[) is an eigen- 
value problem whose solution yields the mass for the 
dressed-quark-core of the A-resonance, tua- It is one 
dimensional because only the axial-vector diquark cor- 
relation contributes to the structure of the A. (The 
nucleon, on the other hand, is constituted from scalar- 
and axial- vector-diquarks and presents a five-dimensional 
eigenvalue problem.) 

With experience, one can look at Eq. (|X.3[) and see that 
increasing the current-quark mass will boost the mass of 
the bound-st ate. This is just one of the Faddeev equa- 
tions in Re f. [164J | . In fact, building on lessons from me- 
son studies [121 ] . a unified spectrum of u, d-quark hadrons 
was obtained therein using the symmetry-preserving reg- 
ularization of a vector x vector contact interac tion that 
we have briefly described herein. Reference [l64| re- 
ports a study that simultaneously correlates the masses 
of meson and baryon ground- and excited-states within 
a single framework. In comparison with relevant quan- 
tities, the computation produces rms=13%, where rms 
is the root- mean-square-relative-error/degree-of freedom. 
As evident in Fig. IX.2l the prediction uniformly overesti- 
mates the PDG values of meson and baryon masses [168| . 
Given that the employed truncation deliberately omitted 



36 




Location of the 3-dressed-quark core 
predicted by the DSE study 
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FIG. X.2. Comparison between DSE-computed ha dron 
masses (filled circles) and: bare baryon masses from Ref. [24 ll ] 
(filled diam onds) and Ref. [2421 ] (filled triangles); and exper- 
iment 168], filled-squares. For the coupled-channels models 
a symbol at the lower extremity indicates that no associated 
state is found in the analysis, whilst a symbol at the upper 
extremity indicates that the analysis reports a dynamically- 
generated resonance with no corresponding bare-baryon state. 
In connection with f2-baryons the open-circles represent a 
shift downward in the computed results by 100 MeV. This 
is an estimate of the effect produced by pseudoscalar-meson 
loop corrections in A-like systems at a s-quark current-mass. 



meson-cloud effects in the Faddeev kernel, this is a good 
outcome, since inclusion of such contributions acts to re- 
duce the computed masses. 

Following this line of reasoning, a striking result is 
agre ement between the DSE-computed baryon masses 
[l64| and the bare masses employed in modem coupled- 
channels models of pi on-nu cleon reactions [24ll . |242| . see 
Fig. lX.2l and also Ref. [l43j . T he Ro per resonance is very 
interesting. The DSE study [1641 ] produces an excita- 
tion of the nucleon at 1.82 ± 0.07 GeV. This state is 
predominantly a radial excitation of the quark-diquark 
system, with both the scalar- and axial-vector diquark 
correlations in their ground state. Its predicted mass 
lies preci sely at the value determined in the analysis of 
Ref. [24l| . This is significant because for almost 50 years 
the "Roper resonance" has defied understanding. Dis- 
covered in 1963, it appears to be an exact copy of the 
proton except that its mass is 50% greater. The mass was 
the problem: hitherto it could not be explained by any 

symmetry-preserving QCD-based tool. T hat has now 

changed. Combined, see Fig. lX.3[ Rcfs. |l64l24lj demon- 
strate that the Roper resonance is indeed the proton's 
first radial excitation, and that its mass is far lighter than 
normal for such an excitation because the Roper obscures 
its dressed-quark-core with a dense cloud of pions and 
other mesons. Such feedback between QCD-based theory 
and reaction models is critical now and for the foreseeable 
future, especially since analyses of experimental data on 
nucleon-resonance electrocouplings suggest strongly that 
this structure is typical; i.e., most low- lying A r *-states 
can best be understood as an i nternal qu ark-core dressed 
additionally by a meso n clo ud [243l . |244| | . 

Additional analysis [143| suggests a fascinating new 
feature of the Roper. To elucidate, we focus first on the 
nucleon, whose Faddeev amplitude describes a ground- 
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FIG. X.3. The Excited Baryon Analysis Center (EBAC) ex- 
amined the Pi i -channel and found that the two poles associ- 
ated with the Roper resonance and the next higher resonance 
were all associated with the same seed dressed-quark state. 
Coupling to the continuum of meson-baryon final states in- 
duces multiple observed resonances from the same bare state. 
In EBAC's analysis, all PDG-identified resonances were found 
to consist of a cor e sta te plus meson-baryon components. 
(Adapted from Ref. [HJ].) 



state that is dominated by its scalar diquark compo- 
nent (78%). The axial- vector component is significantly 
smaller but nevertheless important. This heavy weight- 
ing of the scalar diquark component persists in solutions 
obtained with more sophisticated Faddeev equation ker- 
nels (sec, e.g., Table 2 in Ref. [33|). From a perspective 
provided by the nucleon's parity partner and the radial 
excitation of that state, in whi ch t he scalar and axial- 
vector diquark probabilities are [245j 51%-49% and 43%- 
57%, respectively, the scalar diquark component of the 
ground-state nucleon actually appears to be unnaturally 
large. 

One can nevertheless understand the structure of the 
nucleon. As with so much else, the composition of the 
nucleon is intimately connected with dynamical chiral 
symmetry breaking. In a two-color version of QCD, 
the s calar diquark is a Goldstone mode, just like the 
pion [l25| . (Thi s is a long- known result of Pauli-Gursey 
symmetry [246l |247| .) A memory of this persists in 
the three-color theory and is evident in many ways. 
Amongst them, through a large value of the canoni- 
cally normalized Bcthc-Salpcter amplitude and hence a 
strong quark+quark— diquark coupling within the nu- 
cleon. (A qualitatively identical effect explains the large 
value of the ttN coupling constant.) There is no such en- 
hancement mechanism associated with the axial-vector 
diquark. Therefore the scalar diquark dominates the nu- 
cleon. 

With the Faddeev equation treatment described 
herein, the effect on the Roper is dramatic: orthogonality 
of the ground- and excited-states forces the Roper to be 
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FIG. X.4. Vertex which ensures a conserved current for on- 
shell nucleons described by the Faddeev amplitudes, ob- 
tained from the equation depicted in Fig. IX.ll The photon 
probing the nucleon introduces momentum Q and is repre- 
sented by the wiggly line; the single line represents S(p), the 
dressed-quark propagator, and the double line, the diquark 
propagator; and F is the diquark Bethe-Salpeter amplitude. 
The vertex in Diagram 4 (second row, right column) repre- 
sents an electromagnetically induced transition between the 
scalar and axial- vector diquarks; and the vertices in row 3 
are "seagull terms," which appear as partners to Diagram 3 
and arise because binding in the nucleons' Faddeev equa- 
tions is, in general, effected by the exchange of a momentum- 
carryi ng d ressed-quark between nonpointlike diquark correla- 
tions [24811. 



constituted almost entirely (81%) from the axial-vector 
diquark correlation. It is important to check whether 
this outcome survives with a Faddeev equation kernel 
built from a momentum-dependent interaction. 

With masses and Faddeev amplitudes in hand, it is 
possible to compute baryon electromagnetic form factors. 
For the nucleon, studies of the Faddeev equation exist 
that are based on the one-loop renormalisation-group- 
improvcd interactio n th at was used efficaciously in the 
study of mesons [H, 1 1 3Q| . These studies retain the scalar 
and axial- vector diquark correlations, for the reasons ex- 
plained in Sec. IX Al 

A nonpointlike composite nucleon must interact with 
the photon via a sophisticated current, whose form is con- 
strained by vector Ward-Takahashi identities. This con- 
tinues a thread that pervades these notes. For the bound- 
state described by the Faddeev equ ation in Fig. IX. II 
that current is described in Ref. [248| and depicted in 
Fig. IX. 41 Diagrams 4-6 represent eight-dimensional inte- 
grals, which can be evaluated using Monte-Carlo tech- 
niques. 

The pattern of the computation should now be clear. 
In principle, one specifies an interaction kernel for the gap 
equation and solves for the dressed-quark propagator. 
This completes the specification of the Bcthc-Salpctcr 
kernel, so one may compute the masses and amplitudes 
for the diquark correlations. With all these quantities 



determined, the Faddeev equation is defined and can be 
solved for a baryon's mass and amplitude. In combina- 
tion with the current, Fig. IX.4l it is a straightforward 
numerical task to compute the form factors. That is, at 
least, once one has information about the diquark elec- 
tromagnetic couplings which appear in the current. 



Diquark and meson elastic and transition Form 
Factors 



p-meson elastic form factors 
This is an ideal place to detail the properties of 
nonpointlike diquarks, At leading-order in a symme- 
try preserving truncation of the DSEs jfjll . |fj2| . simple 
changes in the equations describing tt- and p mesons 
yield expressions that provide detailed information about 
the s calar and axial- vector d iquarks; e.g., their masses 



164 12361 12371 . 12491 1250| , and electromagnetic elas- 
tic }238| and transition form factors, which are critical 
elements in the computation of a baryon's kindred prop- 
erties via the current in Fig. IX. 41 It is therefore natural to 
elucidate concurrently the properties of 7r- and p-mesons 
and those of the scalar and axial- vector diquark correla- 
tions because it opens the way to a unified, symmetry- 
preserving explanation of meson and baryon properties 
as they are predicted by a single interactio n. The po ten- 
tial of this approach is apparent in Rcfs. |l39l Il7l| but 
it has yet to be fully realised. For the present the best 
connection is provided by the less rigorous approach of 
Ref. HH , which uses more parameters to express features 
of QCD but also predicts a nd describ es simultaneously a 
larger array of phenomena |25ll4253| . 

A wide-ranging study of u/<i-quark meson and diquark 
form factors, using the contact interaction explained in 
Sec. lVIII is described in Ref. [77j . Here we recapitu- 
late on some interesting points that are relevant to the 
study of baryons. To begin, we note that in analogy 
with Eq. (|VIL5[) . one has the following general form of 
the Bethe-Salpeter amplitude for the p-meson obtained 
with a symmetry-preserving regularisation of a contact- 
interaction: 

= llE p {P) + -L<j^P„F p (P) . (X.5) 

However, as a peculiar consequence of the rainbow-ladder 
truncation: 



„ , „. ladder 

F P (P) = 0. 



(X.6) 



One has F p (P) ^ in any symmetry-preserving trunca- 
tion that goes beyond this leading-order [62j. The acci- 
dent expressed in Eq. (|X.6j) has material consequences. 

Owing to the connection between properties of vector 
mesons and axial- vector diquarks, a computation of the 
p-meson elastic electromagnetic form factors is an im- 
portant step along the way to diquark form factors. The 
' p-meson has three elastic form factors and we 
follow Ref. [351 ] in defining them. Denoting the incoming 
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photon momentum by Q, and the incoming and outgoing 
p-meson momenta by p % = K — Q/2 and p* = K + Q/2, 
then K ■ Q = 0, K 2 + Q 2 /4 = — m 2 and the p-j-p vertex 
can be expressed: 



k Xtllv {K,Q) 
Tl tlv (K > Q) = 



Y. T l^Q) F] (Q 2 ) 



"2m 2 , 



Q: 



■pl 



2m 2 



(X.7) 
(X.8) 

(X.9) 



Kx 
m 2 



"2m 2 



Qu+P 



2m 2 



(X.10) 



where fP^ip) = S pv — p^p v /p 2 . A symmetry-preserving 
regularisation scheme is essential so that the following 
Ward-Takahashi identities are preserved throughout the 
analysis: 



QxA x , flu (K,Q) = 0, 



(X.11) 



P^Ax.MK, Q) = = plAx t „ v (K, Q) ■ (X.12) 

The electric, magnetic and quadrupole form factors are 
constructed as follows: 



G E (Q 2 ) = F 1 (Q 2 ) + -r ] G Q (Q 2 ), 



(X.13) 

G M (Q 2 ) = -F 2 (Q 2 ), (X.14) 
Gq{Q 2 ) = Fx(Q 2 ) + F 2 {Q 2 ) + + F 3 (Q 2 ) , (X.15) 



where r\ = Q 2 /[4m 2 ]. In the limit Q 2 — > 0, these form 
factors define the charge, and magnetic and quadrupole 
moments of the p- meson; viz., 



G£(Q 2 = 0) = 1, 



(X.16) 



G P M {Q 2 = 0) = p p , G P Q (Q 2 =0) = Q p . (X.17) 



It is readily seen that Eq. (|X.16[) is a symmetry con- 
straint. One has G E (Q 2 = 0) = Fi(Q 2 = 0) and 



A(K,Q) 



CT-s-0 



2K X (K) <bZ v (K) F, (0) . (X.18) 



Using the vector WTI for the quark-photon vertex and 
the explicit form for Pr{Q 2 ), Fig. IVII.Tl this becomes 

Jf A (0) 



N c E 2 tT B 



?4 c\ 

1 Hv—S(e + K)i ltl S(£). (X.19) 



(2tt 



dK x 



The right-hand-side (rhs) is simply the analogue of 
Eq. (|VII.8[) for the rainbow-ladder vector meson. Hence, 
when the Bethe-Salpeter amplitude is canonically- 
normalised and so long as one employs a symmetry- 
preserving regularisation procedure, the rhs is equal to 
K x Tl v (K) and thus Fi(0) = 1. 



TABLE X.l. Row 1 : Form factor radii (in fm), and magnetic 
and quadrupole moments for the p-meson, G P M (Q 2 = 0) and 
Gq(Q 2 = 0) respectively, computed with the parameter val- 
ues in Table IVII.ll For a structureless vector meson, p = 2 
and a = -1. Experimentally, rv = 0.672 ± 0.008 fm [lol l. 
Row 2: Results for the scalar and axial- vector diquark corre- 
lations. Here the magnetic and quadrupole moments should 
be multiplied by the relevant charge factor; viz., e^ uu y = |, 
e{ud) = | and e{dd} = — § ■ 



rv 


' p 


r M 
' p 


' p 


Hp 




0.45 


0.56 


0.51 


0.51 


2.11 


-0.85 


r o+ 


rf + 




rf + 


Pi+ 


Qi+ 


0.49 


0.55 


0.51 


0.51 


2.13 


-0.81 



Table IX. II reports computed form factor radii, and 
the magnetic and quadrupole moments. An interpre- 
tation of the ratio r w /r p — 0.80 determined from the 
Table is complicated by the fact that rainbow-ladder 
truncation was used consistently; but in this case alone 
F p (P) = 0, whereas F^(P) always and F p (P) ^ in 
all other truncations. It is useful, therefore, to observe 
that ?v = 0.51 fm if one artificially sets F 7r {P) = 0, in 
which case r^/vp = 0.92. Moreover, the DSE computa- 



tion in Ref. [35|, which employs a QCD-based interaction, 
produces r^/vp — 0.90; and in combin ation , the more 
phenomenological DSE studies of Refs. |202l |254| | yield 
r n /r p = 0.92. 

The computed p-meson electric form factor is plotted 
in Fig. IX.5l It displays a zero at Q 2 = 5.0 GcV 2 and 
remains negative thereafter. Given that the deuteron 
is a weakly-bound J = 1 system, constituted from two 
ferm ions, and its electric form factor possesses a zero 
[2551 ], it is unsurprising that G P E (Q 2 ) exhibits a zero. It 
is notable in addition that the dcutcron's zero is located 
at zg := sjQ 2 = 0.8 GeV, so that 



(X.20) 



where td is the deuteron's radius. An interpolation valid 



on Q 2 e [-m 2 , 10 GeV 2 ] is 
G p E (Q 2 ) 



,2\ interpolation 



1 -0.20Q 2 



1 + 1.15Q 2 -0.013Q 4 



(X.21) 



Figure IX. 51 also depicts the magnetic and quadrupole 
form factors of the p-meson, both normalised by their 
values Q 2 = 0. Notably, neither of these two form fac- 



tors change sign: for Q 2 > 



G P M (Q ) is positive 



definite and G P E (Q 2 ) is negative definite. Furthermore, 
over this entire domain of Q 2 , these form factors exhibit 
a very similar (5 2 -depcndcncc, which is made especially 
apparent via the dotted-curve in Fig. IX.5l Interpolations 
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FIG. X.5. Solid curve - p-meson electric form factor, G P E {Q 2 ), 
which exhibits a zero at Q 2 = 5.0 GeV 2 . (It is notable 



FIG. X.6. Solid curve - the full result for G' r7P (Q 2 ); and 



that 1 



\n = for Q 2 



6mi 



5.2 GeV 2 .) The dashed dashed curve - G 7T1P {Q 2 ) obtained with F n (P) = 0. Experi- 



ence, G P M (Q )/V P , and dot-dashed curve, Gq(Q )/Qp, are 
almost indistinguishable, as emphasised by the dotted curve, 
[G P AI (Q 2 ) / h p ]/[Gq(Q 2 ) / Qp]. The charge radii, and magnetic 
and quadrupole moments are given in Table IX. II NB. All 
form factors exhibit a pole at Q 2 
photon vertex is dressed as described in Sec. lVII 



-m p because the quark- 



(X.22) 



valid on Q 2 e [-m 2 p , 10 GeV 2 ] arc 

p 2 interp. 2.11 + 0.021 Q 2 

mW ' 1 + 1.15 Q 2 - 0.015 Q 4 

nP irfl\ int — p ' 0.85 + 0.038 Q 2 

G Q {Q > ~ "1 + 1.17Q2 + 0.014 Q*- (X - 23j 

The similar momentum-dependence of G P M and Gq re- 
calls a prediction in Ref. }256| ; namely, 



G £ (Q 2 ) : Gm(Q 2 ) : G Q (Q 2 ) Q2 =°° 1 - ^ : 2 : -1 

(X.24) 

in theories with a vector- vector interaction mediated via 
bosons propagating as 1/k 2 at large-fc 2 . The computed 
ratio Tm/q : = G p m (Q 2 )/Gq(Q 2 ) conforms approximately 
with this prediction on a large domain of Q 2 ; e.g., 



Q 2 10 10 2 10 3 
r M/Q -2.48 -2.54 -2.38 -2.17 



(X.25) 



However, at Q 2 = 10 4 GeV 2 , r M /Q = -1-28. Moreover, 
the remaining two ratios are always in conflict with the 
prediction; and closer inspection reveals that even the 
apparent agreement for G p m (Q 2 )/Gq(Q 2 ) is accidental, 
since Eqs. (|X.24[) are true if, and only if, 

Fi(Q 2 ) : F 2 (Q 2 ) : Q 2 F 3 (Q 2 ) Q ^°° 1 : -2 : ; (X.26) 

and none of these predictions are satisfied in this compu- 
tation. 

The mismatch originates, of course, with Eq. (|VII.2|) 
and the concomitant need for a rcgularisation procedure 
in which the ultraviolet cutoff plays a dynamical role. If 



mentally [16S| ] . the pa rtial width for p + —¥ % +r y is 68 ± 7keV, 
which corresponds to [1861 ] g^ lp = (0.74±0.05) m p . This is in 



one carefully removes A uv -> oo, Eqs. (IX3Gl) arc recov- 
ered but at the cost of a logarithmic divergence in the 
individual form factors. One concludes therefore that a 
vector-vector contact interaction cannot reasonably be 
regularised in a manner consistent with Eq. (jX.24l) . 

In closing this subsection it is worth reiterating that it 
is only in the rainbow-ladder truncation that F P (P) = 0. 
Therefore in connection with the p-meson's form factors, 
material changes should be anticipated when proceeding 
beyond this leading-order truncation. 

p-ir transition form factor 
This transition is closely related to the 7*777 transi- 
tion form factor, whose behaviour in connection with 
Eq. (|VII.2[) was analysed in Ref. [76[. The corresponding 
interaction vertex defines a single form factor; viz., 

TZd P {k u k 2 ) = SjnL e p „ afS k la k 20 G^ P (Q 2 ) , (X.27) 
m p 

where k\ = Q 2 , k\ = —rn 2 . The coupling constant, g nip , 
is defined such that G^ P (Q 2 = 0) = 1. 

The contact-interaction result for the form factor [77j 
is depicted in Fig. IX. 61 Naturally, because the quark- 
photon vertex is dressed (see Fig. IVII.Tj) . the transition 
form factor exhibits a pole at Q 2 = —m 2 , which is not dis- 
played. An interpolation valid on Q 2 G [—m 2 p , 10 GeV 2 ] 
is 



r ^ P(n 2, intop. 1 + 0.37Q 2 + 0.024Q 4 
W) 1 + 1.29 Q 2 + 0.015 Q 4 



(X.28) 



In the neighbourhood of Q 2 = 0, the form factor is 
characterised by a radius- like length-scale; viz., 



6 W G^m 



(0.46 fm) 2 , (X.29) 



Q 2 =0 



which is almost indistinguishable from both rv = 0.45 fm 
in Table IX. II and the anomaly interaction radius defined 
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FIG. X.7. Solid curve - full result for scalar-diquark elastic 
electromagnetic form factor; and dashed curve - result ob- 
tained without dressing the quark-photon vertex. The com- 
puted mass of the diquark is m qq , = 0.776 GeV and the 
charge radius is given in Table IXlf 



in Ref. |7a |; viz., T o = 0.48 fin. On the other hand 

lim G^ P {Q 2 ) = 0.11, (X.30) 



FIG. X.8. Solid curve - Pseudovector-diquark electric form 
factor, G^j (Q 2 ), which exhibits a zero at Q 2 = 6.5 GeV 2 . 
(In this case 1 - §r? = for Q 2 = 6m 2 + = 6.7 GeV 2 , 
given the computed mass of 1.06 GeV.) The dashed curve, 

G\ i (Q 2 )/h 1 +, and dot-dashed curve, Gq (Q 2 )/Qi+, are al- 
most indistinguishable, as emphasised by the dotted curve, 

[G^(Q 2 )/mp]/[G^ + (Q 2 )/Q.i+]- The charge radii, and mag- 
netic and quadrupole moments are given in Table IXll NB. 
All form factors exhibit a pole at Q 2 
quark-photon vertex is dressed as described in Sec. lVII 



m p because the 



owing to the presence of the pion's pseudovector com- 
ponent, a result in keeping with the pointlike nature of 
bound-states generated by a contact-interaction [7^, [7(| . 

Scalar-diquark elastic form factor 
In the context of the interaction in Eq. (|VII.2j) . 
a detailed discussion of the relationship between 
pseudoscalar- and vector-mesons and scalar- and axial- 
vector-diquark correlations may be found in Ref. |l64j . 
Using the information provided therein, it is straight- 
forward to show that in rainbow-ladder truncation the 
electromagnetic form factor of a scalar diquark is readily 
obtained from the expression for F% m (Q 2 ). Namely, 



qq o+ 



(X.31) 

where the scalar-diquark Bcthc-Salpctcr amplitude is ex- 
pressed via (C = 7274 is the charge-conjugation matrix) 



r 



9?o+ 



(P)Ct 



75 



1 



iE qq AP) + — 1 .PF qq AP) 



J m + V + ) -r M I * * qq + < 

(X.32) 

The result for the scalar diquark elastic electromag- 
netic form factor is presented in Fig. IX.7l An interpola- 
tion valid on Q 2 G [—m% 10 GeV 2 ] is 



pom 



(0 2 ) int 4 rp 



1 1 + 0.25 Q 2 + 0.027 Q 4 
3 1 + 1.27Q 2 + 0.13Q 4 



(X.33) 



The normalisation is different but the momentum- 
dependence is similar to that of F% m . This is indicated, 
too, by the ratio of charge radii; viz., r +/?v = 1.08, 
whic h ma y be compared to the value of 1.09 obtai ned in 
Ref. 12381 and contrasted with the value of 0.8 in 2571. 



In the absence of the scalar-diquark Bethe-Salpeter am- 
plitude's pseudovector component, F qqQ+ = 0, we find 
r + = 0.51 fm; i.e., an increase of 6%. 

Pseudovector-diquark elastic form factors 
From the above observations it will be apparent that 
the rainbow-ladder results for the {ud} axial-vector di- 
quark elastic form factors may be obtained directly from 
those of the p-meson through the substitutions 



F? 



(X.34) 



The momentum-dependence of the form factors for the 
{uu} and {dd} correlations is identical but in these cases 
the normalisations are, respectively, | and — |. 

The axial-vector diquark form factors are depicted in 
Fig. IX. 81 They are similar to but distinguishable from 
those of the p-meson, falling-off a little less rapidly owing 
to the larger mass of the axial-vector diquark. Interpo- 
lations valid on Q 2 e [— m 2 , 10 GeV 2 ] are 



G l+ (g2) interp. 



G^(0 2 ) int ^ rp - 



G 1 +(Q2) interp. 



1 -0.16Q 2 



1 + 1.17 Q 2 + 0.012 Q 4 ' 
2.13-0.19Q 2 

i + 1.07 q 2 - o.iog 4 ' 

0.81 - 0.029 g 2 
~1 + 1.11Q 2 - 0.054 Q 4 



(X.35) 
(X.36) 
(X.37) 



from which the particular pseudovector diquark form fac- 
tors are obtained after multiplication by the appropriate 
charge factors, listed in Table DP] 
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FIG. X.9. Solid curve - momentum-dependence of full re- 
sult for axial-vector-scalar-diquark transition form factor, 
G 0+7l+ (Q 2 ); and dashed curve - result for G" r7P (Q 2 ) in 
Fig, IX. 61 The different rates of evolution are typical of me- 
son cf. diquark form factors. Note that e{ u[ j}<; + 7l +m 99 + = 
e { „ d} 0.74 = 0.25. 



l^-O"*" diquark transition form factor 
Owing to the flavour structure of the scalar diquark, 
this transition can only involve the {ud} axial- vector di- 
quark. It is described by a single form factor, which can 
be introduced through 



T°> 1+ (fci,fc 2 ) 



l go+ 7 i+ 

3 771, 



k 2 p G°^(Q 2 ) 



99i+ 



(X.38) 

and one may readily determine that in rainbow-ladder 
truncation 

g o+ 7 i+ (q2) 

q 7! "/p(q 2 ) \ ( Ei " Fi " E p^V^( E ii + ' F ii + ' E ii 1+ ) (X 39) 



The computed result for the transition form factor [73] 
is depicted in Fig. IX. 91 An interpolation valid on Q 2 e 



-m^lOGeV 2 ] is 



G 0+7l+ (Q 2 ) 



interp. 1 + 0.10Q 2 



(X.40) 
(X.41) 
(X.42) 



1 + 1.073 Q 2 ' 

The associated transition radius is 

^0+71+ = 0.48 fm, 
which is 5% larger than r^ lp in Eq. (|X.29|) . and 

lira G 0+7l+ (Q 2 ) = 0.049, 

Q 2 ->oo 

just under one- half of the value in Eq. (pC30|) . 

C. Nucleon form factors 



As we have seen in connection with the properties 
of mesons including at least one heavy quark, Sec. lIXI 




FIG. X.10. Proton Dirac (upper panel) and Pauli (lower 
panel) form factors, as a fu nction of x — Q 2 /m 2 N . Solid 
curve - result obtained |l43l ] using the contact-interaction, 
Eq. (|VII.2I) . and hence a dressed-quark mass- function and 
diquark Bethe-Salpeter amplitudes that are momentum- 
independent; dashed curve - result obtained in Ref. [33l ]. 
which employed QCD-like momentum-dependence for the 
dressed-quark propagators and diquark Bethe-Salpeter am- 
plitudes in solving the Faddeev equatio n; do t-dashed curve - 
a parametrisation of experimental data [258f ] . 



the computation of form factors can often be simpli- 
fied by employing algebraic parametrisations of one or 
more elements (propagators and/or amplitudes). Such 
parametrisations, based upon numerical solutions of the 
gap- and Bethe-Sal peter equation s, have long been em- 
ployed efficaciously jl38l |202| . l259f . Following that path, 
Ref. [lH produced a comprehensive survey of nucleon 
electromagnetic form factors, the quality of which is char- 
acterised by the results depicted in Fig. IX.10l With these 
results, Ref. [33[ unified the computation of meson and 
nucleon form factors, and also their valence-quark distri- 
bution functions (see Sec. lXDl) . 

Figurc JX.lOl also depicts nucleon form factors ob- 
tained [l43| with a symmetry-preserving DSE-treatment 
of the contact interaction in Eq. (|VII.2[) and hence the 
diquark form factors described in Sec. lXBI These form 
factors characterise a nucleon that is constructed from di- 
quarks whose Bethe-Salpeter amplitudes are momentum- 
independent and drcsscd-quarks with a momentum- 
independent mass-function, which inputs to the Faddeev 
equation yield a bound-state described by a momentum- 
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FIG. X.ll. Neutron Dirac (upper panel) and Pauli (lower 
panel) form factors, as a fu nction of x = Q 2 /m%. Solid 
curve - result obtained |l43l ] using the contact-interaction, 
Eq. (|VH,2|) . and hence a dressed-quark mass- function and 
diquark Bethe-Salpeter amplitudes that are momentum- 
independent; dashed curve - result obtained in Ref. [33| . 
which employed QCD-like momentum-dependence for the 
dressed-quark propagators and diquark Bethe-Salpeter am- 
plitudes in solving the Faddeev equatio n; do t-dashed curve - 
a parametrisation of experimental data |258l |. 



independent Faddeev amplitude. This last is the hall- 
mark of a pointlikc composite particle and explains the 
hardness of the computed form factors, which is evident 
in Figs. IX. 101 The hardness contrasts starkly with results 
obtained from a momentum-dependent Faddeev ampli- 
tude produced by dressed-quark propagators and diquark 
Bethe-Salpeter amplitudes with QCD-like momentum- 
dependence; and with experiment. As we have remarked 
above, evidence for a connection between the momentum- 
dependence of each of these elements and the behaviour 
of QC D's /3 -function is accumulating; e.g., Refs. [3ol [7a - 
El HM [Till- The comparisons in Figs. [XTlO] add to 
this evidence, in connection here with readily accessible 
observables, and support the view, presented herein, that 
experiment is a sensitive probe of the running of the /3- 
function to infrared momenta. 

This returns us to the ratio G P E {Q 2 ) / G P M {Q 2 ) de- 
scribed in the introduction to this section and depicted 
in Fig. IX.12l One observes that, independent of the na- 
ture of the interaction, computations of this ratio exhibit 
a zero. Its location, however, is very sensitive to the in- 
teraction, although not to the electromagnetic size of the 
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FIG. X.12. Upper panel: Normalised ratio of proton Pau li and 
Dirac form factors. Solid curve - co ntact interaction |143| ]; 
long- dashed curve - result from Ref. |l59l |. which employed 
QCD-like momentum-dependence for the dressed-quark prop- 
agators and diquark Bethe-Salpeter amplitudes; long-dash- 
dotted curve - draw n from parametrisation of experimen- 
tal data in Ref. 258]; and dotted curve - softened contact- 
interaction result, described in connection with Eq. (|X.43[) . 
Lower panel: Normalised ratio of proton Sachs electric 
and magnetic form factors. Solid curve and long-dashed 

curve, as above; dot- d ashed curve - linear fit to data in 

Refs. [234, 252, 260 2621 ] . constrained to one at Q 2 = 0; short- 
dashed curve - [1, 1]-Pade fit to that data; and dotted curve 
- softened contact-interaction result, described in connection 
with Eq. (|X.43[) . In addition, we ha ve re presente d a selec- 
tion of d ata e xplicitly: filled-squ ares [26(il |; circles [262f | : up- 
triangles [252l |: and open-squares [263l |. 



diquark correlations [33[ . (We note that the DSE predic- 
tion of the same ratio fo r the neutron has been confirmed 
by recent experiments [253[ , as also is the trend of the 
nucleoli's flavo ur separated Dirac and Pauli form factors 
[Til l25l . [2641 . ) 

In order to assist in explaining the origin and location 
of a zero in the Sachs form factor ratio, in the top panel 
of Fig. IX.12l we depict the ratio of Pauli and Dirac form 
factors: both the actual contact-interaction result and 
that obtained when the Pauli form factor is artificially 
"softened;" viz., 



F 2p (Q 2 ) 



F2 P (Q 2 ) 



1 + Q 2 /(4m 2 w ) 



(X.43) 



As observed in Ref. 265], a softening of the proton's Pauli 
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form factor has the effect of shifting the zero to larger 
values of Q 2 . In fact, if F 2p becomes soft quickly enough, 
then the zero disappears completely. 

The Pauli form factor is a gauge of the distribution 
of magnetisation within the proton. Ultimately, this 
magnetisation is carried by the dressed-quarks and in- 
fluenced by correlations amongst them, which are ex- 
pressed in the Faddccv wave-function. If the dressed- 
quarks are described by a momentum-independent mass- 
function, then they behave as Dirac particles with con- 
stant Dirac values for their magnetic moments and pro- 
duce a hard Pauli form factor. Alternatively, suppose 
that the dressed-quarks possess a momentum-dependent 
mass-function, which is large at infrared momenta but 
vanishes as their momentum increases. At small mo- 
menta they will then behave as constituent-like parti- 
cles with a large magnetic moment, but their mass and 
magnetic moment will drop toward zero as the probe mo- 
mentum grows. (N.B. As described in Sec. I VI Bl mass- 
less fermions do not possess a measurable magnetic mo- 
ment.) Such dressed-quarks will produce a proton Pauli 
form factor that is large for Q 2 ~ but drops rapidly 
on the domain of transition between nonperturbativc 
and perturbative QCD, to give a very small result at 
large-Q 2 . The precise form of the Q 2 -dcpendcnce will 
depend on the evolving nature of the angular momen- 
tum correlations between the dressed-quarks. From this 
perspective, existence, and location if so, of the zero in 
H p Ge p (Q 2 )/Gm p (Q 2 ) are a fairly direct measure of the 
location and width of the transition region between the 
nonperturbativc and perturbative domains of QCD as 
expressed in the momentum-dependence of the dressed- 
quark mass-function. 

One may expect that a mass-function which rapidly be- 
comes partonic - namely, is very soft - will not produce 
a zero; has seen that a constant mass-function produces 
a zero at a small value of Q 2 , and knows that a mass- 
function which resembles that obtained in the best avail- 
able DSE studies [HI, [57} and via lattice-QCD simulations 
[ll| , produces a zero at a location that is consistent with 
extant data. There is plainly an opportunity here for 
very constructive feedback between future experiments 
and theory. It is anticipated that experiments at JLab in 
the 12 GeV era will establish conclusively whether or not 
this ratio possesses a zero. The result will assist greatly 
in refining understanding of the dressed-quark mass func- 
tion and therefrom QCD's /3-function. 



D. Valence-quark distributions at x = 1 

Before closing this section we would like to exploit a 
connection between the Q 2 = values of elastic form 
factors and the Bjorken-x = 1 values of the dimensionless 
structure functions of deep inelastic scattering, F^' p (a;). 
First recall that the x = 1 value of a structure function is 
invariant under the evolution equations, Scc. lVIIll Hence 



the value of 



d v (x) 



u v (x) 



where 



d v {x) 
u v {x) 



aE£M 



1 



F«(x) 



(X.44) 



is a scale-invariant feature of QCD and a discriminator 
between models. Next, when Bjorken-x is unity, then 
Q 2 + 2P ■ Q = 0; i.e., one is dealing with elastic scat- 
tering. Therefore, in the neighbourhood of x = 1 the 
structure functions are determined by the target's elas- 
tic form factors. The ratio in Eq. (|X.44|) expresses the 
relative probability of finding a e?-quark carrying all the 
proton's light-front momentum compared with that of a 
w-quark doing the same or, equally, owing to invariance 
under evolution, the relative probability that a Q 2 = 
probe either scatters from a d-quark or a u-quark; viz., 



d v (x) 



UyiX) 



(X.45) 



In constituent-quark models with 5'f7(6)-symmctric 
spin-flavour wave-functions the right-hand-side of 
Eq. (|X.45[) is 1/2 because there is nothing to distinguish 
between the wave-functions of u- and d-quarks, and the 
proton is constituted from M-quarks and one <i-quark. 
On the other hand, when a Poincarc-covariant Faddccv 
equation is employed to describe the nuclcon, 
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(X.46) 



where we have used the notation of Ref. [33J . Namely, 
pp,s _ pi JQ 2 = 0) is the contribution to the proton's 
charge arising from diagrams with a scalar diquark com- 
ponent in both the initial and final state: u[ud](3"/(3u[ud\. 
The diquark-photon interaction is far softer than the 
quark-photon interaction and hence this diagram con- 
tributes solely to u v at x = 1. Pf a = F? p (Q 2 = 0), 
is the kindred axial- vector diquark contribution; viz., 
2d{uu} (g) 7 ® d{uu} + u{ud} (8 7 ® u{ud}. At x = 1 
this contributes twice as much to d v as it does to u v . 
p p,m _ p™(Q 2 = 0), is the contribution to the pro- 
ton's charge arising from diagrams with a different di- 
quark component in the initial and final state. The 
existence of this contribution relies on the exchange of 
a quark between the diquark correlations and hence it 
contributes twice as much to u v as it does to d v . If 
one uses the "static approximation" to the nucleon form 
factor, Eq. (|X. 1|) . as wi th the treatment of the contact- 
interaction in Ref. |l43| . then P^' m = Q. It is plain from 
Eq. (|X.46j) that d v /u v — in the absence of axial- vector 
diquark correlations; i.e., in scalar-diquark-only models 
of the nucleon, which were once common and, despite 
their weaknesses, still too often employed. 

Using the probabilities presented in Refs. [33L Il43j . re- 
spectively, one obtains: 





pP,s 

r i 


Pf° 


pp,m 


a 1 


F 2 " 
Fj 


M(p 2 ) 


0.60 


0.25 


0.15 


0.28 


0.49 


M=constant 


0.78 


0.22 





0.18 


0.41 



(X.47) 
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Both rows in Eq. (|X.47|) arc consistent with d v /u v = 
0.23 ± 0.09 (90% confidence level, F 2 "/F| = 0.45 ± 0.08) 
inferred recently via con sideration of electron-nucleus 
scattering at x > 1 |266l |. On the other hand, this is 
also true of the result obtained through a naive consid- 
eration of the isospin and helicity structure of a proton's 
light- front quark wave function at x ~ 1 , which leads one 
to expect that d-quarks are five-times less likely than u- 
quarks to possess the same heli city as the proton they 
comprise; viz., d v /u v = 0.2 [197| . Plainly, contemporary 
experiment-based analyses do not provide a particularly 
discriminating constrain t. Fu ture experiments with a tri- 
tium target should help |267| , emphasising again the crit- 
ical interplay between experiment and theory in elucidat- 
ing the nature of the strong interaction. 



XI. STRONGLY-COUPLED QUARKS 
IN-MEDIUM 

In this section we review a small selection of the recent 
progress made using the DSEs on QCD's phase transi- 
tions at nonzero temperature (T) and chemical potential 
(/i) , transitions which are the f ocus of much experimental 
and theoretical attention j268| . In this regard, owing to 
the absence of a probability measure at /i ^ in lattice- 
QCD; i.e., the fermion sign problem [269], the DSEs are 
a valuable tool because of their ability t o co nnect con- 
finement and DCSB in the continuum [21 |270| |. 



A. Critical endpoint and phase coexistence 

Concerning the phase diagram in terms of temperature 
(T) and chemical potential (/j,), besides the general phase 
separation line, there are some other interesting possibil- 
ities, such as: a domain of color superconductivity; and 
the existence of a critical endpoint (CEP) and an associ- 
ated coexistence region. In the chiral-limit theory, a CEP 
marks the end of a line of second-order chiral-symmctry- 
restoring (and possibly deconfining) transitions, origi- 
nating on the temperature axis in the (/x, T)-plane, and 
the beginning of a line of first-order transitions. Such 
a critical endpoint would have observable consequences 
27ll . l272 |. so it is imperative to demonstrate its existence, 



determine its location and demarcate the subsequent do- 
main of phase coexistence. 

Despite the fermion sign problem, attempts have been 
made with lattice-QCD, using mat hematical devices to 
extrapolate from /Lt = 0. They yield [III [Hi]: ^/T c = 
0.73-1.06 and T E /T c « 0.95, and a signal for a mate- 
rial phase coexistence region |274j . However, it is not 
yet certain whether the existence of a CEP survives in 
simulations with lattice pa rame ters that more closely re- 
semble the physical world |275j . 

A wide variety of models have also been used to search 
for a CEP . The Nambu-Jona-Lasinio (NJL) type yield 
[22l,[Iz3: Li E /T c « 1.7, T E /T C w 0.4; their Polyakov- 



loop extensions produce [278^2811 : fi E /T c = 1.5-1.8, 
T E /T C = 0.3-0.8; and a chiral quark model gives [HI 
(H E ,T E )/T C = (2.0,0.4). On the other hand, by adding 
six- and eight-fermion interactions to the NJL model, one 
can obtain (fi E ,T E )/T C = (1.1,0.8), which is similar to 
the result produ ced b y a Polyakov-loop-augmented chi- 
ral quark model [HI (fi E ,T E )/T c = (0.9,0.8). The for- 
mer, mutually consistent results for the CEP's location 
conflict markedly with those obtained from lattice-QCD: 
/i E /T c is significantly larger and T E /T C , much smaller. 
If they are nevertheless correct, then finding the CEP in 
experiment will be difficult because modern colliders are 
restricted to exploration of the small- /i domain. Given 
this observation, it is unsurprising that an analysis of 
flow data fro m th e relativistic heavy ion collider leads to 
the estimate [Hi]: fi E /T c > 1.0 and T E /T C < 1.0. 

Simple DSE truncations have also been applied to the 
CEP problem. A confining zero-width momentum-space 
interaction, the antithesis of the NJL-modcl, produces 
241 u E /Tr = 0, T E /T C = 1; and a separable-interaction 
285j : /U E /T C = 1.09, T E /T C = 0.78. However, neither 
study described a region of coexisting phases. Notwith- 
standing that, in this chain of remarks about the model 
results, there is a hint that the length-scale characterising 
confinement in the quark-antiquark interaction markedly 
influences the location of the CEP. Another notable fea- 
ture of the studies is that the size of the coexistence re- 
gions associated with the CEP is very model dependent. 

A novel principle for locating the CEP and demarcat- 
ing the coexistence region was introduced in Ref. [HI- 
The basic tools are the chiral susceptibility and the gap 
equation. As we have seen, the gap equation's kernel is 
defined by a contraction of the dressed-gluon-propagator 
and -quark-gluon-vertex. In order to exemplify the 
method, Ref. [H| used a propagator that can interpo- 
late between models of the non-confining NJL-type and 
the confining interactions used in efficacious DSE studies 
of hadron observables [1, 0, [H, [12] , whilst always provid- 
ing a super-renormalisable interaction. For the vertex, 
the studied used either the rainbow-truncation; i.e., the 
leading-order term in a symmetry-preserving scheme [62| , 
or a dressed-vertex Ansatz. The capacity to draw the 
phase diagram derived from an arbitrary dressed-vertex 
is an essentially new feature of the method. 

At T 7^ /i, the gap equation is (uj n = uj n + ifi) 

S(p,w„) _1 = ij -p + ij^n +m + £(p,u) n ) , (XI. 1) 

1— — QO 

X y7^(<f,Wz)yr^(<f,W ; ,p,W n ) , (XI.2) 

where: tu n = (2n + 1)ttT is the fermion Matsubara fre- 
quency; fl n i = ui n — u>i; D^ v is the dressed-gluon prop- 
agator; and r„ is the dressed-quark-gluon vertex. (As 
an ultraviolct-finitc model is employed, rcnormalisation 
is unnecessary and m = in Eq. (jXI.ll) defines the chiral 
limit.) 
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In-medium, the gap equation's solution can be ex- 
pressed as 

+ i l4 Q n C(p 2 , S>1) + B{p\ul) , (XI.3) 

with, e.g., B(p 2 ,uj 2 )* = B(p 2 ,ui 2 _ n _ 1 ). The dressed- 
gluon propagator has the form 

g 2 D fW (k, n nl ) = PlD T (k 2 ,n 2 u ) + p^D L (k 2 ,n 2 u ) , 

. (XL4) 

where Pfy L are, respectively, transverse and longitudinal 
projection operators, Eq. (|XI.23|) . Whilst for T^O^/i 
it is generally true that Dt 7^ Dl, there are indications 
j286l | that for T < 0.2 GeV, the domain relevant here, 
it is a good approximation to treat Dt = Dl =: Dq. 
For the in-vacuum interaction, a simplified form of the 
interaction in Ref. [99[ was used; viz., with k = k 2 



O 2 



4tt 2 1 2 

D (k) = D— kc- k ^ 



(XI.5) 



The parameters in Eq. (|XI.5[) are D and a but it has long 
been known that they are not independent: a chan ge in 
D can be compensated by an alteration of a [215| |. For 
(7 e [0.3, 0.5] GeV, using Eq. (|XL6p below, ground-state 
pseudoscalar and vector-meson observables are roughly 
constant if oD = (0.8 GeV) 3 . A value of a = 0.5 GeV was 
usually employed in Ref. [28|. N.B. Eq. (|XI.5p was used 
for illustrative simplicity, not out of necessity. The most 
up-to-date interaction is described in Ref. [57] and the 
status of propagator and vertex studies can be tracked 
from Ref. [8|. 

As explained above, the gap equation is complete once 
the vertex is specified. Here results obtained using the 
rainbow-truncation: 



T v (q,u>i,p,uj n ) = 7„ . 



(XI.6) 



the first term in a symmetry- preserving scheme [621 ] , are 
compar ed with those produced by the Ball-Chiu (BC) 
Ansatz fl36l [287l . [288| : 

iT u (q,u} h p,ui n ) = i 7 j£ A + i^E c 



+ (Pn + gi)f. 



:lI(Pn + qi)a^A 



+ ^{Pn+qi) a ^C + ^B 



(XI.7) 



p n = (p, ui n + ifi), qt = (q, u>i + ifi), with (F = A, B, C) 



Z F (q 2 ,oj?,p 2 ,LJ 2 n ) = -[F(q 2 ,oj?) + F(p 2 ,uj 2 n )}, (XI.8) 

. ^ , F (q 2 ,uj?) - F(p 2 ,uj 2 ) 
A F (q 2 ,u H ,p 2 ,cj n )= yq % ^ " ; , XI.9 

qf - pI 

where, defining u = (0, 0, 0, 1), 7 J = j u - u u j a u a , 7^ = 
u u "f a u a . The comparison is natural because vertices of 
the type in Eq. (|XI.6[) have widely been used in studies 
of hadron observables 

n n n m, and the bc Ansatz 
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FIG. XI. 1. Panel (a) - Temperature dependence of the chiral- 
symmetry order parameter in Eq. (|XI,11[) . Chiral susceptibil- 
ity computed in the Wigner phase - Panel (b), and in the 
Nambu phase - Panel (c). In all panels: the Ball-Chiu vertex 
was used, Eq. (|XI.7|) ; and p = 0, D = 0.5 GeV 2 , m = 0. 



provides a semi-quantitatively accurate representation of 
latticc-QCD results for important terms in T n at T = 
= n (see Sec. fVTB]) . (It is plain from Eq. (|XT2|) that 
with the BC vertex the effective interaction strength is 
D = DA(0, 0). In contrast, owing to Eq. (|XL6|) . D = D 
in rainbow-ladder truncation.) 

The gap equation formulated above is readily solved 
and in Fig. IXI.il w e depict the T-dependence of a chiral 
susceptibility [289j and a chiral-symmetry order param- 
eter [Hi 



X (0,wo) = ^B(6,lo 2 ). 



(XI.10) 



-(qq)° = N C T tr D J J^L S m=0 (p>„) . (XI.ll) 



n— — 00 



For T < T E the behaviour of the order parameter is typ- 
ical of models wit hout long-range correlations in the gap 
equation's kernel [2891 ] . Namely, initially slow evolution 
from its T = value: (qq)° = (-0.258 GeV) 3 , which sig- 
nals chiral symmetry realised in the Nambu mode; i.e., 
dynamically broken chiral symmetry, and this followed 
by a mean-field transition to a phase with chiral symme- 
try restored; i.e., realised in the Wigner mode. In the 
Wigner phase the dressed-quarks are not confined. 

The lower panels in Fig. IXI.il show the chiral suscep- 
tibility of the Wigner and Nambu phases, which corre- 
spond to gap equation solutions that are, respectively, 
within t he do main of attraction of the B = or B ^ 
solution [290l ]. A phase is unstable in response to fluc- 
tuations if the susceptibility is negative, but stable and 
realisable otherwise. With \i = 0, one sees the Nambu 
phase completely replaced by the Wigner phase at T = 
124 McV. 

This should be contrasted with the behaviour in 
Fig. IXI.21 At T — the order parameter remains con- 
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FIG. XI. 2. Analogue of Fig. IXI.il displaying evolution with 
chemical potential at T = 0. 



stant with increasing fi until fi a = 0.30 GeV, which is 
the upper bound on the domain of analyticity for this 
gap equation's kernel [26|. On a small domain beyond 
this; viz., fi £ (fi a ,fi^), with /j£ = 0.314 GeV, the or- 
der parameter diminishes smoothly, an effect that may 
be denominated a partial restoration of chiral symme- 
try. For fi > fi^ the order parameter vanishes so that 
chiral symmetry is completely restored via a first-order 
transition. 

The lower panels of Fig. IXI.2l provide extra informa- 
tion. For n < (j,™ = 0.286 GeV, the Wigner phase is 
unstable. That changes at fi^ , when \w switches sign 
and thereafter, on the domain fi^f < fi < fj,^ , both 
the Wigner- and Nambu-phasc susceptibilities are pos- 
itive. This is the domain of phase coexistence, with 
a metastable Wigner phase for fi^ < fj, < fx a and a 
metastable Nambu phase for fi a < \i < fi^ . The pres- 
sure of the phases is equal at fi = \i a \ and the Nambu 
phase is completely displaced by the Wigner phase for 



fi > fi^ . Notably, with an Ansatz for the dresscd-quark 
gluon vertex, the diagrammatic content of the gap equa- 
tion's kernel is generally unk nown . However, owing to 
the insights provided in Ref. j29l| . one can draw these 
conclusions despite being unable to calculate an explicit 
expression for the thermodynamic pressure. 

At the onset of the coexistence domain one can expect 
pockets of deconfined, chirally symmetric quark matter 
to appear in the confining Nambu medium. Their num- 
ber and average volume will increase with fi. The oppo- 
site situation occurs at the termination of the domain; 
i.e., it is the Nambu phase which exists only in pockets. 
For fi £ (fi a , fj,^ ) , which is the domain of Nambu-phase 
metastability, the properties of observed hadrons will be 
affected by the partial restoration of chiral symmetry. 
The coexistence domain is not of a quarkyonic nature 
|270j . 

The computed phase diagrams [28| are drawn in 
Fig. IXI.31 The upper panel was obtained with the rain- 
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FIG. XI. 3. Chiral-limit phase diagram in the tempera- 
ture/chemical-potential plane for strongly-interacting quarks. 
The critical endpoint (CEP) is marked explicitly. Upper 
panel - rainbow vertex, Eg. (|XI.6[) . with D = 1.0 GeV 2 , 
a — 0.5 GeV. The /i = critical temperature is T° = 
0.133 GeV. Lower panel - Ball-Chiu vertex, Eq. ()XI.7[) . with 
D = 0.5 GeV 2 , a = 0.5 GeV. The p = critical temperature 
is T° = 0.124 GeV. 



bow truncation, Eq. (|XI.6[) . in which case the diagram- 
matic content of the gap equation's kernel is known and 
one can thus compute the dressed-quark component of 
the pressure. It is given by the auxiliary-field effective 
action evaluated at its extrema [l], Q . Within the domain 
(fi < fi E ,T > T E ) the system exhibits a mean- field tran- 
sition, which is signalled both by: equality of the Wigncr- 
and Nambu-phasc pressures; and coincident singularities 
in the Wigner- and Nambu-phase chiral susceptibilities. 
The curve tracking the singularity location in the Wigner 
and Nambu susceptibilities bifurcates at the CEP, with 
a domain of phase coexistence opening. Naturally, the 
curve of equal thermodynamic pressure lies within this 
domain. 

The lower panel in Fig. IXI.3l was obtained using the 
dressed vertex, Eq. (|XI.7[) . Its features are similar to 
those displayed in the upper panel but the domain of 
phase coexistence is smaller. Here, one cannot derive 
a form for the dressed-quark pressure. Thus it is only 
an appreciation of the information contained in the chi- 
ral susceptibilities that enables a phase diagram to be 
drawn. 

As explained in Sec lIIIl confinement is expressed in 
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TABLE XI. 1. Parameter- and vertex-dependence of the crit- 
ical endpoint and coexistence region. Ac is the width of the 
coexistence region on the /i-axis. Each parameter set pro- 
duces similar results for the in- vacuum values of the so-called 
vacuum quark condensate (qq) and the pion's leptonic decay 
constant f n Lattice-QCD suggests (fi E ,T E )/T c = (0.73- 
1.06,0.95). (Dimensioned quantities in GeV.) 



model 


result 


vertex 


£1/2 


a 


Tc 


Ac 


(^E,r E )/r c 




BC 


0.7 


0.50 


0.124 


0.026 


(1.13,0.89) 


1.27 


BC 




0.45 


0.128 


0.048 


(0.69,0.92) 


0.75 


BC 




0.40 


0.139 


0.076 


(0.16,0.96) 


0.17 


Bare 


1.0 


0.50 


0.133 


0.220 


(0.98, 0.90) 


1.08 


Bare 




0.45 


0.136 


0.280 


(0.81,0.89) 


0.91 


Bare 




0.40 


0.148 


0.360 


(0.17,0.95) 


0.18 



the analytic structure of the dressed-quark propagator. 
Measured in this way, it is significant that the models 
considered in this subsection are members of a class in 
which chiral symmetry restoration is accompanied by a 
coincident dressed-quark deconfinement transition in the 
chiral limit. 

Tabic IXLll illustrates the response of the CEP's loca- 
tion to changing the vertex or the parameters. Defin- 
ing a confinement length-scale r CT = l/cr, it is appar- 
ent that the CEP rotates toward the temperature axis 
as r a is increased. The extreme case is r„ = oo, 
which was computed in Ref. [24| and reported above: 
(/i E ,T E )/T c = (0, 1). Models of the NJL-type, as they 
have usually been used in the current context, represent 
the opposite limiting case: they are expressed via a gap 
equation in which the confinement length-scale vanishes. 
From this perspective, it is unsurprising that they pro- 
duce a CEP whose angular separation from the /x-axis is 
significantly smaller. 

This subsection reviews a method, based on the chi- 
ral susceptibility, which enables one to draw a phase di- 
agram in the chemical-potential/tcmperature plane for 
quarks whose interactions are described by any sensibly- 
constructed gap equation. Thus, in attempting to chart 
the phase structure of QCD using the methods of con- 
tinuum quantum field theory, one is no longer restricted 
to the simplest class of mean-field kernels: sophisticated 
quark-gluon vertices can be used. The method is gen- 
eral and potentially useful in all branches of physics that 
explore the properties of dense fermionic systems. 

A class of models that successfully describes in- vacuum 
properties of tt- and p-mesons, exhibits a critical endpoint 
(CEP) in the neighbourhood (m E ,T e ) - (1.0,0.9)T C . 
The CEP's angular separation from the temperature axis 
is a measure of the confinement length-scale: the separa- 
tion decreases as the confinement length-scale increases. 
Furthermore, a domain of phase coexistence opens at the 
CEP. It's size depends on the structure of the gap equa- 
tion's kernel but, other aspects being equal, it increases 



in area as the confinement length-scale increases. Wc 
are hopeful that illumination of the CEP and its conse- 
quences is within the reach of modern colliders. 



B. Quark spectral density and a sQGP 

It is believed that a primordial state of matter has been 
recre ated by the relativistic heavy-ion collider (RHIC) 
[292l | . This substance appears to behave as a nearly- 
perfect fluid on some domain of tem perature, T, above 
that required for its creation, T c [293]. An ideal fluid has 
zero shear- viscosity: 77 = 0, and hence no resistance to 
the appearance and growth of transverse velocity gradi- 
ents. A perfect fluid with near-zero viscosity is the best 
achievable approximation to that ideal . Gr aphene might 
provide a room temperature example [294j . From New- 
ton's law for viscous fluid flow; viz., dv/dz = (l/r))(F/A), 
it is apparent that in near-perfect fluids a macroscopic ve- 
locity gradient is achieved from a microscopically small 
pressure. Strong interactions between particles consti- 
tuting the fluid are necessary to achieve this outcome. 
Hence the primordial state of matter is described as a 
strongly-coupled quark gluon plasma (sQGP). 

We have explained that at T = 0, QCD is characterised 
by confinement and dynamical chiral symmetry breaking 
(DCSB). These phenomena are represented by a range of 
order parameters which all vanish in the sQGP. Under- 
standing the sQGP therefore requires elucidation of the 
behaviour and properties of quarks and gluons within this 
state. Perturbative techniques have been developed for 
use far above T r \ viz., the hard thermal loop (HTL) ex- 
pansion |295l . |296| , which has enabled the computation of 
gluon and quark thermal masses mj ~ gT and damping 
rates 7r ~ 9 2 T, with g = g(T) being the strong running 
coupling. It also suggests the existence of a collective 
plasmino or "abnormal" branch to the dressed-quark dis- 
persion relation, which is charac teris ed by antiparticle- 
like evolution at small momenta 
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Owing to asymptotic freedom, the running coupling in 
QCD increases asT-> T+. Therefore, a simple interpre- 
tation of the HTL results suggests the plasmino should 
disappear before T c is reached because 7t increases more 
rapidly than and "fx/'m-T ~ 1 invalidates a quasi- 
particle picture. On the other hand, lattice-regularised 
quenched-QCD suggests that the plasmino branch per- 
sists in the vicinity of T c |298| | . It is necessary to resolve 
the active degrees of freedom in the neighbourhood of 
T c because the spectral properties of the dressed-quark 
propagator are intimately linked with light-quark con- 
finement (Sec. lIIip and it is the long-range modes which 
might produce strong correlations. 

We saw in the preceding section that below T c dressed- 
gluons and -quarks are confined and chiral symmetry is 
dynamically broken; and beyond T c , deconfinement and 
chiral symmetry restoration occur via coincident second- 
order phase transitions. The sQGP seems to occur in the 
neighbourhood of T c , so in this section we review how the 
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gap equation may be used to elucidate the active fermion 
quasiparticles for T >T C [2^] . 

On the domain T > T c the chiral-limit dressed-quark 
propagator can be written 

S(iu n ,p) = -17 ■ pa A (u> n ,p 2 ) - ij4Uj n a c (uJn,p 2 ) ■ 

(XI.12) 

There is no Dirac-scalar part because chiral symmetry is 
realised in the Wigner mode but this does not mean that 
nonperturbative phenomena are excluded, as is appar- 
ent, for example, in the discussion of no vel Wigne r-mode 
solutions to the gap equation in Refs. [2901 |299( . The 
retarded real-time propagator is found by analytic con- 
tinuation 

S r (lo,p) = S(iu n ,?)\ iUn ^ +irl+ (XI.13) 
and from this one obtains the spectral density 

p(u,p) = -2ZS R (u J ,p). (XI.14) 



Equations (|XI.13|) and (|XI.14I) arc equivalent to the state- 
ment: 



b{iu} n ,p) = — / dui — 
2tt J_ 00 u' 



(XI.15) 



NB. If one requires a nonnegative spectral density, then 
Eq. (|XI.15|) is only valid for T > T c ; i.e., on the decon- 
flned domain (Sec. lIII|) . Furthermore, including a light 
current-quark mass has no material impact on the pri- 
mary conclusions of Ref. [2!| . 

For an unconfined dressed-quark propagator of the 
form in Eq. (|XI.12[) , the spectral density can be expressed 

p(u,p) = P+ (uj,p 2 )P+ + p_(u,p 2 )P_ , (XI.16) 



where P± = (74 ± ij ■ u p )/2, 



P 



\p\, are operators 



which project onto spinors with a positive or negative 
value for the ratio H := helicity/chirality: H = 1 for a 
free positive-energy fermion. The spectral density is in- 
teresting and expressive because it reveals the manner by 
which interactions distribute the single-particle spectral 
strength over momentum modes; and the behaviour at 
T ^ shows how that is altered by a heat bath. As with 
many useful quantities, however, it is nontrivial to eval- 
uate p(u), \p\). Nonetheless, if one has at hand a precise 
numerical determination of the dressed-quark propaga- 
tor in Eq. (|XI.3|) . then it is possible to obtain an accurate 
approximation to the spe ctra l den sity via the maximum 
entropy method (MEM) |300ll30l} . 

In the MEM the entropy functional of the spectral 
function is defined as 



/+00 
dui 
-00 



p(uj) — m(oj) — p{uj) lo 



p(uj) 



'm(oj) 



(XI.17) 

where m(uj) is the "default model" of the spectral func- 
tion. One typically adopts a simple form; viz., 



m(w) = m 9(A 2 - u 2 ) , 



(XI.18) 
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FIG. XI. 4. Dependence of the spectral density on the cutoff, 
A, in the default MEM model, Eq. (|XI.18|) . at T = 1.1 T c (up- 
per panel) and T — 3.0 T c (lower panel), where T c is discussed 
in Eq. l|XI,26p . 



which is a uniform distribution with no a priori struc- 
ture assumption. It should be borne in mind that one can 
only claim a reliable result from the MEM if the spectral 
function produced is insensitive to mo and A. The sensi- 
tivity to mo is uniformly weak in Ref. [29j ] . However, on 
a material domain, the sensitivity to A is strong. Nev- 
ertheless, for any given temperature, it was found that 
there is always a value A above which a stable spectral 
function is obtained. This is illustrated in Fig. IXI.4l 

In this instance one is interested in the chiral-limit 
dressed-quark propagator, which is obtained from (F = 
A,C) 



S(iu n ,p) 1 = Z 2 ij ■ p + Z 2 i'f4^n + T,'(iui n ,p) . 



(XI.19) 



£' (iuj n , p) = «7 ■ P'E'a ( iu} n , P) + «74 u n T,' c (iuj n , p) , 



(XI.20) 



T,' F (iuj n ,p) = 



d 3 q 
(2tt) s 



g 2 D f _ lu (uj n -uJi,p-q) 



x -tr D £P F 7 / _ 1 S , (ia;i,(7)r ; ,(w„,a;;,p,g), 



(XI.21) 



where: T A = -Z x tf-p/p , T c = -Zii^/uj n ; and Z 1>2 , 
Zf 2 are, respectively, the vertex and quark wave function 
renormalisation constants. The re gular isation and renor- 
malisation procedures of Refs. 0, 287 1 are followed. 
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Reference [29| used the rainbow-ladder trunca- 
tion and the phenomenologically-efRcacious one -loop 
renormalisation-group-improved interaction of Ref. |287| . 
Namely: 

g 2 D llv {uj n - oj h p- q)T v {u) n ,u) U p,q) 
= [P^(k n )D T (k n )+P^(k n )D L (kn)}lu, (XI.22) 

where k n := (fi, k) = (oj„ - 0Ji,p — q); 



5.0 A 



fi and/or v = 4 , 



fj,,v = 1,2,3. 



(XI.23) 



with P£" + P£" = 5^ - and 



D r (fcn) = ©(4,0), D £ (fcn) = ©(Ag.mJ) , (XI.24) 

+ , f -7^TT2 (XI.25) 
ln[r+(l + sj2/A^ CD ) 2 ] 

with f (sq) = (l-exp(-sn/4m|)/sQ, s n = ft 2 +fc 2 +m 2 , 

r = e 2 - 1, m t = 0.5 GeV, lrn = 12/25, and A^ C = A = 
0.234 GeV. For pseudoscalar and vector mesons with 
masses < 1 GeV, this interaction provides a uniformly 
good description of their T = properties Q when 
u = 0.4 GeV, D = (0.96 GeV) 2 . In generalising to T ^ 0, 
Ref. [2!| have followed perturbation theory and included 
a Debye mass in the longitudinal part of the gluon prop- 
agator: to 2 = (16/5)T 2 . This kernel produces coincident 
second-order deconfmement and chiral symmetry restor- 
ing transitions for two massless flavors at 

T c = 0.14 GeV , (XI.26) 

which is 10% smaller than that obtained in Ref. |302| . 
(As mentioned above, this interaction has only recently 
been superseded (57j.) 

One insufficiency of the interaction defined above is 
that -D, the parameter expressing its infrared strength, is 
assumed to be T-independent. Since the nonperturbative 
part of the interaction should be screened for T >T C , one 
may remedy that by writing D — > D(T) with 



D(T) = 



D . 



T <T, 



p ■ 



b + HT/A Q cd] 



T > TU 



, (XI.27) 



where T p is a 



persistence" temperature; i.e., a scale be- 
low which nonperturbative effects associated with con- 
finement and dynamical chiral symmetry breaking are 
not materially influenced by thermal screening. Loga- 
rithmic screening is typical of QCD and with a = 0.028, 
b = 0.56 the modified interaction yields tot = 0.8 T for 
T > 2T C ; vi z., a thermal quark mass consistent with 
lattice-QCD J298(. Unless stated otherwise, T p = T c 
herein. 

The spectral density can be obtained by employing the 
MEM in connection with the solution of our gap equa- 
tion. Notably, the behaviour changes qualitatively at 
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FIG. XI. 5. Left panel - Temperature-dependence of the 
dressed-quark thermal masses. Notably, spectral strength is 
associated with a massless-mode. Right panel - T-dependence 
of the residue associated with that zero mode. 



T c . Indeed, employing a straightforward generalisation 
of the inflexion point criterion introduced in Refs. [8j, [20J, 
one can readily determine that reflection positivity is vio- 
lated for T <T C . This signals confinement. On the other 
hand, the spectral function is nonnegative for T > T c . 

Figure IXLol depicts the T > T c -dependence of the loca- 
tions of the poles in p{uj,p = 0); i.e., the thermal masses. 
As conventionally anticipated, spectral strength is lo- 
cated at uj+(p = 0) and 0J-(p = 0) = — oj+{p = 0), corre- 
sponding to the fermion's normal and plasmino modes 
at nonzero temperature. However, it is striking that 
on a measurable T-domain, spectral strength is also 
associated with a quasiparticlc excitation described by 
uo(p = 0) = 0. The appearance of this zero mode is 
an essentially nonperturbative effect. It is an outgrowth 
of the evolution in-mcdium of the gap equation's T — 
Wigner-mode solution and analogous to this solution's 
pers istence at nonzero current-quark mass in vacuum 
[290j . This feature is stable for sufficiently- large A in 
Eq. (|XI.18|) . as evident in Fig.|XL! 

The spectral density possesses support associated with 
this zero mode on T e [0, T s ]. In fact: all the Wigncr- 
phase spectral strength is located within this mode at 
T = 0; it is the dominant contribution for T > T c ; and, 
while it is dominant, it is the system's longest wavelength 
collective mode. On the other hand, as evident in the 
right panel of Fig. IXI.51 the mode's spectral strength di- 
minishes uniformly with increasing T and finally vanishes 
at T s fn 1.35 T c . Then, for T > T s the quark's normal and 
plasmino modes exhibit behavior that is broadly consis- 
tent with HTL calculations. This is apparent in Fig. IXI.51 
and in a comparison between the upper and lower panels 
of Fig. IXI.61 Given these observations, one may reason- 
ably argue that the system should be considered a sQGP 
for T £ [T C ,T S ], whereupon it contains a long-range col- 
lective mode. 

It should be noted that the HTL approach is pertur- 
bative and only applicable for T » T c . Hence it could 
not have predicted the zero mode's existence. Numeri- 
cal simulations of lattice-QCD, on the other hand, are 
nonperturbative. However, it is practically impossible 
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FIG. XI.6. Upper panel, left - Quasiparticle dispersion re- 
lations, uj±fi(p), at temperature T = 1.1T C . The diagonal 
dotted line is the free-fermion dispersion relation at this T. 
Upper panel, right - The residues associated with these quasi- 
particle poles. Lower panel - Same information for T = 3T C , 
whereat the zero mode has vanished. 



in contemporary computations to exactly preserve chiral 
symmetry. This can plausibly explain the absence of the 
zero mode in lattice simulations because any source of 
explici t chi ral symmetry breaking heavily suppresses the 
mode [29fjj . 

The upper- left panel of Fig. lXLol depicts the dispersion 
relations for all dressed-quark modes that exist for T < 
T s and the behavior of their associated residues. On this 
sQGP domain the dispersion relations are atypical, with 



.,p~o -0.2 $, 
uj±(\p\) = m T ' 
+0.3 1$ , 

uj (\p\) p =°0.80 \p\. 



(XI.28a) 
(XI.28b) 



Notwithstanding this, all realise free-particle behaviour 
for \p\ 3> T. Moreover, the usual spectral sum rules are 
satisfied. Indeed, the identity 



zi 



duj' 



(XI.29) 



assists in understanding the momentum-dependence in 
the upper-left panel of Fig. IXI.6l The upper-right panel 
displays the momentum-dependence of the pole residues: 
spectral support is located completely in the normal 
mode for \p\ 3> T\ i.e., on the perturbative domain. 

The lower panel of Fig. IXI.6l characterises the be- 
haviour of p(w, \p\) on t he high-T domain. In agreement 
with HTL analysis 296] , expected to be valid thereupon, 
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FIG. XI. 7. Numerical check of the momentum sum rule in 
Eq. (|XI,29[) : the curves would be indistinguishable with a per- 
fect reconstruction of the spectral density. 

one finds only normal and plasmino modes, with 



uj±(\p\) P = m T ±0.33|$ 



(XI.30) 



The plasmino dispersion law exhibits the expected min- 
imum, in this case at |p|/T ~ 1; and both w±(|p|) ap- 
proach free-particle behaviour at \p\ ^> T, with that of 
the plasmino approaching this limit from below. The 
lower-right panel shows that the contribution to the spec- 
tral density from the plasmino is strongly damped and 
contributes little for p > 2T. These results are i n-lin e 
with those obtained via simulations of lattice-QCD j298| . 

It is anticipated from HTL analyses that T ^ prop- 
agators exhibit branch cuts whose appearance can be 
attributed to the op ening of scattering channels that 
are absent at T — [303j |. However, such branch cuts 
do not materially contribute to the nonperturbatively- 
determined spectral density described herein. This is 
plausible because a branch point is a lower-order nonan- 
alyticity than a pole; i.e., in numerical studies, poles are 
features with large height, small width and significant 
spectral strength, whilst branch points are low, broad 
features with lesser spectral strength. Thus, compared 
with poles, branch points can be invisible to a numerical 
procedure. Uncovering them requires fine tuning within 
the MEM, or any other method. 

One should nevertheless check whether the numerical 
result omits sig nificant spectral strength, and this was 
done in Ref. [231 i e -§-> ^ was verified that, to a high degree 
of accuracy, the spectral density satisfies the textbook 
relation: 



-t-oo 



du' 
27 



p±(u/,p) « Z 2 



E 

=0,+,- 



Z Q » 1 (XI.31) 



where Zq are the quasiparticle residues. Other, more de- 
manding sum rules, such as that in Eq. (|XI.29|) . were also 
checked with the result depicted in Fig. IXI.7l These tests 
indicate that no significant contribution to the dressed- 
quark's spectral strength was overlooked in applying the 
MEM. 
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Equation (|XI.27|) is a model and it is natural to enquire 
after its influence. None of the results in this subsection 
are qualitatively altered by varying T p but, as one would 
expect, the width of the sQGP domain expands slowly 
with increasing T p ; e.g., a 50% increase in T p produces a 
30% increase in T s . 

The appearance of a third and long-wavelength mode 
in the dressed-fermion spectral density on a material 
temperature domain above T c has also been observed in 
one-loop computations of the fermion self-energy, irre- 
spective of t h e na ture of the boson which dresses the 
fermion }304 |305| . This mode appears for T > tog, 
where tog is an infrared mass-scale associated with the 
boso n, and persists for fermion bare masses < 0.2tog 
|306l |. (NB. In the present case, tog « 0.17 GeV, hence 
0.2 tog ~ 34MeV, appreciably greater than the light- 
quark current-masses.) Where a comparison is possible, 
the dependence of the spectral density described in this 
subsection on (u,\p\,T) is similar to that seen in the 
one-loop analyses of model gap equations. In analogy 
with a sim ilar effect in high-temperature superconduc- 
tivity j307| . that behaviour has been attributed to Lan- 
dau damping, an interference phenomenon known from 
plasma physics. Indeed, Landau damping is typical of 
in-medium self-energy corrections when the thermal en- 
ergy of the fermion is commensurate with the mass-scale 
which characterises the dispersion law of the dressing bo- 
son. 

A feature that unifies Ref. H3| and Rcfs. HH [306j 
is the strictl y nonper turbative phenomenon of DCSB. 
Whilst Rcfs. |304l . l306j extract spectral densities via one- 
loop estimates in Nambu-Jona-Lasinio- or Yukawa-like 
models, the couplings are tuned to mimic a world in 
which chiral symmetry is dynamically broken. In this 
connection we re-emphasise that the zero mode appears 
only when the strength of the interaction is capable of 
producing DCSB in-vacuum; viz., when the gap equa- 
tion's T = kernel has sufficient support at infrared mo- 
menta. Hence, it is an essentially nonperturbative phe- 
nomenon. Fur thermore , given the qualitative differences 
between Refs. |304 l306l | and Ref. [29j . the similarity be- 
tween results suggests the possibility that the appearance 
of zero modes is model- independent. 

It is notable that a coupling to meson-like correlations 
in the gap equation is not a precondition for the appear- 
ance of the zero mode because such correlations are ab- 

I It 

sent in the rainbow truncation [289j . On the other hand, 
the kernel of the gap equation used in this subsection is 
characterised by an interaction that features an infrared 
mass-scale tog ^ T c and supports dynamical chiral sym- 
metry breaking at T = 0. We anticipate that the zero 
mode will markedly affect colour-singlet vacuum polari- 
sations on T e [T c , Tg\. This could be explicated using 
the methods of Refs. |63| . 

Experiment has presented us with the fascinating pos- 
sibility that a near-perfect fluid might have been a key 
platform in the universe's evolution. This subsection re- 
views the possibility that a zero mode exists just above 



the critical temperature associated with the strong phase 
transition, T c . The mode contains the bulk of the spec- 
tral strength for T > T c and so long this mode persists, 
the system may reasonably be described as a strongly- 
interacting state of matter. If the existence of this long- 
wavelength mode is model-independent, then it is natural 
to anticipate that a strongly-interacting state of matter 
should precede the QCD phase transition. 



XII. EPILOGUE 

QCD is the most interesting part of the Standard 
Model and Nature's only example of an essentially non- 
perturbative fundamental theory. Whilst confinement re- 
mains a puzzle, it is recognised that dynamical chiral 
symmetry breaking (DCSB) is a fact in QCD. It is man- 
ifest in dressed-propagators and vertices, and, amongst 
other things, it is responsible for: the transformation of 
the light current-quarks in QCD's Lagrangian into heavy 
constituent-like quarks, in terms of which order was first 
brought to the hadron spectrum; the unnaturally small 
values of the masses of light-quark pseudoscalar mesons 
and the rj-rj' splitting; the unnaturally strong coupling 
of pseudoscalar mesons to light-quarks - g-^qq ~ 4.3; and 
the unnaturally strong coupling of pseudoscalar mesons 
to the lightest baryons - g^NN ~ 12.8 « 3g W q q . 

Herein we have illustrated the dramatic impact that 
DCSB has upon observables in hadron physics. A 
"smoking gun" for DCSB is the behaviour of the 
dressed-quark mass function. The momentum depen- 
dence manifest in Fig. III.ll is an essentially quantum 
field theoretical effect. Exposing and elucidating its 
consequences therefore requires a nonperturbative and 
symmetry-preserving approach, where the latter means 
preserving Poincare covariance, chiral and electromag- 
netic current-conservation, etc. The Dyson-Schwinger 
equations (DSEs) provide such a framework. We have 
explained the nature of some of the experimental and 
theoretical studies that are underway which can poten- 
tially identify observable signals of M(p 2 ) and thereby 
confirm and explain the mechanism responsible for the 
vast bulk of visible mass in the Universe. 

Along the way we have described a number of exact 
results proved in QCD using the DSEs, amongst them: 

• Light-quark confinement is a dynamical phe- 
nomenon, which cannot in principle be expressed 
via a potential; 

• Goldstone's theorem is fundamentally an expres- 
sion of equivalence between the one-body prob- 
lem and the two-body problem in the pseudoscalar 
channel; 

• quarks arc not Dirac particles - they possess 
anomalous chromo- and electro-magnetic moments 
which are large at infrared momenta; 
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• and gluons are nonpcrturbativcly massive, being 
described by a mass-function which is large in the 
infrared but diminishes with power-law behaviour 
in the ultraviolet. 

Numerous items could be added to this list, some of 
which arc described above. 

There are many reasons why this is an exciting time 
in hadron physics. These lecture notes emphasise one. 
Namely, through the DSEs, one is unifying phenomena 
as apparently diverse as: the hadron spectrum; hadron 
elastic and transition form factors, from small- to large- 
Q 2 ; parton distribution functions; the physics of hadrons 
containing one or more heavy quarks; and properties of 
the quark gluon plasma. The key is an understand- 
ing of both the fundamental origin of visible mass and 
the far-reaching consequences of the mechanism respon- 
sible; namely, DCSB. Through continuing feedback be- 
tween experiment and theory, these studies should lead us 
to an explanation of confinement, the phenomenon that 
makes nonperturbative QCD the most interesting piece 
of the Standard Model. They might also provide an un- 
derstanding of nonperturbative physics that enables the 
formulation of a realistic extension of that model. 
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